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PREFACE 

This  effort  was  conducted  by  The  University  of  Kentucky  under  the 
sponsorship  of  the  Rome  Air  Development  Center  Post-Doctoral  Program  for 
RADC's  Compatibility  Eianch.  Mr.  Jim  Brodock  of  RADC  was  the  task  project 
engineer  and  provided  overall  technical  dlrectiona  and  guidance. 

The  RADC  Post-Doctoral  Program  is  a cooperative  venture  between  RADC 
and  some  sixty-five  universities  eligible  to  participate  in  the  program, 
Syracuse  University  (Department  of  Electrical  Engineering) , Purdue  University 
(School  of  F-lectrical  Engineering),  Georgia  Institute  of  Technology  (School 
of  Electrical  Engineering) , and  State  University  of  New  York  at  Buffalo 
(Department  of  Electrical  Engineering)  act  as  prime  contractor  schools  with 
other  schools  participating  via  sub-contracts  with  the  prime  schools.  The 
U.S.  Air  Force  Academy  (Department  of  Electrical  Engineering),  Air  Force 
Institute  of  Technology  (Department  of  Electrical  Engineering),  and  the 
Naval  Post  Graduate  School  (Department  of  Electrical  Engineerings  also 
participate  in  the  program. 

The  Post-Doctoral  Program  provides  an  opportunity  for  faculty  at 
participating  universities  to  spend  up  to  one  year  full  time  on  exploratory 
development  and  problem-solving  efforts  with  the  post-doctorals  splitting 
their  time  between  the  customer  location  and  their  educational  institutions. 
The  program  is  totally  customer-funded  with  current  projects  being  undertaken 
for  Rome  Air  Development  Center  (RADC),  Space  and  Missile  Systems  Organi- 
zation (SAMSO),  Aeronautical  Systems  Division  (ASD) , Electronics  Systems 
Division  (ESD) , Air  Force  Avionics  Laboratory  (AFAL) , Foreign  Technology 
Division  (FTD) , Air  Force  Weapons  Laboratory  (AFWL) , Armament  Development 
and  Test  Center  (ADTC) , Air  Force  Communications  Service  (AFCS) , Aerospace 
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Defense  Command  (ADC).  Hq  USAF,  Defense  Communications  Agency  (DCA) , Navy, 
Army,  Aerospace  Medical  Division  (AMD),  and  Federal  Aviation  Administration 
(FAA). 

Further  information  about  the  RADC  Post-Doctoral  Program  can  je 
obtained  from  Mr.  Jacob  Scherer,  RADC/RBC,  Griffiss  AFB,  NY,  13443, 
telephone  Autovon  587-2543,  commercial  (315)330-2543. 

The  author  of  this  report  is  Clayton  R.  Paul.  He  received  the  USEE 
degree  from  Tho  Citadel  (1963) , the  MSFF  degree  from  Georgia  Institute  of 
Technology  (1964),  and  the  Ph.D.  degree  from  Purdue  University  (1970).  He 
is  currently  an  Associate  Professor  with  the  Department  of  Electrical 
Engineering,  University  of  Kentucky,  Lexington,  Kentucky  40506. 

The  author  wishes  to  acknowledge  the  capable  efforts  of  Ms.  Donna  To /n 
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The  problem  of  determining  the  currents  Induced  In  termination  networks 
at  the  ends  of  a multiconductor  transmission  line  by  an  Incident  electro- 
magnetic field  Is  obviously  quite  Important  In  determining  the  electromagnetic 
compatibility  of  electronic  systems.  The  digital  computer  program  described 
in  this  report  Is  Intended  to  be  used  for  this  purpose. 

The  specla]  case  of  a transmission  line  consisting  of  two  wires 
(cylindrical  conductors)  Immersed  In  a general,  nonuniform  field  was  consid- 
ered by  Taylor,  Satterwhlte  and  Harrison  [3].  The  equations  for  the  terminal 
currents  obtained  in  [3]  were  placed  In  a more  convenient  form  by  Smith  [4]. 

The  special  case  of  a uniform  plane  wave  incident  on  a three-wire  line 
(the  three  wires  lie  in  a plane)  In  the  transverse  direction  (perpendicular 
to  the  transmission  line  longitudinal  (x)  axis)  with  the  electric  field 
intensity  vector  polarized  parallel  to  the  line  axis  was  obtained  by 
Harrison  in  [5].  Paul  has  extended  these  special  case  results  to  (n+1) 
conductor  (multiconductor)  lines  for  an  arbitrary  Incident  electromagnetic 
field  [1,6]. 

This  report  describes  a digital  computer  program,  WIRE,  which  is 
designed  to  calculate  the  sinusoidal .steady  state .terminal  currents  induced 
at  the  ends  of  a uniform,  multiconductor  transmission  line  which  Is  illuminated 
by  an  incident  electromagnetic  field.  Three  types  of  transmission  line 
structures  are  considered.  TYPE  1 structures  consist  of  (n+1)  parallel 
wires.  TYPE  2 structures  consist  of  n wires  above  an  infinite  ground  plane. 
iYPE  3 structures  consist  of  n wires  within  an  overall , cylindrical  shield. 

For  each  structure  type,  one  of  the  conductors  is  designated  as  the 
reference  conductor  for  the  line  voltages.  For  TYPE  1 structures,  the 
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reference  conductor  is  one  of  the  (n+1)  wires.  For  TYPE  2 structures,  the 
reference  conductor  is  the  ground  plane.  For  TYPE  3 structures,  the  refer- 
ence conductor  is  the  overall,  cylindrical  shield. 

All  of  the  transmission  lines  are  considered  to  be  uniform  in  the  sense 
that  there  is  no  variation  in  the  cross-sections  of  the  (n+1)  conductors 
along  the  transmisi^ion  line  axis  and  all  (n+1)  conductors  >.re  parallel  to 
each  other.  All  conductors  are  considered  to  be  perfect  conductors  and  the 
surrounding  medium  is  considered  to  be  homogeneous,  linear,  isotropic  and 
lossless. 

The  incidert  field  can  be  in  the  form  of  a uniform  plane  wave  for  TYPE 
1 and  TYPE  2 structures  or  a nonuniform  field  for  all  structure  types.  The 
uniform  plane  wave  excitation  is  specified  by  data  entries  describing  the 
magnitude  of  the  electric  field  intensity  vector,  the  orientation  of  this 
vector  and  the  direction  of  propagation.  These  quantities  will  be  made  pre- 
cise in  the  following  chapters.  For  the  nonuniform  field,  the  data  entries 
are  the  values  of  ti'.e  incident  electric  field  intensity  (magnitude  and 
phase)  at  points  along  the  axes  of  the  conductors  and  along  contours  between 
the  wires  at  the  two  ends  of  the  line.  Piecewise-linear  behavior  of  the 
fields  (magnitude  and  pliase)  is  assumed  between  these  data  points. 

The  primary  restrictions  on  the  program  are  that  the  cross-sectional 
dimensions  of  the  line,  e.g.,  conductor  separations,  are  much  smaller  than 
a wavelength  at  the  frequency  in  question  and  the  ratios  of  conductor 
separation  to  wire  radii  are  greater  than  approximately  5.  The  first  re- 
striction is  imposed  to  insure  (in  a qualitative  fashion)  that  only  the  TEM 
mode  of  propagation  is  significant,  i.e.,  the  higher  order  modes  are  non- 
propagating. This  requirement  that  the  cross-sectior.al  dimensions  of  the 
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line  are  electrically  small  must  also  be  imposed  to  insure  that  the  definition 
of  voltage  is  independent  of  path  if  the  incident  field  is  not  curl  free  in 
the  line's  cross-sectional  plane.  (See  Chapter  II,)  The  second  restriction 
is  necessary  to  insure  the  validity  of  the  entries  in  the  ner-unit-length 
transmission  line  inductance  and  capacitance  matrices.  The  entries  in 
these  matrices  are  derived  by  assuming  that  the  per-unit-length  charge  dis- 
tributions on  the  wires  are  essentially  constant  around  the  peripheries  of 
the  wires,  i.e,,  the  wires  are  separated  from  each  other  sufficiently  to 
insure  that  proximity  effect  is  not  a factor. 

General  termination  structures  are  provided  for  at  the  ends  of  the 
transmission  line.  These  terminations  are  assumed  to  be  linear. 

Chapter  II  contains  the  derivation  of  the  equations  for  general  field 
excitations.  Chapter  III  contains  a derivation  of  the  equivalent  sources 
induced  in  the  structure  types  by  uniform  plane  waves  as  well  as  jionuniform 
fields.  Chapter  IV  contains  a discussion  of  the  contents  of  the  program. 
Chapter  V contains  a User's  Manual  and  Chapter  VI  contains  examples  which 
are  used  to  check  the  program  operation. 


II.  MODEL  DERIVATIONS 


Cioss-sections  of  the  three  basic  types  of  structures  considered  by  Che 
program  are  shown  in  Figure  2-1.  The  axis  of  the  line  is  the  x coordinate 
and  the  (n+1)  conductors  are  perpendicular  to  the  y,z  plane  as  indicated  in 
Figure  2-1.  The  TYPE  1 structure  consists  of  (n+1)  wires  in  which  one  of 
the  wires  is  designated  as  the  reference  conductor  for  the  line  voltages. 

The  TYPE  2 structure  consists  of  n wires  above  an  infinite  ground  plane 
where  the  ground  plane  is  the  reference  conductor  for  the  line  voltages. 

The  TYPE  3 structure  consists  of  n wires  within  an  overall  cylindrical 
shield.  In  this  case,  the  shield  is  the  reference  conductor. 

All  conductors  are  considered  to  be  perfect  conductors  and  the  sur- 
rounding medium  is  considered  to  be  homogeneous,  linear,  isotropic  and 
lossless.  The  surrounding  medium  (homogeneous)  is  characterized  by  a 


permittivity  e and  a permeability  p.  Throughout  this  report,  the  perme- 
ability and  permittivity  of  free  space  will  be  denoted  by  = 4it  x 10  ^ 
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and  s (l/36n)  x 10  , respectively, and  the  permeability  and  permittivity 

of  the  medium  are  related  to  the  free  space  values  by  the  relative  perme- 
ability, and  relative  permittivity  (relative  dielectric  constant),  e^, 
as  p = p^  p^  and  e = e^,  respectively.  For  structure  TYPE  1 and  TYPE  2, 

a logical  choice  for  and  p^  would  be  1 (free  space).  For  structure  TYPE  3, 

a logical  choice  for  the  relative  permeability,  p^,  would  be  1 as  is  typical 
of  dielectrics.  The  program,  however,  allows  for  any  and  p^  for  all 
structure  types. 

The  n wires  are  labeled  from  1 to  n and  the  radius  of  the  i-th  wire  is 
denoted  by  r ^ , The  reference  conductor  is  designated  as  the  zero-th  con- 

Wx 

ductor.  For  TYPE  i structures,  the  reference  wire  has  radius  r _ and  the 
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center-to-center  separation  between  the  i-th  and  j-tli  wires  is  designated  as 
djj . For  TYPE  2 structures,  the  i-th  wire  is  at  a height  h^  about  the  ground 
plane  with  a center-to-center  separation  between  the  i-th  and  j-th  wires  of 
d^^.  For  TYPE  3 structures,  the  interior  radius  of  the  cylindrical  shield 
is  designated  by  r^,  the  i-th  wire  is  at  a distance  r^  from  the  shield  center 
and  the  angular  separation  between  the  i-th  and  j-th  wires  is  designated  by 


Implicit  in  the  following  is  the  requirement  for  the  transmission  line 
to  be  uniform.  Transmission  lines  considered  here  are  uniform  in  the  sense 
that  all  (n+1)  conductors  have  uniform  cross-sections  along  the  line  axis 
and  all  n wires  are  parallel  to  each  other  and  the  reference  conductor. 

2,1  Derivation  of  the  Multiconductor  Transmission  Line  Equations 

The  distributed  parameter  transmission  line  equations  for  multiconductor 
lines  with  incident  field  illumination  can  be  derived  and  are  similar  (with 
matrix  notation  employed)  to  the  familiar  equations  for  two-conductor  lines 
[1,2, 6, 7, 8],  Assuming  sinusoidal  excitation  at  a radian  frequency  u)=2TTf, 
the  electric  field  intensity  vector,^ (x,y,z,t),  and  the  magnetic  field 
intensity  vector,  (x,y,  z,  t) , are  written  as  ^{x,y,z,t)  = E(x,y,z)e^“‘^  and 
•^(x,y,z,t)  = H(x,y,z)e^*^^.  The  complex  vectors  E(x,y,z)  and  H(x,y,z)  are 
the  phasor  quantities.  Line  voltages,2/ j^(x,t)  = V^(x)e^*^*',  of  the  i-th 
conductor  with  respect  to  the  zeroth  conductor  (the  reference  conductor)  are 
defined  as  the  line  integral  of^  between  the  two  conductors  along  a path  in 
the  y,z  plane.  V^(x)  is  the  complex  phasor  voltage.  The  line  current  ,xSI^(x,t)>= 
I^(x)e^*^*'  associated  with  the  i-th  conductor  and  directed  in  the  x direction 
is  defined  as  the  line  integral  of^  along  a closed  contour  in  the  y,z  plant 
encircling  only  th e i-th  conductor  and  Ij(x)  is  the  complex  phasor  current. 


's: ' 

II 
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Ttie  current  in  the  reference  conductor, ^^(Xjt)  = satisfies  1^  = 

It  is  convenient  to  consider  the  effects  of  the  spectral  components  of 
the  incident  field  as  per-unit-length  distributed  sources  along  the  line. 
The  sources  appear  as  series  voltage  sources  and  shunt  current  sources  as 
indicated  in  Figure  2-2  for  an  "electrically  small"  Ax  section  of  the  line. 
The  multiconductor  transmission  line  equations  may  then  be  derived  for  the 
Ax  subsection  in  Figure  2-2  in  the  limit  as  Ax-H)  as  a set  of  2n  coupled, 
complex,  ordinary  differential  equations  [1], 

V(x)  + jwLI(x)  = Vg(x)  (2-la) 

i(x)  + ja3Cy(x)  = Ig(x)  (2-lb) 


A matrix  M with  m rows  and  n columns  is  denoted  as  mxn  and  the  element 
in  the  i-th  row  and  j-th  column  is  denoted  by  V(x)  and  I(x)  are 

nxl  vectors  of  the  line  voltages  and  currents,  respectively.  The  elements  In 
the  i-th  rows  are  [y(x)]^  =■  V^(x)  and  [I(x)]^  = [V(x)]^  = 

(d/dx)V. (x).  The  nxn  real,  symmetric,  constant  matri  es  L and  C are  the 
per-unit-length  inductance  and  capacitance  matrices,  respectively.  From 
Figure  2-2  one  can  derive  (2-1)  and  the  entries  in  L and  C become  [1] 


and 


^I^^ll  “ ^i  ^0  " ^"“iO 

(2-2a) 

’ ^0  ""ij  ~ ™i0  ■ "jO 

(2-2b) 

i^5 

■ 'iO  'ij 

j=l 

(2-3a) 

[Sjij  ■ -'ij 

(2-3b) 

1 1*  J 
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The  entries  in  V (x)  and  I (x)  are  the  per-unit-length  distributed  sources 
along  the  line  induced  by  the  incident  field,  i.e.,  [V^Cx)]  = Vg^(x)  and 
ll^(x)]  = Ig^(x),  as  shown  in  Figure  2-2. 

In  order  to  consider  general  termination  networks  (and  allowing  inde- 
pendent sources  in  these  networks)  we  may  characterize  these  as  generalized 
Thevenin  equivalents  [1].  For  a line  of  total  length^,  the  equations  for 
the  termination  networks  at  x = 0 and  x = ^ are 

y(0)  = Vq  - ZqI(O)  (2-4a) 

y(^)  = I(/)  (2-4b) 

where  and  are  nxl  vectors  of  equivalent  open  circuit  port  excitation 
voltages,  [Vg]^  = and  and  Zg  and  Z^are  nyn  symmetric  im- 

pedance matrices  as  shown  in  Figure  2-3.  This  is,  of  course,  a completely 
general  and  arbitrary  characterization  of  these  linear  termination  networks. 
The  entries  in  these  termination  equations  can  be  easily  determined  for  a 
given  network  by  considering  V^(0)  and  V^(i)  (the  termination  port  voltages) 
as  independent  sources , and  writing  the  loop  current  equations  for  each 
network  where  ij^(O)  ^^d  I^(3L)  are  subsets  of  the  loop  currents  in  each 
network.  (See  Section  2,6.) 

With  the  line  immersed  in  a homogeneous  medium  with  permittivity 
e and  permeability  u,  the  product  of  L and  C becomes  [1] 

LC  » CL  *=  pel  (2-5) 

,.n 

where  1^  is  'the  nxn  identity  matrix  with  ones  on  the  main  diagonal  and  zeros 

elsewhere,  i.e.,  [T  ]..  = 1,  and  [1  ] . . = 0,  iii^j . Fcr  this  case,  the 
’ ’ \.n  ii  ’ ~n  ij 

solution  to  (2-1)  and  (2-4)  is  in  a simple  form  [1] 


I 


[cos(kX)/ZQ  + Z^l  + j sln(k/)  j Z^  + Z^Z^h^]]  1(0) 
= -V^+  Ij  sin(k^)  Z^Z~^  + cos(kjf)l^]  Vq 


l(i)  = -j  sln(Ki:)2^"\Q 
+ [cos(kZ)l^  + j sln(k;^)Zp"^ZQ]  1(0)  + IgC?) 


(2-6a) 


(2-6b) 


where  the  wavenumber  is  k = 2Tr/A,  A ® v/f,  v = 1//^  = v^Z/ji^e^,  Vq  = 1/ 
and  the  nxn  characteristic  impedance  matrix,  Zq,  is  [1] 


Zc  =*  V L 


(2-7) 


A 

The  inverse  of  an  nxn  matrix  M is  denoted  by  M and  V (x)  and  I (X)  in 

«w  «•  *S  ““S 

(2-6)  are  given  by  [I] 


(2-8a) 


(2-8b) 


V (s()  * j {cos(k(>^  - x))  V (x) 
-s  Jq  -s 


-j  sin(k(;f  - x))  ZqI^(x)}  dx 


I (jf)  “ j {cos(k(/-  x))  I (x) 
-s  Q -8 


-j  sin(k(/  - x))  Zp  (x)}  dx. 

^ *"S 


Solution  of  (2-6a)  for  the  current  vector,  1(0),  requires  the  solution  of  n 
complex  equations  in  n unknowns  (1^(0)),  Once  (2-6a)  is  solved,  (2-6b) 
yields  the  currents  I(jf)  directly. 

In  this  report,  no  Independent  excitation  sources  in  the  termination 
networks  will  be  considered.  The  program  XTALK  described  in  Vol,  VII  of 
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this  series  [2]  can  be  used  to  compute  the  contribution  to  the  response  due 

to  these  sources.  Thus  the  source  vectors  in  the  generalized  Thevenin 

equivalent  representations  in  (2-4)  will  be  zero,  i.e.,  = V / = 0 where 

-u  -d- 

the  r.i)(p  zero  matrix,  0 , has  zeros  in  every  position,  i.e.,  [ 0 1 =0  for 

m-p  m^p^ 

i = 1;  ...,  m and  j=l,  ...,  p.  Thus  the  generalized  Thevenin  equivalent 

representation  becomes 


y(0)  = -Zq  1(0)  (2-9a) 

V(t)  = Z^I(;^)  (2-9b) 


end  the  equations  for  the  terminal  currents  in  (2-6)  become 

fcos(kj0  {Zq  + Z^}  + j sin(k;t)  IZ^  + Z^"^  Z^}  ] 1(0)  = 

(2-lOa) 

UtO  = (cos(ki)l^  + j sin(kX)  Z^"^  Zq]  1(0)  + 1^(0  (2-lOb) 


As  an  alrernat^,'  formulation,  a generalized  Norton  equivalent  represent- 
ation may  be  usei  to  (haracterize  the  termination  networks.  It  we  define 
Y-  = Z_  ^ and  Yy  =•  7. , ^ ‘•he  generalized  Norton  equivalent  representation 
becomes 


I'O)  = -Yq  y(0) 
U/)  - Y^  V(/) 

Equations  (2-10)  can  then  be  written  as 


(2-lla) 

(2-llb) 


[cos(k?:)  {Yq  + Y^}+  j sin(kr)  'i^Z^  Yp  ’■  Z^"^)  ] V(0)  = 
I3(Z)  - Y^y,(zr) 


(2-12a) 


lU) 


-[cos(k:f)  Yq  + j sin(ki) 


Z(,”^]  y(0)  + l^(rt  (2-12b) 


where  1(0)  can  be  recovered  from  V(0)  via  (2-lla). 

There  remain  two  basic  problems:  determining  the  entries  in  the  per- 
unit-length  inductance  and  capacitance  matrices,  L and  C,and  determining 

A A 

the  equivalent  source  vectors,  V^(;C)  and  which  are  induced  by  the 

incident  electromagnetic  field.  The  derivations  of  L and  C for  the  three 
Structure  types  have  been  given  previously  [1,9]  and  will  be  summarized  in 
the  following  sections.  It  will  become  clear  in  the  following  section  that 

A A 

once  the  equivalent  source  vectors,  and  are  determined  for 

the  TYPE  1 structure,  they  can  be  immediately  obtained  for  the  TYPE  2 and 
TYPE  3 structures  with  a parallel  development.  Thus  the  basic  problem  is 
the  determination  of  these  equivalent  source  vectors  for  the  TYPE  1 structure. 


2 , 2 Derivation  of  the  Equivalent  Induced  Source  Vectors,  V. (^)  and 


for  TYPE  1 Structures 

In  order  to  determine  the  equivalent  induced  sources,  and  I^^Cx), 

consider  Figure  2-4,  The  method  used  in  [3]  can  be  adapted  here  in  a 
similar  fashion.  Faraday’s  law  in  integral  form  becomes 


/ E . dCj  • -Jw  j H • ; dSj 


(2-13) 


where  is  a flat,  rectangular  surface  in  the  x,y  plane  between  wire  i and 
wire  0 and  between  x and  x + Ax  as  shown  in  Figure  2-4.  The  unit  normal 
ii  is  n ” z where  z is  the  unit  vector  in  the  z direction,  dS^  = dx  dy  and 
is  a contour  encircling  in  the  proper  direction  (counter-clockwise 
according  to  the  right-hand  rule).  Equation  (2-13)  becomes  for  the  indicated 
integration^ 

■^In  integrating  from  y=0  to  implicitly  assuming  that  the  wires 

are  sufficiently  separated  so  that  they  may  be  replaced  by  infinitesimally 
small  filaments  of  current  (charge). 
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[E^j^(y,x  + Ax)  - E^^(x,y)]  dy 


^^111^*^10’^^  ” E^^(0,x)]  dx 


r x+Ax  f d . 

■ -'“i  J, 


H^i(y,x)  dy  dx 


(2-14) 


where  E Is  the  component  of  the  total  electric  field  (incident  plus 

Li 

scattered)  transverse  to  the  line  axis  and  lying  along  a straight  line 
joining  the  two  conductors  i.e.,  = E^;  is  the  component  of  the  total 

electric  field  along  the  longitudinal  axis  of  the  line,  i.e,,  E , ■=  E ; 

J.  X 

and  H . is  the  component  of  the  total  magnetic  field  perpendicular  to  the 
nl 

plane  formed  by  the  two  wires,  i.e,, 

Defining  the  voltage  between  the  two  wires  as 


d 

i^U)  - -j 


E^^(y,x)  dy 


(2-15) 


dV^(x) 


lim  J 
Ax+0 


1 

nr  J 


:^j^(y»x  + Ax)  - 


Ej.^(y,x)]  dy 


(2-16) 


The  total  electric  field  along  the  wire  surfaces  is  zero  since  we  assume 
perfect  conductors.  (One  can  straightforwardly  include  finite  conductivity 
conductors  through  a surface  impedance  as  was  done  in  [3]).  Therefore 
(2-14)  becomes  in  the  limit  as  Ax  ->■  0 


I 


H 


i 


i 


dV. (x) 
1 


dx 


. f iO 

= jwp 


J 

0 


Hj^i(y,x)  dy. 


(2-17) 


Tlie  total  magnetic  field  is  the  sum  of  an  incident  and  a scattered  field 


H^.  (y.x)  = H^(y,x) 


(scat)  (inc) 


(2-18) 


= H (y,x)+  H (/,x) 

^ Z 


scattered  Incident 

and  the  scattered  field  here  is  considered  to  be  produced  by  the  transmission 
line  currents.  The  scattered  flux  passing  between  the  two  conductors  per 
unit  of  line  length  is  directly  related  to  the  scattered  magnetic  field  and 
the  per-unit-length  inductance  matrix,  L,  as 


(scat)  f 

>^(x)  - - J 


(scat)  f “iO  (scat) 
0 


^^il’^i2’”*’^in^ 


Ij^(x) 


t2(x) 


In(x) 


(2-19) 


where  Substituting  (2-19)  and  (2-18)  into  (2-17)  and  arranging 


for  i = 1,  ..o,  n yields 


V(x)  + jmLI(x) 


/iO  (inc) 

H(y,x)  dy 

0 : 


(2-20) 
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I 


and  the  source  vector  V^(x)  In  (2-1)  is  easily  identified  by  comparing 
(2-20)  and  (2-1). 

For  transmission  line  theory  to  apoly,  the  cross-sectional  dimensions 

of  the  line  (wire  spacing^  etc.)  must  be  electrically  small,  i.e.,  <<  1. 

Thus  the  result  indicates  that  the  voltage,  V , induced  in  the  loop 

si 

between  the  Jth  conductor  and  the  zeroth  conductor  and  between  x and 
X + Ax  is  equal  to  the  rate  of  change  of  the  incident  flux  penetrating  this 
"electrically  small"  loop  which,  of  course,  makes  sense. 


Ampere's  law  yields 


, 3H  3H 

g _ 1 X _ z 

y jme  8z  3x 


(2-21) 


Ey  will  consist  of  scattered  and  incident  field  components  and  is  written 


E (y*x) 


(scat)  (inc) 

- E (y,x)  + E (y,x). 

scattered  incident 


(2-22) 


Substituting  (2-21)  into  (2-15)  we  have 


'l(x)  - - J 


Ey(y.x)  dy 


---  f 

jwe  J 
(sea?) 


( (scat)  (inc) 

j f iO  j 3H^(y,x)  9H^(y,x) 


(inc) 


3H^(y,x^  3H^(y,x) 


(2-23) 


Utilizing  (2-19)  we  obtain 


V^(x) 


jilie  L 


--^  f 

3we  j 


(scat) 

dio  3H^(y.x) 


3z 


dy 


/ 


^io 

^^E^i(y,x)  dy. 


(2-24) 


If  we  assume  that  the  currents  on  the  wires  are  directed  only  in  the  x 


direction,  i.e.,  tiiere  are  no  transverse  components  of  the  currents  on 

(scat) 

the  wire  surfaces  , then  H^(y,x)  = 0 and  (2-24)  becomes 


1 


I 


m 


i 


I (x)  = - joiC 


/ 


iO 


(inc)  I 
E^^(y,x)  dy 


(2-27b) 


The  shunt  current  sources  in  I (x)  are  therefore  a result  of  the  line 

-s 

voltage  induced  by  the  incident  electric  field  being  applied  across  the 
per-unit-length  line-to-line  capacitances  which,  of  course,  satisfies  our 
intuition. 

A 

The  final  problem  remaining  is  to  obtain  simplified  versions  of  V^and 

A 

in  (2-8)  to  be  directly  used  in  (2-10)  and  (2-12).  First  consider  the 

A 

determination  of  (i) , Substituting  (2-27)  into  (2-8a)  yields 


S!.«) 


COS  (k(*f  - x)) 


I 


I' 


I 


X[/  c 


iO  (inc)  t 

(y,x)  dy 


J dX 


- k 


sin  (k(;Z'-  x)) 


X 


iO 


(inc)  I 

E . (y,x)  dy 
tl  • 


? dx. 


From  Faraday’s  law  we  obtain 


(inc) 

«ni 


1 

jwp 


(inc) 

9y 
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(inc) 

5E 

—Ci. 

9x 


(2-28) 


(2-29) 


I 


Substituting  this  into  (2-28)  yields 


Y.,(iC) 


■I'f 

■(■■ 


(inc)  I 

os(k(^-  x)) 


(inc)  j 
E^i(0,x)  > dx 


cos  (k(jf  - x)  ) 


^ (inc): 

/[  ,-■‘10  aEti(y.^) 

X i 


- k I < sin  (k(<^-  x)) 


dy  r dx 


(2-30) 


XI 


•dio  (inc)  : I 


Utilizing  Leibnitz's  rule  (see  [10,  p.  219]),  (2-30)  is  equivalent  to 


I = ^ * cos  (k(;^  - x)) 

• T 

(inc)  : (inc) 

I ’ 

- I <cos(k(3:-  x)) 

J ax 


cos(k(3:-  x)) 


(2-31) 


E^j^(y,x)  dy  dx 


m 


I 

7 

I 


P 


i' 


I 

I 


I 


and  this  may  be  written  as 


iu)  -f 


[ 


cos  (k(3!  - x) ) 


X 


(inc)  ; (inc) 
^)li^‘^iO’*^  " Ej^^(0,x) 


dx 


O'^iO  (inc)  : 

1^  Eti(^,X)  dy 


+ cos(k3!) 


d (inc) 

^ E^^(y,0) 


dy. 


Similarly  1 (jC)  may  be  obtained  as 


.if, 

IsW  = -iZc  Jo  <sin  (k(/-  X)) 


X 


(inc)  ’ (inc) 


r dx 


- iz 


-u 


sin  (k^) 


f *^10 

J E^^(y,0)  dy 
- *0 


(2-32) 


(2-33) 


The  important  quantity  in  (2-lOa)  is  V (X)  ' Z^l  (jC).  Combining  (2-32)  and 

*^S  *,x**s 

(2-33) , this  becomes 


■I 


[cos  (k(;^-  x))  + j sin  (k(/ - x)) 


X 


(inc)  I (inc) 


dx 


/ 
•-  0 


10  (inc) 

E^i(y.3f) 


dy 


(2-34) 


+ [cos(lof)  1 + j sin  (ki)Z..Z  "■■■] 
^n  ,0 


(inc) 

E^i(y,0) 


dy. 


Note  that  the  equivalent  forcing  function  on  the  right-hand  side  of  (2- 10a), 

A A 

i given  in  (2-34)  is  simply  determined  as  a convolution  of 

differences  of  the  incident  electric  field  vector  along  the  wire  axes, 

(inc)  (inc) 

E^i(dio,x)  - E^^(0,x),  and  a linear  combination  of  integrals  of  components 

of  the  electric  field  vectors  at  the  endpoints  of  the  line  which  are  trans- 
(inc)  ’nc) 

verse  to  the  line,  F^^(y,3C)  and  E^^(y,0),  This  is,  of  course,  precisely 
the  result  obtained  by  Smith  [4]  for  two  conductor  lines.  Substituting 
(2-34)  into  (2-lOa)  one  can  verify  that  the  result  reduces  for  two  con- 
ductor lines  (n  = 1)  to  the  result  given  by  Smith  [4]  since  Z^,Zjj^,Zq  become 
scalars  for  two  conductor  lines  and  (2-lOa)  becomes  one  equation  in  only 
one  unknown  1(0). 
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The  final  equations  for  the  line  currents  then  become  (substituting 
(2-34)  into  (2-10) 

[cos(kX^)  ^Zo  + Zjj}  + j sin  (k;^)  {Z  + Z-Z  "^Z^}]  1(0)  = 

- ~C  -<*-C  “ 

/ ^fcos(k(2'-  x))  1^  (2-35a) 

0 

^ _ -1.  (inc) 


+ j sin  (k(Z'-  x))  ^Z^  E^(x)}  dx  - E^(£) 

1 (inc) 

+ {[cos(k;r)l^  + j sin(k;r)^Zj,~-^]  £^.(0)} 


I(i)  » [cos(k:0l^  + j sin(kt’)^“^ZQ]  1(0) 


-1  r*f  ^ 

1 {sln(k(2'-  x))  (x)}  dx 


^3:  f 

•'O 


(2-35b) 


(inc) 

- ,1Z^  {sln(k5;)  Ej.(0)} 

(Inc)  (inc)  (inc) 

where  Ej,(x),  E^(jf),  and  E^(0)  are  nxl  column  vectors  with  the  entries  in 
the  i-th  rows  given  by 

(inc)  (inc)  (inc) 


= 1 

(2-36a) 

(inc) 

r'^iO  (inc) 

- J 

Iq  ^ti^Pi*^^  ‘^i 

(2-36b) 

(inc)  j 

r djQ  (inc) 

E^l(Di.O)  dp^ 

0 
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(2-36c) 

i 


for  i = 1,  . . . , n. 

A word  of  caution  in  the  Interpretation  of  the  notation  is  in  order. 

Although  it  should  be  clear  from  the  derivation,  the  reader  should  never- 

(inc) 

theless  be  reminded  that  the  integration  path  for  the  component  is 

in  the  y direction  when  the  1-th  conductor  is  concerned.  When  other 
conductors  are  concerned,  the  integration  path  is  a straight  line  in  the 
y,2  plane  which  joins  the  conductor  and  the  zeroth  conductor  and  is  per- 
pendicular to  these  two  conductors.  This  is  designated  as  in  (2-36) 
and  replaces  the  y variable  for  the  path  associated  with  conductors  i and 
0,  The  notation  may  be  cumbersome  but  the  idea  and  the  implementation  are 
quite  simple. 

Defining  the  vectors 


(inc) 


M 


= J cos  x))  (x)  dx 


(2-37a) 


(inc) 


IT  ^ 

N = sin  (k(^  - x))  E.  (x)  dx 

" Jo 


we  may  write  (2-35)  as 


(2-37b) 


Icos(W)  {Zq  + + j sin(W)  {Z^  + ZjfZ  }]  1(0) 


-1. 


(inc) 


= M + j "N  - £^.(•5^) 


+ [cos(k3:)  1 + j sin(kaf)  Z.Z„  ] 

-.n  -t 


(inc) 

(0) 


-1. 


I(sf)  = [cos(k:t)  1 + j sln(k:^)  Z„  Zj  1(0) 
*"  .„u  *• 

, , (Inc) 

- j ^ N - j Z^,  ^ {sln(k/)  E^(0)} 


(2-38a) 


(2-38b) 
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For  the  generalized  Norton  equivalent  representation,  equations 
(2-38)  can  be  written  as 

[cos(ka:){YQ  + Y + j sin(k;f)  {Y^Z^Y^  + Z'^}][-V(0)]  = 

Y^M  + j Z"^  N - Y^  (Z)] 

+ [cos(ki)  Y^+  j sin(k<)  Z"^][E^^"''\0)] 


(2-39a) 


l(.-t)  = [cos(kZ)  Yq  + j sin(k<)  Z^^][-V(0)] 
-j  Z“^  N - j Z"^  {sin(k?) 


(2- 39b) 


and  1(0)  is  obtained  from  1(0)  = -Y^  V(0)  = Yq[-V(0)], 

2 . 3 Determining  the  Per-Unit-Length  Inductance  Matrix,  L,  for  TYPE  1 
Structures 

For  TYPE  1 structures,  one  final  calculation  remains;  the  determination 
of  the  per-unit-length  inductance  matrix,  L,  which  is  related  to  the 
characteristic  impedance,  Z^,  via  (2-7),  Ordinarily  this  is  a difficult 
calculation  [11],  However,  if  we  assume  that  the  wires  are  separated  suf- 
ficiently such  that  the  charge  distribution  around  the  periphery  of  each 
wire  is  constant,  then  the  wires  can  be  replaced  by  filamentary  lines  of 
charge.  Typically,  this  will  be  accurate  if  the  smallest  ratio  of  wire 
separation  to  wire  radius  is  greater  than  approximately  5 [11],  In  this 
case,  the  entries  in  L for  TYPE  1 structures  are  given  by  [1,9] 

wi  wO 
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(2-40b) 


‘^<11  - >'=  >?llj  ■ 2^  ‘”  <^> 


wO  Ij 


For  closer  wire  spacings,  proximity  effect  will  alter  the  charge  distribution 
from  constant  ones  and  numerical  approximations  must  be  employed  to  find 
C and  L [11].  Although  the  entries  in  L have  been  derived  elsewhere,  we 
shall  show  a direct  derivation  which  relates  the  scattered  flux  passing 
between  the  wires  to  the  wire  currents  as  was  used  in  (2-19). 

The  matrix  L relates  the  scattered  flux  $ passing  between  the 
wires  to  the  wire  currents  as 


(scat) 




^1 

« 

• 

(scat) 

s 

L^n  J 

11 


In 


nn 


(2-41) 


The  respective  entries  are  determined  as 

(scat) 


ii 


T • ♦ • I T 

^1’  » i-1’^ 


i+T 


.ln=0 


(2-42a) 


ij 

i»‘j 


(scat) 

A 

^j 


^1’ 


j+i’ 


(2-42b) 


and  Large  wire  separations  are  assumed  so  that  the  wires  may 

be  replaced  by  filaments  of  current.  VJhen  the  wires  are  not  widely 
separated,  accurate  values  for  L can  be  obtained  by  numerical  methods  [11]. 
Consider  Figure  2-5 (a).  The  magnitude  of  the  magnetic  field  intensity 
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vector  due  to  I.  on  wire  i at  a distance  r>r  . away  from  wire  i is 
1 wi  ^ 


H = 


I. 

1 


r 2TTr 


(2-43) 


and  the  total  flux  passing  between  wire  i and  wire  0 due  to  both  currents 
is 


(scat) 


P 

Itt 


I r r 

(I  i dr  -f  J 1 

TT  J r Jr 


iO 


dr  } 


wi 


wO 


(2-44) 


fcn  ( - — - — ) . 


2n 


Thus  is  easily  identified  as  in  (2-40a). 


(scat) 


Consider  Figure  2-5 (b).  The  portion  of  the  flux  passing  between 


wire  i and  wire  0 due  to  -1^  in  the  reference  conductor  is  as  above 


(scat)  pl^  d^Q 

*^10  ""  2x 


£n  (^) 


(2-45) 


wO 


and  the  portion  of  the  flux  passing  between  wire  i and  wire  0 due  to  1^  in 


the  jth  conductor  can  be  found  to  be 


(scat) 

'll 


f 

= -P  J 
0 

2 IT  3 J 


dp 
p=d 


(2-46) 


iO 


(p  - Pj) 


p=0 


l«0^  + (»  - Po>'l 


dp  ) 


■ ^ I.  (i  ' 
2 TI  J 2 


^ Po'> 


r] 


2 2 M 

^'^iO  - 
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Combining  (2-45)  and  (2-46)  we  obtain 


(scat)  (scat)  (scat) 


♦lo  * 


ah.  i„  (fialio , 

V»o 


since 


^ <‘^10  - 


J 2 _ 2 ^ 2 

‘^jO  “ ^0  ^0 


(2-47) 


(2-48a) 


(2-48b) 


and  Is  easily  Identified  as  In  (2-40b). 

2.4  Determination  of  the  Equivalent  Induced  Source  Vectors  and  the  Per- 
Unlt-Length  Inductance  Matrix  for  TYPE  2 Structures 

Consider  the  system  of  n wires  above  an  Infinite  ground  plane  shown 
In  Figure  2-1  (b).  The  result  for  (n+1)  wires  given  In  (2-35)  - (2-39) 
can  be  extended  to  this  case  with  the  following  observations.  Consider 
Figure  2-6.  Clearly  we  may  apply  Faraday's  law  in  the  previous  development 
to  the  flat,  rectangular  surface  in  the  x,y  plane  shown  in  Figure  2-6(b) 
between  the  ground  plane  and  the  1-th  wire  and  between  x and  x-t^x.  This 
flat,  rectangular  surface  lies  In  the  x,y  plane.  Equations  (2-35)  - 
(2-39)  will  again  be  obtained.  Equations  (2-36)  become  for  this  case 

(inc)  (inc)  (inc) 

[Ej^(x)lj^  » " E4i(0,x)  (2-49a) 

(inc)  /’i  (inr'* 

lEt(X)li  = j E^i(p^,3t)  dp^  (2-49b) 

''c 
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(2-49c) 


(inc) 


(inc) 


[E^(0)li  -J  Xj(p,.0) 


where  is  a straight-line  cont'^ur  in  the  y,z  plane  between  the  position 

of  the  ground  plane,  y=0,  and  the  i-th  wire,  and  is  perpendicular  to  the 

ground  plane,  i.e.  p . = y.  This  is  indicated  in  Figure  2-6 (a). 

(inc) 

Epj  (h  x)  is  the  component  of  the  incident  electric  field  parallel 
^ (inc) 

to  the  axis  of  the  i-th  wire  at  y=h^  and  (0,:c)  is  the  component  of 
the  incident  field  parallel  to  the  ground  plane  directly  beneath  the  i-th 
wire.  In  the  program,  it  assumed  that  the  net  incident  electric  field 
(the  vector  sum  of  the  incident  field  in  the  absence  of  the  ground  planf? 


and  the  portion  of  this  field  which  is  reflected  by  the  ground  plane)  is 

(inc) 

obtained.  Therefore  F (0,x)  = 0. 

X<!l  Cl 


is  the  component  of  the 


incident  electric  field  parallel  to  p^  and  directed  in  the  +y  direction. 

The  per-unit-length  inductance  matrix,  L,  can  be  obtained  in  a fashion 
similar  to  Section  2.3  by  determining  the  scattered  magnetic  flux  passing 
through  the  surface  between  the  1-th  wire  and  the  position  of  the 
ground  plane  (the  ground  plane  is  replaced  by  image  wires)  and  is  given 
by  fl.9] 

2h, 


U 

‘wi 


(2-50a) 


f^^ij  = IT 


d * 
ij 


(2-50b) 


for  i,  j=l,  , n where 


d.  .*  .^  + 4h,h. 

ij  ' ij  i j 


(2-51) 
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2 . 5 Determination  of  the  Equivalent  Induced  Source  Vectors  and  the 


Per-Unit-Length  Inductance  Matrix  for  TYPE  3 Structures 


Consider  the  system  of  n vjires  within  an  overall,  cylindrical  shield 

shown  in  Figure  2-1 (c).  Obviously,  a parallel  development  to  that  of 

Section  2-2  and  2-A  can  be  used  to  obtain  the  equations  (2-35)  - (2-39). 

The  image  of  the  i-th  wire  (assuming  the  i-th  wire  can  be  replaced  by 

2 

a filament)  is  located  at  a distance  of  r , from  the  shield  center  as 

s / 

shown  in  Figure  2-7  [1],  Equations  (2-36)  become  for  this  case 


(inc)  (inc)  (inc) 

[Ej^(x)]^  = - E^^(0,x)  (2-52a) 

(inc)  ^s  ’^i  (inc) 

[Et(:^)]i  = / ^2-52b) 


(inc)  *^1  (inc) 

[Ej.(0)]^  “ / Eti(p^»0)  dp^  (2-52C) 

0 


where  is  a straight-line  contour  in  the  y,z  plane  between  the  i-th  wire 

and  its  image  and  beginning  at  the  Interior  of  the  shield.  This  contour 

is  on  a line  between  the  1-th  wire  and  its  image. 

(inc) 

F^.(r  - r ,x)  is  the  component  of  the  incident  electric  field  parallel 

^ (inc) 

to  the  axis  of  the  i-th  wire  at  y=r^-r^  and  E^^(0,x)  is  the  component 

of  the  incident  field  parallel  to  the  axis  of  the  shield  on  the  interior 

of  the  surface.  In  tiie  program,  it  is  assumed  that  the  net  Incident  electric 

(inc)  (inc) 

field  is  obtained  so  that  E (0,x)  • 0.  E is  the  component  of  the 

X- 1 w X 

incident  electric  field  parallel  to  and  directed  in  the  +y  direction. 

The  per-unlt-length  inductance  matrix,  L,  can  be  obtained  in  a fashion 
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similar  to  Section  2,3  by  determining  the  scattered  magnetic  flux  passing 
through  a surface  between  the  1-th  wire  and  Interior  of  the  shield  in  the 


y,x  plane  and  is  given  by  [1,9] 


t!;^ii  = 27 


2 2 
r - r . 

tn  i-) 

^s  ’^wi 


(2-53a) 


•*  s 


(r.r.)^  + r ^ - 2r.r.r  cos  0 , . 
i j s 1 1 s ii 


(rir^)2  + cos 


(2-53b) 


where  0^^  is  the  angular  separation  between  the  i-th  and  j-th  wires  (see 
Figure  2-1). 


2,6  Determining  the  Entries  in  the  Termination  Network  Impedance 


(Admittance)  Matrices 

In  order  to  implement  this  method,  one  is  required  to  determine  the 
entries  in  the  nxn  terminal  Impedance  (admittance)  matrices,  Zq  and 
(Yq  and  Y^),  which  characterize  the  termination  networks  at  the  two  ends 
of  the  line  as: 


V(0)  = -Zq  1(0) 
V(r)  = Z;f  I(i) 


Thevenin  Equivalent 


(2-54a) 


1(0)  = -Yq  V(0) 

i(/)  = Yjf  y(r) 


Norton  Equivalent 


(2~54b) 


In  these  matrix  equations,  the  entries  in  the  1-th  rows  of  the  nyl  vectors 
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Clearly,  the  entries  in  and  Y f or  this  case  become 


!o- 


»«oi>  


0^ 

I 

I 

I 

0 


_0 


(2-57a) 


(1/Z^,)  0- 0 


0 (1/Zyo)  I 

1 ^ ^ "N) 


(2-57b) 


Note  that  Y^  = ^ and  Yj^=  case,  determining  the  entries 

in  the  terminal  impedance  (admittance)  matrices  is  a trivial  matter  and  the 
terminal  impedance  (admittance)  matrices  are  diagonal. 

The  more  difficult  case  occurs  when  each  wire  is  not  connected  directly 
to  the  reference  conductor  by  a single  impedance.  Two  examples  which 
illustrate  this  situation  are  shown  in  Figure  2-9.  First  consider  the 
situation  in  Figure  2-9 (a).  Here  it  is  obviously  not  possible  to  obtain 
terminal  Impedance  (admittance)  matrices  which  are  diagonal.  The  termination 
impedance  matrices  can,  however,  be  obtained  by  defining  loop  currents  In 
which  two  of  the  loop  currents  so  defined  are  the  terminal  currents 
and  12(0).  Writing  the  required  three  loop  equations  we  obtain 
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A,.  -g. 


V^iO)  = -Z^  - Z^Clj^  - I^)  - Z^(I^  + I^  - I^) 

V^(0)  - -Z2  - I^) 

0 = Z,  I + Z,(I  - L)  + Z,  (I  - I - I.) 

5c  3c  b 4c  3 b 

The  objective  is  to  eliminate  the  current  from  these  equations  leaving 

I = 1,(0)  and  I.  = I„(0)  as  a function  of  V (0)  and  V-(0).  The  third 
a 1 b z 1 / 

equation  yields 

(Z„  + Z, ) I.  + Z.  I 
i 4 b U a 

c " (Z^  + Z4  + Z5) 

Substituting  this  result  for  into  the  first  two  equations  eliminates 
the  current  from  these  equations  and  leaves 

Vi(0)  = Z^  1^(0)  + Zj^  I^CO) 

V^CO)  = Z^  1^(0)  + Z^  I^CO) 

where  Z^,  Z^,  Z^,  Z^  are  the  resulting  combinations  of  Z^^,  Z2,  Z^,  Z^,  Z^ 

and  we  have  substituted  1,(0)  = I , I„(0)  = I,. 

1 a Z b 

This  technique  can  obviously  be  generalized  for  any  number  of  wires  and 
additional  extraneous  loops  in  the  termination  networks.  Treating  the  n 
line  voltages  as  independent  sources  and  writing  the  required  number  of 
loop  equations  for  the  terminal  network,  we  may  obtain 
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-fr,^r>  ■ 


Vj(0) 


V„(0) 


nxn  j nxm 


nixn  I mxm 


I^CO) 


I (0) 
n 


(2-58) 


7'' 

te 


I 


In  (2-58)  we  may  eliminate  the  extraneous  loop  currents  I I bv 

1 m 

solving  the  second  set  of  equations  to  yield 


■ . “1 

- 

^1 

I;^(0) 

• 

, -1 

= - D C 

4 

• 

« 

f 

m 

• 

I„(0> 

. 

(2-59) 


Substituting  this  result  into  the  first  set  of  equations  we  obtain 


v^(0) 

• 

= (A  - B D~^  C) 

“ - 

1^(0) 

• 

• 

v^(0) 

« 

• 

I„(0> 

u «i 
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Clearly,  then  we  may  identify 


?0 


= - (A  - B d"^  C) 


(2-61) 


The  extension  of  this  technique  to  obtain  the  Norton  Equivalent  char- 
acterization employs  a dual  technique.  Here  we  define  node  voltages  (with 
respect  to  the  reference  conductor  at  either  x=0  or  x=I)  of  all  nodes  of 
the  termination  network  (including  the  n nodes  connected  to  the  line)  and 
write  the  node  voltage  equations  of  the  network  treating  the  line  currents 
as  independent  sources.  Therefore  we  write  (for  x=X) 


“ - 

1 

- 

* 

I^(^) 

{ 

1 

V^(rf) 

• 

• 

1 

A 

A 

• 

• 

• 

A 1 

B 

• 

I (t) 
n 

1 

nxn  . 

nxm 

JuflL 

8 

— ~ ~ ~ 

A 

0 

1 

V, 

A 1 

1 

• 

c 

1 

A 

D 

• 

• 

1 

“ 

• 

A 

0 

mxn  1 

mxm 

V 

in 

1. 

— 

A A 

Eliminating  the  extraneous  node  voltages  V^, we  obtain 


A 

"l 

• 

A 1 A 

= - d"  c 

• 

• 

A 

- 

• 

V 

vjt) 

mj 

n 

» ^ 

(2-62) 


(2-63) 


Substituting  we  obtain 


(2-64) 


A A A ^ A 

(A  - B d"  C) 


and  the  terminal  admittance  matrix  is  identified  as 


A A A *1  A 

= (A  - B D • C) 


(2-65) 


As  an  example  of  a Norton  Equivalent  formulation,  consider  the  termina- 
tion network  in  Figure  2-9  (b).  Here  we  may  write 


13(0  = (1/Z^)  VjU) 


= (I/Z2)  [V2(;()  - V^«)] 


I^(3J)  - (1/22)  [V^«)  - V2«)] 


Thus  we  may  identify  by  writing 


l2i£) 

I3(X) 


1/22 

-1/Z„  1/2^ 

^ 40 


V^(!f) 

V2«) 

V3(t) 


Note  for  this  example,  it  is  not  possible  to  obtain  the  Thevenln  Equivalent 
characterization,  Zj^,  since  Y^is  an  obviously  singular  matrix. 
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II T . DERIVATION  OF  THE  EXCITATION  SOURCES  FOR  UNIFORM  PLANE 


WAVE  AND  NONUNIFOKM  FIELD  EXCITATIONS 

In  the  previous  Chapter,  equations  for  the  terminal  currents  of  the 
line  were  derived  for  general  forms  of  the  excitation  field.  In  this  Chapter, 
we  will  derive  explicit  formulas  for  the  equivalent  induced  source  vectors 
for  uniform  plane  wave  excitation  of  TYPE  1 and  TYPE  2 structures.  The  co- 
ordinate system  and  reference  directions  for  the  incident  field  which  are 
assumed  by  the  program  will  be  indicated.  The  formulas  for  nonuniform  field 
excitation  which  assume  a spatial  piecewise  linear  characterization  of  the 
incident  field  will  also  be  derived. 

In  the  following,  some  confusion  may  arise  concerning  the  use  of  the 
word  "vector".  A spatial  or  physical  vector  will  be  denoted  as  E.  A 
matrix  or  column  array  vector  is  denoted  by  E.  These  two  "vectors"  are 
obviously  quite  different  quantities  however  the  word  "vector"  will  be  used 
for  both  with  the  distinction  between  the  two,  although  generally  obvious  from 
the  context,  being  denoted  by  an  arrow,  •>,  over  the  symbol  or  a bar,  -,  under 
the  symbol. 

The  equations  for  the  terminal  currents  of  the  line  for  all  structure 
types  are  repeated  here  for  convenient  reference.  If  the  Thevenin  equivalent 
characterization  of  the  terminal  networks  is  chosen: 

V(0)  = -Zq  1(0)  (3-la) 

va)  = ZjfI«)  (3-lb) 

then  the  equations  for  the  terminal  currents  are 
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-1 


[cos(ki)  {Z„  + Zi  + j sin(k«)  {Z^  + Z/ Z.)  ] 1(0) 
f (inc)  .C  4:  -0  - 

= M f j Z^Z-"  N - E^(gC)  + [cos(kr)  + j sin(k;C)  Z^  ] E^(0) 


I (a:)  = -j  Z^“^[n  + sln(k/r)  E^(0)] 

+ [cos(kX)  1 + j sin(kX)  Z~'^  Z„]  1(0) 


(inc) 


(3-2b) 


If  the  Norton  equivalent  characterization  of  the  terminal  networks  is  chosen: 


1(0)  = -Yp  V(0)  (3-3a) 

I(JC)  = Y^V(;f)  (3-3b) 


then  the  equations  for  the  terminal  currents  are 


(cos(k2)  {Yq  + U + j sin(kr)  (Y^Z^  Y + Z^  ^}  ] [-V(O)] 

^ ~ “vine)  “ 

« Y^M  + j Z_‘^  N - Y .E^Cf)  (3-4a) 

- .I-t 

+ [cos(k;<)  Y£+ j sin(k/)  z ' ] E (0) 

•0  ^ ^ ^ 


. (inc) 

!(>:)  » -j  Z ^ (N  + sin(kir)  E (0)  } 

• <*L*  *“v 

+ [cos  (ki)  Yq  + j sin(k<)  Z^’^]  I-V(O)] 


(3-4b) 


(inc)  (inc) 

The  nxl  induced  source  vectors  M,  N,  E^(0),  E^(;C),  in  these  equations 
are  defined  in  the  previous  Chapter  for  the  various  structure  types  and  the 
entries  in  these  vectors  are  due  to  the  incident  field.  It  is  the  purp^.s  "'f 
this  Chapter  to  derive  the  entries  in  these  vectors  for  uniform  plane  wave 
illumination  of  TYPE  1 and  TYPE  2 structures  and  nonuniform  field  illumination 
of  ail  structure  types. 


3.1 


Basic  Integrals 


(inc)  (inc) 

The  entries  in  the  Induced  source  vectors,  M,  N,  E^(0),  E^(/),in  (3-2)  and 


(3-4)  all  involve  Integrals  of  components  of  the  incident  electric  field 
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intensity  vector  along  certain  spatial  contours.  It  is,  of  course,  highly 
desirable  for  computer  implementation  to  obtain  closed  form  solutions  for 
these  integrals.  Throughout  the  following  derivations,  we  will  encounter  two 
fundamental  integrals  which  must  be  evaluated.  These  are  designated  as 
El(a,b,k)  and  E2(a,b,k)  and  are  given  by 


I 

»,b,k)  = J 


jkx 

X e dx 


(3-5a) 


E2(a,b,k)  =J  e^*^^  d:: 


(3-5b) 


The  straightforwcird  solutions  of  these  integrals  are 


El(a,b,k)  —as ) + (S ^ ) 


(3-6a) 


jkb  jka 

E2(a,b,k)  = (- ) 


(3-6b) 


Note  til, It  .v'lien  k~0,  evaluation  of  the  solutions  in  (3-6)  will  result  in 
obvious  prolilems.  Of  course,  the  Integrals  in  (3-5)  have  well  defined 


solutions  lor  k=0  and  these  are  quite  obviously 


El(a,b,0) 


k2  2 

b - a 


(3-7a) 


E2(a,b,0)  = b-a 


(3-7b) 


For  values  of  the  argument  k equal  to  zero,  the  program  evaluates  (3-7). 

The  following  solutions  for  the  entries  in  the  induced  source  vectors 
(inc)  (inc) 

M,  N,  E^(0),  E|.(jJ)  in  (3-2)  and  (3-4)  will  be  written  in  terms  of  these  integrals 
and  the  fundamental  Integrals  are  stored  in  the  program  as  function 
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subprograms.  (See  Section  4.2.) 


3.2  Derivation  of  the  Source  Vectors  for  Uniform  Plane  Wave  Illumination  and 

/ 

TYPE  1 Structures 

(inc)  (inc) 

The  basic  source  vector  quantities,  M,  N,  E^(0),  E^(X)>  involved  in 
the  equations  for  the  terminal  currents  in  (3-2)  and  (3-4)  for  TYPE  1 


structures  are  given  in  (2-36)  and  (2-37)  which  are 


(Inc) 

(inc)  (inc) 

lEj,(x)] 

1 ""  ~ hi 

(3-8a) 

(inc) 

f ^10  (inc) 

[Ej.(^)] 

i=J 

(3-8b) 

0 

(inc) 

/iO  (inc) 

[Ej.  (0)] 

i = J E^^(p^,0)dp^ 

(3-8c) 

0 

f* 

f (inc) 

M.| 

cos(k(jf-  x))  Ej^(x)  dx 

(3-8d) 

0 

f 

^ (inc) 

sln(k(;i^-  x))  Ejj^(x)  dx 

(3-8e) 

0 

(inc) 

(inc) 

where 

Ej^i(dio.x)  and 

E|^^(0,x)  are  the  components  of 

the  incident  electric 

field 

in  the 

x direction  (along  the  line  axis)  along 

the  1-th  wire  and  along 

(inc)  (inc) 

the  reference  wire,  respectively.  The  quantities  E^^(o^,3!)  and  E^^(p^,0) 
are  the  components  of  the  incident  electric  field  along  a straight-line 
contour  joining  the  i-th  wire  and  the  reference  wire  in  planes  (y,z) 
transverse  or  perpendicular  to  the  line  axis  at  x=Jf’and  x=0,  respectively. 
This  contour  is  denoted  by  p^. 

The  coordinate  system  used  to  define  the  wire  positions  and  shown  in 
Figure  2-4  is  used  to  define  the  angle  of  arrival  of  the  uniform  plane  wave 
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and  polarization  of  the  Incident  electric  field  intensity  vector.  In  de- 
fining the  wire  positions  for  ITPE  1 structures,  an  arbitrary  rectangular 
coordinate  syt  tern  is  established  with  the  reference  wire  at  the  center 
(y=0,  z-0)  of  this  coordinate  system  as  sho^m  in  Figure  3-1.  The  i-th  wire 
has  coordinates  y=y^,  z=z^,  relative  to  this  coordinate  system.  The 
direction  of  propagation  of  the  incident  wave  is  defined  in  Figure  3-2  by 
the  angles  0^  and  (J>^.  The  angle  0^  is  the  angular  orientation  of  the  electric 
field  intensity  vector,  E,  in  the  plane  containing  E (which  is  perpendicular 
to  the  direction  of  propagation)  and  measured  from  the  projection  of  the  y 
axis  onto  this  plane.  The  zero  phase  reference  is  taken  at  the  origin  of 
the  coordinate  system,  i.e.,  x=0,  y=0,  z=0. 

The  electric  field  intensity  vector  can  be  written  in  terms  of 
components  as  [12| 


J+E  y 
xm  ym 


, „ , -j  (k  X + k,,  y + k z) 

+ E z]e'’  x y-'  z 


The  Items  E , E , E are  the  magnitudes  of  the  projections  of  |(inc) 

in 

xm  ym  zm 

the  x,y  and  z directions,  respectively  and  x,y,z  are  unit  vectors  in  the 

x,y,  and  z directions,  respectively.  The  quantities  k , k , and  k are  the 

X y z 

components  of  the  propagation  constant,  k,  in  the  x,y  and  z directions, 
respectively.  To  determine  these  quantities,  note  that  the  electric  field 
intensity  vector  can  be  most  directly  related  to  a spherical  coordinate 
system  in  terms  of  the  unit  vectors  r.O,(j)  as  shown  in  Figure  3-2.  In  this 
spherical  coordinate  system,  we  may  write  [12] 


E = (E  r + E^  0 + E,  I) 

rm  0m  (t>m 


-3^ 


■> 

R 


(3-10) 
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TYPE 


I 


1 Yj  ^T>  H-r 

I 

I 

I 

I 

I 


1 (0,0) 


Figure  3-1.  The  TYPE  1 structure. 
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where  from  Figure  3-2 


E =0  (3-lla) 

rm 

= -E  cos  0_  (3-llb) 

0ra  m E 

'’e 

l!  - k r 

R=rr-xx+yy+zz  (3-lld) 

R is  a vector  from  the  origin  to  a point  P and  E is  the  magnitude  of  the 

m 

electric  field  intensity.  To  detenuine  the  components  E , E , E , k , 

xm  ym  zm  x 

k and  k in  (3-9)  we  simply  need  the  transformation  from  a spherical  co- 

y 2 

ordinate  system  to  a rectangular  coordinate  system  (see  reference  [12],  p.9). 
Employing  this  conversion  of  coordinate  systems,  we  find 


E = -E  cos  0„  cos  0 sin  di  - E sin  0^  cos 
xm  m E ppm  E 


E = E cos  0„  sin  0 
ym  m E p 


E = -E  cos  0_  cos  0 cos  <p  + E sin  0„  sin 
zram  E p pm  E 


k = k sin  0 sin  A 
x p "^p 


k = k cos  0 

y p 


k = k sin  0 cos  A 
z p ^p 


Calculation  of  the  quantities  in  (3-8)  proceeds  as  follows. 


entry  in  the  nyl  vector  M is 


r""  (Inc)  (inc) 

[M]^  = J cos(k(/-  x)){Ej^(d^Q,x)  - Ej^^(0,x)}  dx 


(inc) 


(3-12a) 

P 

(3-12b) 

(3-12c) 

P 

(3-12d) 

(3-12e) 

(3-12f) 

The  i-th 

(3-13) 
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where 


(inc) 


(inc) 

E^^(O.x)  = E 


- E I r 


^.y=yi 

z=z. 


(3-14) 


E e 
xm 


■jk^x  (^-3(ky  "i)  - 1 ) 


end  the  i-th  wire  has  y and  z coordinates  of  and  z^,  respectively.  Sub- 
stituting (3-14)  into  (3-13)  one  can  obtain 


{e  ^ ^^y  ^i  ^ ^z  ^i^  ~ J x))  e dx 

0 


f 


-jk  X 


-j(k.  y,  + k_  z,)  + k,)x 


= E {e  y ^i  ■ z "i^  -1}{„ 
xm  2 


(3-15) 


x"  dx 


-jki 


e"  f j(k-k)x,  } 
I e X dx  ■’ 


+ 2 J e 

1) 


This  result  can  be  written  in  terms  of  the  basic  intergral  E2  in  Section 
3 J as 


[M]^  = ^ le-^S  ^i  ‘^z  "i^  E2(0./.-(k  + k^)) 

+ e'^'^  E2(0,t»(k  - k ))  } 


(3-16) 


Similarly  the  entries  in  the  nxl  vector  N become 


[N]^  = 


I sin(k(3( 

J n 


(inc) 


(inc) 


- x))  {E„^(d^„,x)  - E^^(0,x)}  dx 


U'  iO 


= -J 


xm 


{e"^^Vi  -1}  E2(0,^-(k  + k^) 


(3-17) 


- E2(0,ar,  (k-k^))  } 
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The  calculation  of  the  entries 
proceeds  as  follows.  The  i-th  entry 


(inc) 

in  the  vectors  E (iL) 
(inc)  ^ 

in  E^(jt.)  is  given  by 


(inc) 
and  E^(0) 


(inc)  , ‘^iO  (inc) 
0 


where  is  a straight-line  contour  in  the  y,z  plane  (at  x=/)  joining  the 

reference  wire  and  the  i-th  wire.  The  i-th  wire  is  located  at  y=y,,  2=z. 

i 1 

and  the  reference  wire  is  located  at  y=0,  z=0.  The  center-to-center 

r~i  2 

separation  between  the  reference  wire  and  the  i-th  wire  is  d.»  = /y,  +z.  . 

lO  -^i  1 

Consider  Figure  3-3  which  shows  this  contour  and  the  appropriate  components 
of  the  electric  field  along  this  contour.  For  this  situation,  (3-18) 
becomes 


(inc)  f “^iO  ^ 

[E  (X)].=  / (E  cos  0 + E sin  0)  e ^ x*^ 
-t  i / ym  zm 

0 

r -j  (k  cos  0 + k sin  0)  p.  -i  , 
le-^'y  z "^  ijdp 


(3-19) 


where 


cos  0=3— 
'‘iO 


sin  0 


(3-20a) 


(3-20b) 


Therefore,  we  obtain 
(inc) 


-j(k,  y.  + k,  z.) 


^ ’ -ik  re  J ^ V ■’’i  7 I''  - U 

[E  tf)].=  (E  y , + E zj  e ^ x — 77^ — ll"  \ ^ (3-21) 

^-t'*"'i  'ym^i  zm  i"  ^’^y  >^1 


This  result  may  be  written  equivalently  in  terms  of  fundamental  Integral  E2 
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as 


(inc) 

[E  (<) ] .=  [E  y . + E 
-t  1 ym  1 zm 


I “jk  ■ 
z^J  e-  X 


E2(0,ly(k  y.  + k z.)) 
’ ’ y 1 z 1 ' 


(3-22) 


(inc) 

The  1-th  entry  in  the  vector  E^(0)  is  given  by  (3-22)  with  ^ = 0. 


r 

I 

I 


I 


I 


I 


i- 


3. 3 Derivation  of  the  Source  Vectors  for  Uniform  Plane  Wave  Illumination 
and  TYPE  2 Structures 

(inc)  (inc) 

Derivation  of  the  source  vectors,  M,  N,  E^(0),  E^(;C),  in  (3-2)  and 

(3-4)  for  uniform  plane  wave  illumination  of  n wires  above  a ground  plane 

(TYPE  2 structures)  proceeds  similarly.  Here,  we  will  cetermine  as 

the  net  electric  field  which  is  the  vector  sum  of  the  incident  wave  and 

the  wave  reflected  by  the  perfectly  conducting  ground  plane.  In  this  case, 

the  net  electric  field  tangent  to  the  ground  plane  will  be  zero.  Therefore 
(inc) 

Ej^^(0,x)  in  (2-49a)  will  be  zero.  Again,  an  arbitrary  rectangular  coordinate 
system  is  used  to  define  the  cross-sectional  positions  of  the  wires.  The 
ground  plane  forms  the  x,z  plane  (y=0)  as  shown  In  Figure  3-4.  The  zero 
phase  reference  for  the  incident  field  will  be  taken  to  be  at  the  origin  of 
this  coordinate  system,  i.e.,  x=0,  y=0,  z=0.  The  various  angles  defining  the 
direction  of  propagation  of  the  incident  wave  and  polarization  of  the  electric 
field  intensity  vector  are  the  same  as  for  TYPE  1 structures  and  are  shown 
in  Figure  3-2. 

The  primary  problem  here  is  to  determine  the  net  electric  field  parallel 
to  the  wire  axes  and  between  the  i-th  wire  and  the  ground  plane  along  a 
contour  perpendicular  to  the  ground  plane.  This  net  electric  field  is  the 
vector  sum  of  the  incident  field  (in  the  absence  of  the  ground  plane)  and 
the  portion  reflected  by  the  ground  plane. 


I 

I 

( 

I 

st 

I 

I 

f 
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We  may  write  the  incident  electric  field 


E ^ = (E^ 
xm 


-»•  1 ->•  1 
X + E y + E 
ym  zm 


-i  (k  X + k y + k z) 
z)  e X y ^ z 


(3-23) 


The  angle  of  reflection  between 
equal  to  the  ang3e  of  incidence 
mediately  write  the  form  of  the 


the  reflected  wave  and  the  ground  plane  is 
by  Snell's  Law  [12],  Therefore  we  may  im- 
reflected  wave  as 


t r 


. r ->■  r r 

(E  X + E y + E z) 
xm  ym  zm 


-j  (k  X 
e X 


k y + k z) 

y z 


(3-24) 


At  y=0,  continuity  of  the  tangential  components  of  the  electric  field  re- 
quire that 


r i A 

E = -E  = -E 
xm  xm  xm 


(3-25a) 


e’^  = -E*  = -E 
zm  zm  zm 


(3-25b) 


Consequently,  the  net  x component  ot  the  electric  field  is  given  by 

I-  = -j  (k  X + k y + k z)  , „r  -j  (k  x - k y + k z) 

=E  e-'  x y ^ z +E  e-'  x y z 

Xxotai  xm  ^ xm  ^ 


^ -i  (k  X + k z)(  -jk  y jk  y ? 
t e ' X z <e  - e-^  y-^  r 

xm  ( ' J 


(3-26) 


-2j  E sin(k  y)  e ^ ^^x  ^ ^z 
xm  y 


where  E^^  is  the  magnitude  cf  the  x comoonent  of  the  incident  electric 
field,  i.e.,  E^^  = E^.  Similarly,  one  may  show  that  the  net  y component 
of  the  eJectric  field  is  given  by  [12] 
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E „ „ /,  V -1  (k  X + k z) 

y„  ^ = 2 E cos(k  y)  e x z 

Total  ym  y 


(3-27) 


The  components  of  the  nyl  vector  M become 
.t  (inc) 


[M]^  = J"  cos(k(^'-  x))  [E^(x)]^  dx 


cos(k(/  - x))  E 


Total 


(3-28) 


= -2jE  e sin(k  y.)  f cos(k(Z'-  x))e  ^^x^  d> 

xm  y i I 

^ 0 

= -JE  e-iVl  sl..(k  y.)  f 

xm  y ^ ^0 

. -jkZ*  jkx  -jk  X,  , 

+ e-'  e"*  e'’x}dx 


-jE  e sin(k  y 

xm  y 


-jE^^  e'-'^'z^i  sin(ky  y^)  E2(0,4  -(k  + k^)) 

+ e"^^^E2(0,2’,  (k  - k^))} 


Similarly  calculation  of  the  entries  in  N yields 


• (inc) 

sin(k(;?’-  x))  [Ej^(x)]^  dx 


sin(k(/  - x))  E 


(3-29) 


Total 


= -E  e ^^z^i  sin(k  y.)  {e^^^E2(0,>f,-(k+k  )) 
xm  y i X 
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(3-29) 


-e’^'^  E2(0,t,  (k  - k^))} 


(inc) 

The  entries  in  E^(T<i)  are  given  by 


(inc)  /i 

IitW]i=  E 

^ i ^Total  ' 


E dy 

I ’z=z. 

x=X 


(3-30) 


or.  /I  \ ~j  (k  Jf  + k z . ) . 

2E  cos(k  y)  e X z i dy 

ym  ' y 


2E  e'^^x 


< -jk  z.  f 

e z i J 


1 jk  y , -jk  y 
e"^  V + e V 


e y + e 
2 


_ -jk 
• E e x 


y 

^ -jk  z.  r r ^i  jk  y f ^-jk  y , •, 

e z i { I e y dy  + I e dy) 


E e'’^*^x^e"’^*^z*l  { | e^^y^  dy  + f e^*^y^  dy} 
ym  J J 


E e ^^x^  e ^^z*i  {E2(-y. , y . , k )} 
■ym  11' 


(inc) 

The  entries  in  E^(0)  are  those  of  (3-30)  withX®  0, 

It  should  be  noted  that  the  above  quantities  can  be  determined  in  an 
alternate  fashion.  Rather  than  determining  the  net  electric  field  as  the 
sum  of  an  incident  and  a reflected  wave,  simply  replace  the  ground  plane 
with  image  wires  as  shown  in  Figure  3--5.  The  entries  in  the  source  vectors 
can  then  be  obtained  by  using  only  the  incident  field  and  treating  the  image 
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of  the  i-th  wire  as  the  "reference"  for  the  i-th  wire  as 


(Ml.  . I 


r (Inc) 

(inc) 

J cos(k(  x))  {E^ 

1 - E 

' \7=\r 

} dx 


0 


" ' 1 
z=z , 
1 


y=-y^ 

2=Z  . 


= cos  (k(^^  - 


wir-  “j(k  x + k y.+k  z.) 
x))  tE  e-^'x  y-'i  zi 

xm 


-j  (k  x - k y.+k  z ) , 

-E  e-^  x y i zi  ldx 

xm 


= E e ^'^z^i  (e  '^'^y^i  -e'^'^y^i)  f cos(k(/-  x))  e ^"^x^  dx 
xm  J 

0 


«-2J  E 

xm 


i sin(k^  y^)  [ cos(k(2’ - x))  e ■^^x^  dx 


H E_  e'^'^z^i  sin(k  y ) {e^‘'^E2(0,Jf,-(k  + k )) 

y 1 X 


xm 


+ e"^'''^E2(0,Jf,(k  - k^))} 


[N]  - f sln(k(;^  - x))  fE 

^ A 


(Inc) 

X 


(Inc) 


- E 


y=yj 

z=z . 


) dx 


y=-yj 
z=z . 


.2 

: I sln(k<2-  X))  (E  \ >i  + ''x  'P 

J xm 


E -j(k  x k y.+k  z.),, 
xmex  yi  zi)dx 


-E  e ^^z^i  sin(k  y.)  E2(0,Z,-(k  + k )) 

xm  y 1 ' > > 


- e"^‘'^E2(0,2:,(k  - k^))} 


(3-31a) 


(3-31b) 


-bO- 


(inc) 


(3-31c) 


z=z. 

1 


=/ 


„ ~j  (k  + k y + k z ) 
E e*^  x Y"^  zi 


ym 


dy 


„ -jk  -jk 
E Q X e z 
‘‘ym 


jk  y . 

y dy 


-y. 


= V ( f dy  + j dy} 


0 


/u 

jk  y , 
e y dy } 


-y^ 


" V e*j''z''l  {K2(-y^,y^,ky)} 


which  are  precisely  the  results  obtained  previously. 

3,4  Calculation  of  the  Source  Vectors  for  Nonunlform  Fields 

Nonuniform  field  excitation  can  be  specified  for  all  structure  types. 
The  problem  here,  again,  is  to  evaluate  equations  of  the  form  in  (3-8). 
This  requires  tliat  we  specify  values  (magnitude  and  phase)  of  the  electric 
field  intensity  vector  along  the  wires  and  reference  conductor  and  between 
each  wire  and  the  reference  conductor  at  the  endpoints  of  the  line.  To 
accomplish  this,  we  will  specify  values  at  a finite  number  of  points  along 
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I 


the  appropriate  contours  and  assume  piecewise-linear  variation  of  the  electric 
field  (magnitude  and  phase)  between  the  specified  points.  This  is  illustrated 
in  Figure  3-6. 

^ (m) 

For  TYPE  1 strucutres,  the  values  of  E,  E^^,  along  the  reference  wires 

(in  the  +x  direction)at  N. „ + 1 points  will  specified  as  shown  in  Figure 

. (">) 

3-6 (a).  The  values  of  E,  E . , along  the  i-th  wire  at  N , + 1 points  will 

(m)  (m) 

be  specified.  The  values  of  E,  and  E^,  at  x=0  and  x=.z'along  a straight- 

line  contour  in  the  y,z  plane  joining  the  reference  wire  and  the  i-th  wire  at 
+ 1 and  1 points,  respectively,  will  be  specified.  Similar  quan- 

tities will  be  specified  for  TYPE  2 and  TYPE  3 structures  as  shown  in 
Figure  3-6 (b)  and  Figure  3-6 (c),  respectively,  with  the  exception  that  fi 
is  taken  to  be  zero  along  the  reference  conductor  for  these  two  cases. 

Piecewise-linear  variation  of  the  electric  field  (magnitude  and  phase) 
is  assumed  between  these  specification  points  as  shown  in  Figure  3-7  where 
the  magnitude  of  the  appropriate  component  of  the  electric  field  is  denoted 
by  I* I and  the  angle  is  denoted  by • Thus  the  problem  is  the  determination 
of  quantities  of  the  form  in  (3-8)  for  this  piecewise-linear  variation  of  the 
field.  The  technique  is  to  write  linear  equations  representing  the  piece- 
wise-linear variation  of  the  magnitude  and  phase  of  the  field  between  suc- 
cessive specification  points  and  add  the  appropriate  integrals  over  the 
adjacent  regions. 

The  first  problem  then  becomes  to  characterize  the  linear  magnitude 

and  phase  variation  between  two  successive  data  points.  Consider  two 

successive  data  points,  x^  and  which  specify  the  magnitude  of  the 

electric  field,  E and  E , , and  phase,  0 and  0 respectively.  Knowing 
m m+l  m m+1  > ^ 

Che  end  points,  one  can  write  linear  equations  characterizing  the  linear 
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Figure  3-6.  Nonuniform  field  specification  for  TYPE  1 
Structures. 
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onunlform  field  specif icacior  for  TYP 
tructures • 


Ground  Plane 


Dehavior  between  these  successive  points  as 


|E„(x) 

m 


0^(x) 

m 


= a X + b 
m ra 


c X + d 
m m 


E - E 
f m+1  m 1 

a ={  } 

m X , , - X 

m+1  m 


E X - E X 
, m m+1  m+1  m •, 

b =1 i 

m X , , ~ X 

m+1  m 


0,-0 

c -f  ^ ] 

m X . 1 - X 
m+1  m 


0 X . - 0 , , X 

, r m m+1  m+1  m, 

d =l } 

m X , , - X 

m+1  m 


(3-32a) 


(3-32b) 


(3-33a) 


(3-33b) 


(3-33c) 


(3-33d) 


The  electric  field  is  then  characterized  by 


E (x)  = |e  (x)|  e^®m^^^ 


(3-34) 


The  quantities  of  the  form  in  (3-8)  which  must  be  evaluated  involve 
certain  integrals  involving  the  form  (3-32),  The  first  type  is  of  the 
form  in  (3-8d)  which  can  be  evaluated  as 


(3-35) 


m+1 


m 


cos(k(;^-  x))  {a  X + b } e'^  ^ ^ 


m m 


dx 


X 


m+1 


Jd 


e“  m {a  / cos(k(;^-  x))  x dx 

m J 

X 

m 


/ 


m+l 


+ I cos(k(2:'-  x))  dx) 


e'^  m ; jk/  , , , 

= -T  {a  El(x  ,x  ,c  -k) 

2 m m m+1  m 


+ a e El  (x  ,x  , , ,c  +k) 

m m m+1  m 


+ b E2(x  ,x  ,c  -k) 

m m m+i  m 


+ b e^^^E2(x,x,,,c  +k)} 
m m m+1  m 


The  second  form  is  similar  to  (3-8e)  which  can  be  evaluated  as 

f sin(k(>f-  x))  {a  x + b } e^  ^*^m  ^ ^ '^m^  c 

J mm 


m 


e ^'^m  {a  El(x  ,x  ,c  -k) 

m m m+1  m 


m m m+I  m 


(3-36) 
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+ b E2(x  ,x  ,i,c  -k) 

m m m+1  m 


- b e ^^^E2(x  ,x  ,,  ,c  +k)} 
m m m+1’  m 


The  third  forms  are  similar  to  (3-8b)  and  (3-8c)  which  can  be  evaluated  as 


/ ..^3(cx  + d), 

(a  X + b ) e m m dx 

m m 


(3-37) 


= e"^  m {a  El (x  , x , , , c ) 
m m m+1  m 


+ b E2 (x  ,x  ^ ,c  )} 
m m m+1  m 


The  program  compC-es  the  items  in  (2-36)  and  (2-37)  for  TYPE  1 
structures;  (2-49)  and  (2-37)  for  TYPE  2 structures;  and  (2-52)  and  (2-37) 
for  TYPE  3 structures  by  breaking  up  the  appropriate  integrals  and  adding 
the  integrals  between  each  pair  of  successive  data  points.  The  data  points 
need  not  be  equally  spaced  along  the  appropriate  contours  so  that  one  can 
model  localized,  extreme  variations  in  the  fields  without  using  an  in- 
ordinately large  number  of  data  points.  In  specifying  the  sequence  of 
electric  field  components  along  each  contour,  one  must  insure  that  the 
first  and  last  specification  points  are  at  the  two  ends  of  the  contour  and 

m m+1 


IV.  COMPUTER  PROGRAM  DESCRIPTION 


The  contents  and  operation  of  the  code  will  be  described  in  this 
Chapter.  The  cards  in  the  program  deck  are  sequentially  numbered  in 
columns  73-80  with  the  word  WIRE  in  73-76  and  the  card  number  in  77-80. 

The  program  is  written  in  Fortran  IV  language  and  is  double  precision.  A 
listing  of  the  program  is  contained  in  Appendix  B and  a general  flow  chart 
of  the  program  is  given  in  Appendix  3.  In  this  flow  chart,  the  numbers  on 
the  left  and  right  of  the  individual  boxes  denote  the  beginning  and  ending 
card  numbers  of  the  corresponding  portion  of  the  code  listing,  respective >> . 
Changes  in  the  program  to  convert  to  single  precision  arithmetic  are  in- 
dicated in  Appendix  B.  The  program  has  been  implemented  on  an  IBM  370/165 
digital  computer  at  the  University  of  Kentucky  using  the  Fortran  IV  G-Ievel 
compiler. 

The  program  requires  two  function  subprograms,  El  and  E2,  and  one  sub- 
routine, LEQTIC,  which  must  follow  the  main  program  and  precede  the  data 
cards.  Subroutine  LEQTIC  is  a general  purpose  subroutine  to  solve  a set 
of  n,  complex,  linear  simultaneous  equations  and  is  a part  of  the  IMSL 
(International  Mathematical  and  Statistical  Library)  package  [13],  If  the 
IMSL  package  is  not  available  on  the  i .ex's  system,  other  appropriate 
general  purpose  subroutines  may  be  used,  (See  Section  4.2  for  a discussion 
of  LEQTIC  and  its  argument  list.)  Listings  of  function  subprograms  El 
and  E2  are  contained  in  Appendix  C. 

4.1  Main  Program  Description 

A listing  of  the  WIRE  program  is  contained  in  Appendix  B.  Cards  0001 
through  0055  contain  general  comments  concerning  the  applicability  of  the 
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program.  Cards  0057  through  0060  contain  the  array  dimension  Information. 
All  vectors  and  matrices  should  be  dimensioned  to  be  of  size  N where  N 
is  the  number  of  wires  exclusive  of  the  reference  conductor.  These  matrices 
and  vectors  must  be  dimensioned  appropriately  for  each  problem.  Cards  0062 
through  0069  declare  double  precision  real  and  complex  variables  and 
dimension  the  vector  and  matrix  arrays. 

Cards  0071  through  0080  define  certain  constants, 


CMTM  (conversion  from  mils  to  meters)  - 2.54  x 10 


MU02PI 


ONE80 


ZEROC 


V(velocity  of  light  in  ';ee  space) 


1.  4 jO. 

0.  + jO. 

0.  + jl. 

2.997925  x 10®  m/sec 


RADEG (conversion  from  radians  to 
degrees) 


180. /ir 


Note  that  n is  computed  to  the  user’s  machine  precision  by  using  the 


relationship 


tan  (ir/4)  = 1 


Cards  0086  through  0145  read  and  print  portions  of  the  input  data  and 
perform  certain  primitive  error  checks.  These  cards  read  the  structure 
type  (TYPE“l,2,or3),  the  load  specification  option  (LS0=ll,12,21,or22) , 
the  field  specification  option  (FS0=l,or2),  the  number  of  wires,  N,  the 
relative  permittivity  of  the  medium,  ER,  the  relative  permeability  of  the 
medium,  MUR,  and  the  line  length,  L.  In  addition,  for  TYPE  1 structures, 
the  radius  of  the  reference  wire,  RWO,  and  for  TYPE  3 structures,  the 
interior  radius  of  the  overall,  cylindrical  shield,  RS,  are  read. 

Cards  0150  through  0197  read  the  radii,  r^^,  and  the  and  y^  coordinates 
(r^  and  0^  for  TYPE  3)  for  the  N wires  and  compute  the  entries  in  the 
characteristic  Impedance  matrix.  The  z and  y coordinates  are  stored  in  the 
real,  nx".  vectors  V3  and  V4,  respectively; 

'z  for  TYPE  1,21 
, f V3(I) 

r^  for  TYPE  3 ) 

’ y for  TYPE  1,2  7 

. > -*■  V4(I) 

0^  for  n'PE  3 I 

In  addition,  the  entries  in  the  nyn,  real  characteristic  impedance  matrix, 

Zj,,  are  temporarily  stored  in  the  real  parts  of  the  nxn  complex  matrix 
Ml,  This  is  done  to  minimize  the  required  array  storage  in  the  program 
since  the  matrix  Ml  will  be  needed  later  as  a complex  matrix. 

Matrix  Array 

Zj,  Ml  (real  part) 
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Cards  0202  through  0211  compute  the  inverse  of  Z^,  Z„  which  Is 
stored  in  the  real  parts  of  the  nxn,  complex  arrays  M2  and  M3.  Subroutine 
LEQTIC  computes  this  inverse  by  solving  the  system  of  equations  AX  ■ 1 
where  1 is  the  nyn  Identity  matrix.  The  solution  X is  therefore  A 

«w  m0 

(See  Section  4.2  for  a more  complete  discussion  of  subroutine  LEQTIC.) 

Since  the  real  part  of  Ml  contains  the  characteristic  impedance  matrix,  the 
real  part  of  the  nxn,  complex  solution  matrix,  M2,  will  contain 
Therefore 


Matrix 


Cards  0222  through  0254  read  and 
Impedance  (admittance)  matrices  at  x 
These  are  stored  in  the  nxn,  complex 

Thevenin  Equivalent 

!o 

U 


Array 


M2  (real  part) 

M3  (real  part) 

print  the  entries  in  the 

- 0,  Zq  (Yq),  and  at  x • 

arrays  as 

terminal 

Norton  Equivalent 

Array 

!o 

YO 

It 

YL 

Cards  0259  through  0267  Interchange  the  entries  in  arrays  Ml  and  M2 
if  the  Thevenin  equivalent  characterization  is  chosen.  Thus 


Thevenin  Equivalent 


Norton  Equivalent  Array 


!c 

Ml 

M2 

-1 

M3 

«C 
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Cards  0273  through  0292  compute  the  quantities: 


Thevenin  Equivalent 

5c  * 5f  5c’'  5o 


Norton  Equivalent  Array 

Yfl  + M2 


The  contents  of  this  array  M2  will  be  retained  throughout  any  frequency 
iteration  so  that  these  matrix  products  need  be  computed  only  once. 

Cards  0293  and  0294  compute  the  phase  constant  for  a frequency  of  one 
hertz  and  its  product  with  the  line  length: 


Quantity 


Variable 


2ii 


'l  Hertz 


(vq/zht;) 


BB 


k<| 


Hertz 


BBL 


To  obtain  the  propagation  constant  at  each  frequency,  BB  is  multiplied  by 
the  appropriate  frequency. 

Cards  0300  through  0323  read  the  input  data  describing  the  uniform 
plane  wave  if  FSO  *•  1,  i,e,,  E^,  0^,,  0^,  and  compute  the  x,  y,  z 
components  of  the  electric  field  and  propagation  constant  (for  one  Hertz) 
as  shown  in  (3-12), 

Cards  0327  through  0333  read  the  frequency  on  the  fltst  frequency  card 
and  compute  the  propagation  constant,  k,  k^,  sln(kX),  cos(k/): 
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Quantity 

Variable 

k 

BETA 

k t 

BETAL 

sin (k<) 

DS 

cos(k/) 

DC 

If  uniform  plane  wave  field  specification  Is  selected  (FSO  1),  cards 

0342  through  0351  compute  certain  preliminary  quantities  to  be  used  in 

(inc)  (inc) 

computing  the  induced  source  vectors  M,  N,  E.  (0),  ^(jC).  If  TYPE  1 
structures  are  selected,  cards  0357  through  0370  compute  the  Items  In 
these  induced  source  vectors  as: 

Vector  Array 


M (3-16) 

VI 

N (3-17) 

V2 

(inc) 

E (0)  (3-22,i  = 0) 

ETO 

(inc) 

E^{£)  (3-22) 

ETL 

If  TYPE  2 structures  are  selected,  the  Items  In  the  Induced  source  vectors 
are  computed  in  cards  0375  through  03R6  as; 


Vector 

Array 

M (3-28) 

VI 

N (3-29) 

V2 

(Inc) 

E (0)  (3-30, 0) 

ETO 

Tine) 

E^(£)  (3-30) 

ETL 
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If  nonuniform  field  excitation  (FSO=2)  is  selected,  cards  0399  through 

(Inc)  (inc) 

05A2  compute  the  entries  in  the  Induced  source  vectors  M,  N,  E (0),  E.  (rf). 

Cards  0399  through  0439  read  the  magnitude  and  phase  of  the  incident  electric 

field  at  specification  points  along  the  reference  wire  and  compute  the 

portions  of  M and  N along  the  reference  wire  if  TYPE  1 structures  are 

selected.  If  TYPE  2 or  TYPE  3 structures  are  selected,  this  computation  is 

bypassed  since  for  these  types  of  structures  it  is  assumed  that  the  total 

electric  fields  tangent  to  the  ground  plane  and  the  interior  wall  of  the 

cylindrical  shield  are  zero.  Cards  0444  through  0479  read  the  magnitude 

and  phase  of  the  incident  electric  field  at  specification  points  along  the 

wires  and  compute  the  entries  in  M and  N for  each  wire  which  are  stored  in 

arrays  VI  and  V2,  respectively.  (The  l~th  entries  contain  the  results  for 

the  i-th  wire.)  Cards  0484  through  0510  read  the  magnitude  and  phase  of 

the  incident  electric  field  at  specification  points  along  contours  in  the 

y,  z plane  between  the  reference  conductor  and  the  1-th  wire  at  x*0  and 

(inc) 

compute  the  entries  in  £,.(0)  which  are  stored  in  array  ETO.  (The  electric 

field  is  tangent  to  these  contours.)  Cards  0515  through  0542  repeat  this 

(inc) 

calculation  for  x^Z  and  compute  the  entries  in  E Of)  which  are  stored  in 

— t 

array  ETL. 

Cards  0548  through  0592  form  the  equations 

(cos(k<)  { Z-  + Z.]  + j sin(k/)  { Z_  + Zy  Z Z_}  ] 1(0)  - 

(inc)  -C  .l.C  ,0  - 

- M + j Zjj  N - Ej.(Z)  + [cos(k<)  + j sin(k;C)  Z^^Z^.  £^.(0) 

for  the  Thevenin  Equivalent  specification  (LS0=11,12)  or 
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[cos(kr)  { Yq  + Y^}  + j sln(k/)  { Yq  + Z~^  }]  [-V(O)]  = 

1 (Inc)  ^ (inc) 

= Y^M  + j ^ N - YjjE^C^)  + [cos(kX)  Y^+  j sin(k<)  £^.(0) 

for  the  Norton  Equivalent  specification  (LSO=21,22),  The  coefficient  matrix 
is  stored  in  array  A and  the  right-hand  side  vector  of  the  equations  is 
stored  in  array  B.  The  arrays  VI,  V2,  ETO,  ETL  contain 


Vector 

Array 

(Inc) 

M - Ej.(i) 

, (Inc) 

?c’  StW' 

VI 

ETL 

(Inc) 

Et(0) 

ETO 

N 

V2 

The  main  diagonal  entries  in  the  array  Ml  contain  Z-  ^ N 

••V 

Vector  Array 

Z^  ^ N Ml  (cn  main  diagonal) 

Subroutine  LEQTIC  is  called  for  the  solution  of  these  equations  in  card 
0596  and  the  solution  vector  (1(0)  for  LSO^ll.lZ  or  -V(0)  for  LSO«21,22) 
is  returned  in  array  B. 

The  terminal  currents  are  computed  in  cards  0610  through  0641.  Cards 
0610  through  0618  compute  the  quantities: 


Array 

WA 

Cards  0619  through  0631  compute  the  terminal  currents  at  x=J^  for  the 
Thevenin  Equivalent  from 

, (inc)  ^ 

Ui)  = -j  f N + sln(k<)  E^(0)}  [cos(kt)  1^  + j sin(k<)  Zq]I(0) 

and  for  the  Norton  Equivalent  from 

1 (inc)  ^ 

I(Tf)  • -j  ■"{  N + sin(k^)  £^.(0)  }+  [cos(k^)  + j sin(k^)  Z^'^]  r-V(O)] 

A. 2 Subroutine  LEQTIC 

Subroutine  LEQTIC  is  a general  subroutine  for  solving  a system  of  n 
simultaneous  complex  equations.  The  program  is  a part  of  the  IMSL 
(International  Mathematical  and  Statistical  Library)  package  113). 

The  subroutine  solves  the  system  of  equations 

A X = B (4-1) 

where  A is  an  nxn  complex  matrix,  B is  an  nxm  complex  matrix  and  X is  an 
nxm  complex  matrix  whose  columns,  X^,  are  solutions  to 

A 5l  = 2l 

where  B^  is  the  1-th  column  of  B. 

The  calliny  statement  is 

CALL  LEQT1C(A,N,N,B,N,M,WA,IER) 
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where 


A -»■  A 
B B 
N n 
M m 

and  WA  is  a complex  working  vector  of  length  n.  lER  is  an  error  parameter 
which  is  returned  as^ 

lER  = 128  ->•  no  solution  error 

lER  =■  129  A is  algorithmically  singular  [13]. 

The  solution  X is  returned  in  array  B and  the  contents  of  array  A are 
destroyed. 

Subroutine  LEQTIC  can  be  used  to  find  the  inverse  of  an  nxn  matrix  by 
computing 

A X = 1 (4-3) 

~ - ~n 

where  1^  is  the  nxn  identity  matrix.  Thus  the  solution  is  X “ A LEQTIC 
is  used  in  numerous  places  to  invert  real  matrices  by  defining  the  real  part 
of  A to  be  the  matrix  and  the  imaginary  part  to  be  zero.  Upon  solution, 
the  real  part  of  X is  the  inverse  of  tho  real  matrix,  A. 

4.3  Function  Subprograms  El  and  E2 

Function  subprograms  El  and  E2  are  used  to  evaluate  (in  closed  form)  the 

^ The  solution  error  parameter  is  printed  out  whenever  A is  singular.  The 
error  is  IER-128  so  that  the  solution  error  will  be  l^when  A is  singular. 
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commonly  occurring  integrals; 


u 

El  (a,b,k)  = / X 
a 

E2  (a,b,k)  d: 


jkx 

e dx 


(4-4) 


(4-5) 


Function  E2  can  be  evaluated  as 


E2  (a,b,k)  = 


j kb  j ka 
e - e 


(4-6) 


which  can  be  written  in  an  alternate  form  as 


jk(b+a) 

E2  (a,b,k)  - e ^ (— 


.ik(b+a) 

2 /e 

a e ( — 


ik  (b-a ) -.ik(b-a) 


(b-a) 


(^7) 


■)  (b-a) 

2 


sin  { j 

(b-a)  k(b-a)  e 
2 


This  form  of  E2  is  more  attractive  from  a computational  standpoint  since  the 

sln(X)/X  expression  in  the  final  result  can  be  computed  quite  accurately 

for  small  values  of  the  argument  whereas  the  form  in  (4-6)  may  suffer  from 

roundoff  errors  when  k is  small.  In  fact, a test  was  conducted  on  the  IBM 

370/165  in  double  precision  by  computing  the  function  sin(X)/X  for  values 

—78 

of  X=l,  .1,  .01,  .001,  , 10  until  exponential  underflow  occurred. 

Tlie  results  converged  to  the  expected  value  of  1.  In  fact, for  values  of 
X from  10  ^ to  10  the  result  was  ip—  ^ . 


15  digits 

Note  that  these  Integrals  can  be  analogously  viewed  as  Fourier  Transforms. 
Although  this  concept  is  interesting,  it  provides  no  significant  help  since 
the  evaluation  of  these  integrals  can  be  easily  obtained  in  a straightforward 
manner  without  resorting  to  a table  of  transforms. 
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The  function  E2  is  computed  in  the  function  subprogrsim  E2  with 


argument  list 


E2  (A,  B,  X) 


X 

The  quantity  DIF=B-A  is  computed  as  well  as  the  quantities  FA=j  DIF  and 
FB=j  (B+A).  If  FA=0,  the  program  evaluates  E2*DIF  since  sin(FA)/FA  “ 1. 

If  not,  the  program  evaluates  E2=DIF  {sin(FA) /FA}  e^^^. 

Finding  a more  suitable  computational  form  for  El  is  considerably  more 
complicated.  El  can  be  evaluated  as 

b 


El  (a,b,k)  = 


J dx 


(4.8) 


, j kb  j ka 
be-*  - ae'' 


gjkb  _ ^jka 


This  result  can  be  separated  into  ..  real  and  imaginary  part,  i.e., 

El(a,  b,  k)  =»  RE  + j IM  (4-9) 


where 


RE  a cos (kb)  + kb  sin (kb)  - co3(ka)  - ka  sin(ka)  (4-lOa) 

k^ 


jU  - sln(kb)  - kb  cos(kb)  - 8in(ka)  4-  ka  cos(ka)  (4-lOb) 

k^ 


The  real  part  can  be  written  as 
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RE  = 


[2  sin  f ^ ] 


(4-11) 


b sin (kb)  - a sln(ka) 


(b  - a ) 


sin  ( '‘-S|=al) 


k(b4-a) 

2 


k(b-a) 

2 


2 sin( 


2 sin(ka) 


Note  that  this  form  of  the  real  part  of  El  contains  only  sin(X)/X  expressions 
and  can  be  computed  very  accurately  for  small  X,  Notice  that  as  k-^O,  the  real 
part  becomes 


RE  » 
k-K) 


k2  2 

b - a 


(4-12) 


which  is  precisely  the  value  of  El (a,  b,  k)  when  k=0.  Therefore  the 
imaginary  part  of  El  must  go  to  zero  as  k goes  to  zero. 

The  Imaginary  part  of  El  can  be  written  as 


sin (kb)  - kb  cos (kb)  - sln(ka)  + ka  cos(ka) 


f „ X sin(ka)-v  , r sin(kb)  , 

at  cos(ka) — ) - b { cos  (kb) } 


(4-13) 


Note  that  as  k->0  th»»re  is  a distinct  possibility  of  roundoff  error  in 


computing  the  function 


cos  (t)) 


sin(0) 


(4-14) 
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1 


RE  = 4 f2  f J 

k 


(4-11) 


b sin (kb)  - a sln(ka) 
k 


/i,2  2, 

- (h  - a ) 
2 


k(b4-a) 

2 


k(b-a) 

2 


¥ 


, , 2 sin(kb)  2 ain(ka) 
kb  ‘ ^ ka 


Note  that  this  form  of  the  real  part  of  El  contains  only  sin(X)/X  expressions 
and  can  be  computed  very  accurately  for  small  X.  Notice  that  as  k->-0,  the  real 
part  becomes 


RE 


k2  2 

b - a 


(4-12) 


which  is  precisely  the  value  of  El (a,  b,  k)  when  k=0.  Therefore  the 
imaginary  part  of  K1  must  go  to  zero  as  k goes  to  zero. 

The  imaginary  part  of  El  can  be  written  as 


IM 


sin (kb)  - kb  cos (kb)  - sin(ka)  * ka  cos(ka) 


(4-13) 


af  cos(ka)  - — ^ - b { cos  (kb)  - 


ka 


sin (kb) 
kb 


Note  that  as  k->0  there  is  a distinct  possibility  of  roundoff  error  in 
computing  the  function 


cos(t))  - 


sin(0) 

0 


(4-14) 
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converges  to  1.  Clearly,  as  k-H),  the  imaginary  part  of  El,  IM,  conver^^es 
to  zero,  as  expected. 

We  will  select  a value  of  |0|  = .01  = .573®  as  the  point  below  which 
we  evaluate  a truncated  portion  of  (4-16).  The  tradeoff  here  is  to  select 
a value  of  l0|  small  enough  so  that  a truncated  portion  of  (4-16)  will  not 
require  many  terms  for  sufficient  accuracy  yet  )o|  is  not  too  small  to  re- 
sult in  round  off  error  when  evaluating  (4-13)  directly.  We  have  selected 
the  value  to  be  |0|  = .01  and  will  truncate  the  series  in  (4-16)  to 


(1  - — ) 
^ 10  ^ 


(4-17) 


For  |0|  “ .01,  the  cos  0and  terras  in  (4-14)  are  Identical  to  only 

4 digits.  This  should  provide  sufficient  accuracy  to  prevent  any  significant 
roundoff  error  in  (4-14).  In  evaluating  (4-16),  we  will  obtain  accuracy 
to  10  digits  by  using  the  truncation  in  (4-17).  (Note,  for  0 « .01,  0^/lO»lO”^, 
0^/280  *3.57  X 10  0^/15120  * 6.61  x 10  and  terras  with  higher 

powers  of  0 affect  only  those  digits  well  to  the  right  of  16  places.)  Thus 
this  criterion  seems  to  provide  sufficient  accuracy  while  limiting  the 
roundoff  error  in  evaluating  IM,  Therefore,  our  result  is 


IM  - IMA  - IMB 


(4-18) 


where 


IMA  = dos(ka)  - } |ka!  >.01 


(4-19a) 


{ 1 - 


I ka  < .01 


i! 

m; 


5.1  Transmission  Line  Structure  Characteristics  Cards.  Group  I 

WIRE  considers  (n+1)  condwctor  transmission  lines  consisting  of  n 
wires  in  a lossless,  homogeneous  surrounding  medium  and  a reference  con~ 
ductor  for  the  line  voltages.  The  n wires  and  the  reference  conductor 
are  considered  to  be  perfect  (lossless)  conductors.  There  are  three 
choices  for  the  reference  conductor  type; 

TYPE  ■ 1;  The  reference  conductor  is  a wire. 

TYPE  ■ 2:  The  reference  conductor  is  an  Infinite 

ground  plane. 

TYPE  * 3;  The  reference  conductor  Is  an  overall, 
cylindrical  shield. 

Cross-sectional  views  of  each  of  these  three  structure  typas  are  shown  in 

Figure  5-1,  5-2  and  5-3,  respectively. 

For  the  TYPE  1 structure  shown  in  Figure  5-1,  an  arbitrary  rectangular 

coordinate  system  Is  established  with  the  center  of  the  coordinate  system 

at  the  center  of  the  reference  conductor.  The  radii  of  all  (n+1)  wires, 

r . , as  well  as  the  Z and  Y coordinates  of  each  of  the  n wires  serve  to 
wi 

completely  describe  the  structure.  Negative  coordinate  values  must  be 
Input  as  negative  data  items,  for  example,  and  Y^  In  Figure  5-1  would 
be  negative  numbers. 

For  the  TYPE  2 struc  re  shown  in  Figure  5-2,  an  arbitrary  coordinate 
system  Is  established  with  the  ground  plane  as  the  Z axis.  The  coordinates 
Y^  and  Y^  ( positive  quantities)  define  the  heights  of  the  i-th  and  J-th 
wires,  respectively,  above  the  ground  plane.  The  necessary  data  are  the 
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TYPE  « I 


TYPE  • 3 


Figure  5-3.  The  TYPE  3 structure. 
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Z and  Y coordinates  and  the  radius,  r^^,  of  each  wire. 

For  the  TYPE  3 structure  shown  in  Figure  5-3,  an  arbitrary  cylindrical 
coordinate  system  is  established  with  the  center  of  the  coordinate  system 
at  the  center  of  the  shield.  The  necessary  parameters  are  the  radii  of  the 
wires,  the  angular  position,  0^,  the  radial  position,  r^,  of  each 

wire  and  the  interior  radius  of  the  shield,  r^. 

The  format  of  the  structural  characteristics  cards.  Group  X,  Is 

shown  in  TABLE  1.  The  first  card  contains  the  structure  TYPE  number 

(1,  2,  or  3),  the  load  stiucture  option  number,  LSO,  (11,  1’'.,  21,  or  22), 

the  field  specification  option  number,  FSO,  (1  or  2),  the  number  of  wires, 

n,  the  relative  dielectric  constant  of  the  surrounding  medium  (homogeneous), 

e^,  the  relative  permeability  of  the  surrounding  medium  (homogeneous),  p^, 

and  the  total  length  of  the  transmission  line,5f,  (meters).  If  TYPE  1 or  3 

is  selected,  a second  card  Is  required  which  contains  the  radius  of  the 

reference  wire,  r^Q,  (mils)  for  TYPE  1 structures  or  the  Interior  radius 

of  the  shield,  r , (meters)  for  TYPE  3 structures.  For  TYPE  2 structures, 
s 

this  card  is  absent.  These  cards  are  followed  by  n cards  each  of  which 
contain  the  radii  of  the  n wires,  r^^,  (mils)  and  the  and  Y^^  coordinates 

of  each  wire  (meters)  for  TYPE  1 and  2 structures  or  the  angular  coordinates 
r^  (meters)  and  0^^  (degrees)  of  the  i-th  wire  for  TYPE  3 structures.  These 
n cards  must  be  arranged  In  the  order  1 ■ 1,  i ■ 2,  — , 1 ■ n. 


5,2  The  Termination  Network  Characterlasation  Cards.  Group  II 

This  group  of  cards  conveys  the  terminal  characteristics  of  the  term- 
ination networks  at  the  ends  of  the  line,  x ■ 0 and  x ■^.  The  termlna  ‘on 
networks  are  characterized  by  either  the  Thevenln  Equivalent  or  the  Norton 
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TM>LK  1 


Format  of  the  Structure  ChnrnctorlHtlcH  CnrdK,  Croup  I 


Card  Group  //I  (total  ■ 1); 

card 

column 

format 

(a) 

TYPE  (1,2,3) 

10 

T 

(b) 

LOAD  STRUCTURE  OPTION, LSO, 

(11,12,21,  or  22) 

19 

- 20 

T 

(c) 

FIELD  SPECIFICATION  OPTION,  FSO, 

(1  or  2) 

30 

1 

(d) 

n (number  of  wires) 

39 

- 40 

I 

(e) 

. relative  dielectric  v 

^constant  of  the  surrounding  medium"' 

41 

- 50 

E 

(f) 

^relative  permeability  of  the 
^r  ^surrounding  medium) 

51 

- 60 

E 

(B) 

it  (line  length  in  meters) 

61 

- 70 

E 

^ ^ ^ //O  /total  - 1 if  TYPE  - 1 or  3 . 

Card  Group  itl  (^g^al  - 0 if  TYPE  - 2 ^ 

(«) 

T/yp  - 1.  /’Radius  of  reference. 

6 

- 15 

E 

(b) 

TYPE  - 2:  absent 

(c) 

- /interior  radius  of. 

TY  E • 3.  ^shield  in  meters  ^ 

6 

- 15 

E 

Card 

Group  113  (total  ■ n) 

(a) 

r .(wire  radius  in  mils) 
wi  

6 

- 15 

E 

(b) 

for  TYPE  1 or  2 in  meters 

r^  for  TYPE  3 in  meters 

21 

- 30 

E 

(c) 

Y,  for  TYPE  ,1  or  2 in  meters 

36 

- 45 

E 

0^  £or  TYPE  3 in  degrees 


Note;  Cards  In  Group  113  must  be  arranged  in  the  order; 
wire  1,  wire  2,  ,,,,  wire  n 
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Equivalent  characterization. 

V(0)  - - Zq  1(0) 
V(3t)  - 2^1(0 


These  characterizations  are  of  the  form 

")  Thevenin 

(5-la) 

Equivalent 


1(0)  =>  - Yq  V(0) 

Hi)  = 


Norton 


Equivalent 


(5-lb) 


(See  Volume  VII  or  Chapter  II  for  a more  complete  discussion  of  determining 
the  entries  In  Z-.,  Zy  , Y.  and  Yy  [2].)  The  transmission  line  consists  of 
n wires  which  are  numbered  from  1 to  n and  a reference  conductor  for  the 
line  voltages.  The  reference  conductor  is  numbered  as  the  zero  (O)  con- 
ductor. Thus  Zq,  Yq,  are  nxn  matrices  which  are  assumed  to  be 
symmetric.  The  n entries  in  each  of  the  nxn  vectors,  V(0)  and  V«),  are 
the  line  voltages  with  respect  to  the  reference  conductor  at  x»0  and  x» 
respectively.  The  n entries  in  each  of  the  nxn  vectors,  1(0)  and  I(jf), 
are  the  line  currents  at  x-0  and  x»3f,  respectively.  The  currents  at  x-0 
are  directed  out  of  the  termination  networks  whereas  the  currents  at  x*  ^ 
are  directed  into  the  termination  networks.  The  entries  in 
these  four  vectors  are  arranged  in  the  order  wire  1,  wire  2,  — , wire  n. 

The  impedance  or  admittance  matrices,  Z-  and  Z^or  Y~  and  Y^, 
respectively,  may  either  be  "full*'  in  which  all  entries  are  not  necessarily 
zero  or  may  be  diagonal  in  which  only  the  entries  on  the  main  diagonals  are 
not  necessarily  zero  and  the  off-diagonal  entries  are  zero.  The  user  may 
select  one  of  four  LOAD  STRUCTURE  OPTIONS  (LSO)  for  communicating  the 
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entries  in  the  vectors  and  matrices  In  (5-1).  These  are: 


LSO  * 11 


LSO  - 12 


LSO  » 21 


LSO  - 22 


f Thevenln  Equivalent  representation;  7 
( diagonal  impedance  matrices,  Zq  and  Z^.  ) 


( Thevenln  Equivalent  representation 
( full  impedance  matrices,  Z^  and  Z 

( Norton  Equivalent  representation;  -j 

c diagonal  admittance  matrices,  and  Y^.  j 

S Norton  Equivalent  representation 
C full  admittance  matrices,  Y^  and 


" I 


The  structure  and  ordering  of  the  data  In  Group  II  are  given  In  Table 
2 and  can  be  summarized  In  the  following  manner.  The  first  group  of  cards 
In  Group  II,  Group  11(a),  will  describe  the  entries  on  the  main  diagonal 
in  Yq(Zq),  Tf  j^(Zjj)»  These  cards  must  be  in 

the  order  from  1 <■  1 to  1 « n.  Each  of  these  entries  is  In  general,  com- 
plex. Therefore  two  card  blocks  are  assigned  for  each  entry;  one  for  the 
real  part  and  one  for  the  Imaginary  part.  For  example,  consider  a 4 con- 
ductor line  (3  wires  and  a reference  conductor).  Here  n would  be  3, 
Suppose  the  Thevenln  Equivalent  c'<'«racterlzation  Is  selected,  with  the 
following  entries  In  the  characterization  ro8ttic''>s: 


Z-  ■ 
-0 


7 + J8  0 

0 j9 

L 0 0 


0 

0 

10  + jll, 


h 


16  0 0 

0 17  + J18  0 

0 0 jl9 
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TABLE  2 


Group  11(b)  ( 


total  = n(n“l)/2  If  OPTION  = 12  or  22 
total  = 0 If  OPTION  = 11  or  21 


) 


card  column 


format 


^Oij  ^Slj  ^ - 

1 

" reax  part 

1-10 

E 

. imaginary  part 

11  - 20 

E 

[ real  part 

41  - 50 

E 

} 

, imaginary  part 

51  - 60 

E 

Note;  If  LSO  ■ 12  or  22,  a total  of  n(n-l)/2  cards  must  be  present  and 
must  follow  Group  II (a) « If  LSO  • 11  or  21,  this  card  group  Is  omitted. 
The  cards  must  be  arranged  so  as  to  describe  the  entries  in  the  upper 
i txlangle  portion  of  Yq(Zq)  and  Y by  rows.i.e.,  the  cards  must 

f contain  the  12  entries,  the  13  entries,  — the  In  entries,  the  23 

r 

i entries,  — , the  2n  entries,  - — etc.  The  ordering  of  the  cards  is 

i 

therefore : 

wires  1,2 
wires  1,3 

• 

wires  l,n 
wires  2,3 
wires  2,4 
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One  wou3.d  have  selected  LSO=ll,  The  n*3  cards  would  be  arranged  (in  this 
order) 


Group 

7.E0 

8.E0 

16.  EO 

O.EO 

11(a)  - 

O.EO 

9.E0 

17.  EO 

18.  EO 

10.  EO 

11.  EO 

O.EO 

19.  EO 

Card 

column 

!o 

^0 

+ 

50 

+ 

60 

If  the  terminal  Impedance  matrices  were  not  diagonal,  e.g*,  LS0*12  is 
selected,  then  n(n~l)/2  additional  cards,  Group  11(b),  would  follow  the 
above  n cards  comprising  Group  11(a).  These  cards  describe  the  entries  in 
the  upper  triangle  portion  of  the  termination  impedance  or  admittance 
matrices  by  rows.  Suppose  the  networks  are  characterized  by  the  Z^nd  Z^r 
matrices : 


The  following  n(n>l)/2  * 3 cards  must  follow  the  above  3 cards  in  the 
order  of  the  12  entries  first,  the  13  entries  next  and  then  the  23  entries s 


I 1 


t 

^ 7 + j8 

20  + j21 

22  + j23' 

i 5o  • 

20  + 131 

J9 

24  + j25 

' 

■ 

,22  + j23 

24  + j25 

10  + jll,. 

1 

mm 

16 

26  + j27 

28' 

1 

26  + J27 

17  + jl8 

J29 

1 

- 28 

j29 

jl9. 

1 
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20.  EO 

21.  EO 

26.  EO 

27. EO  * 

Group 

11(b)  S 

22.  EO 

23.  EO 

28.  EO 

O.EO 

24.E0 

25.  EO 

O.EO 

29.  EO 

Card 

t 

+ 

column 

10 

20 

50 

60 

5.3  The 

Field  Specification  Cards, 

Group  III 

There  are  two  Field  Specification  Options  (FSO)  for  specifying  the  form 
of  the  excitation  fl^ld: 


FSO  - 1 


Uniform  plane  wave  illumination 

of  the  llne(TyPEl  or  TYPE  2 structures 

only) 


1 


FSO  ■ 2 


Nonuniform  field  illumination  of 
the  line  (all  structure  types) 


] 


I 


5.3  .1  Uniform  Plane  Wave  Illumination.  FS0»1 


For  uniform  plane  wave  Illumination  of  the  line,  FS0»1,  the  format  of 

the  data  cards  is  shown  in  Table  3 and  consists  of  two  card  groups.  Card 

Group  #1  consists  of  one  card  containing  the  magnitude  of  the  electric 

field  intensity  vector,  E , the  angle  between  this  vector  and  the  projection 

m 

of  the  y axis  on  the  plane  containing  f (this  plane  is  perpendicular  to  the 
propagation  direction),  the  angle  between  the  y axis  and  the  direction  of 
propagation,  and  the  angle  between  the  z axis  and  the  projection  of  the 
propagation  vector  onto  the  x,z  plane.  (See  Figure  5-4.)  The  x coordinate  is 
parallel  to  the  n wires  and  reference  conductor,  and  the  y,z  plane  foms 
the  cross-section  of  the  line.  The  origin  of  th'_  coordinate  system, 
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TABLE  3 

Format  of  the  Field  Specification  Cards.  Group  III,  for  Uniform  Plane  Wave 


Illumination,  FS0=1.  (See  Figure  5-U) 


Card  Group  //I  (total  = 1); 

card  column 

format 

(a) 

E 

m 

(magnitude  of  the  electric  field 
intensity  vector  in  volts/meter) 

1-10 

E 

(b) 

0*. 

E 

(angle  of  electric  field  Intensity 
vector  in  desrees) 

16  - 25 

E 

(0 

0 

P 

(angle  of  propagation  direction 
from  y axis  in  degrees) 

31  - 40 

E 

(d) 

♦p 

(angle  of  projection  of  propagation 
direction  on  the  x,z  plane  from  z 
axis  in  degrees) 

46  - 55 

E 

Card  Group  #2  (total  ■ unlimited) 


(a)  Frequency  of  incident  wave  in  Hertz 


1-10 


E 


Group  I.  (See  Figure  5-1  and  5-2.)  The  zero  phase  of  the  incident 


! 

wave  is  taken  at  the  origin  of  this  coordinate  system. 


Card  Group  #2  for  FS0=>1  consists  of  ai»  unlimited  number  of  cards  with 


each  frequency  of  the  incident  wave  on  each  card.  More  than  one  frequency 
card  may  be  included  in  this  frequency  card  group.  The  program  will 
process  the  data  provided  by  Groups  I and  II  and  the  wave  orientation  data 
in  Group  #1  in  Table  3 and  compute  the  response  at  the  frequency  on  the 
first  frequency  card.  It  will  then  recompute  the  response  at  each 
frequency  on  the  remaining  frequency  cards.  The  program  assumes  that  the  data 
on  card  Groups  I and  IZ  and  the  wave  orientation  data  in  Group  in  Table 
3 are  to  be  used  for  all  the  remaining  frequencies.  If  this  is  not  intended 
by  the  user,  then  one  may  only  run  the  program  for  one  frequency  at  a time, 
X'hls  feature,  however,  can  be  quite  useful.  If  the  termination  networks 
are  purely  resistive,  l.e.,  frequency  independent,  then  one  may  use  as  many 
frequency  cards  as  desired  in  Group  #2  and  the  program  will  compute  the 
response  of  the  line  at  each  frequency  without  the  necessity  for  the  user 
to  input  the  data  in  Groups  I and  II  and  the  wave  orientation  data  for  each 
additional  frequency.  Many  of  the  time-consuming  calculations  which  are 
Independent  of  frequency  need  to  be  computed  only  once  so  that  this  mode 
of  useage  will  save  considerable  computation  time  when  the  response  at  many 
frequencies  is  desired.  If,  however,  the  termination  network  characteristics 


(in  Group  II)  are  complex-valued  (which  implies  frequency  dependent),  one 
must  run  the  program  for  only  one  frequency  at  a time. 

5.3  .2  Nonuniform  Field  Illumination,  FS0»2 

The  format  of  the  Field  Specification  Cards,  Group  III,  for  nonuniform 


I' 

\ 


I 
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field  illumination,  FSO=2,  is  shown  in  Table  4.  The  first  card  group. 

Group  #1,  cpnslsts  of  one  and  only  one  card  which  contains  the  frequency 
of  the  field. 

The  remaining  cards  contain  the  values  of  the  longitudinal  electric 
field  (magnitude  and  phase)  along  the  n wires  (and  reference  wire  for  TYPE 
1 structures)  which  are  directed  in  the  -he  direction,  and  the  transverse 
electric  field  along  straight  line  contours  joining  the  i-th  wire  and  the 
reference  conductor  at  x»0  and  x«^.  The  directions  of  the  transverse 
field  at  these  specification  points  are  tangent  to  the  contours  and  directed 
from  the  reference  conductor  to  the  l-th  wire.  For  TYPE  1 structures,  the 
precise  location  and  orientation  of  the  transverse  field  specification 
contours  should  be  clear.  For  TYPE  2 structure!,  the  transverse  field  spec- 
ification contours  should  comprise  the  shortest  path  in  the  y, z plane  between 
the  ground  plane  and  the  1-th  wire,  l.e.,  it  should  be  perpendicular  to 
the  ground  plane  or  directly  beneath  the  1-th  wire.  For  TYPE  3 structures, 
the  transverse  field  specification  contours  should  comprise  the  shortest 
path  in  the  y,z  plane  betv;een  the  interior  wall  of  the  cylindrical  shield 
and  the  i-th  wire, (See  Figure  5-5.) 

The  ordering  of  the  card  Groups  ^2-^9  is  quite  logical  but  somewhat 
Involved  to  describe.  The  philosophy  of  the  ordering  is  as  follows.  If 
TYPE  1 structures  are  selected,  we  first  describe  the  longitudinal  electric 
field  (magnitude  and  phase)  along  the  reference  wire  at  (Nj^^-fl)  specification 
points.  This  is  done  in  Groups  #2  and  #3.  (If  TYPE  ■ 2 or  3,  Groups  #2 
and  //3  are  omitted  and  it  is  assumed  that  the  net  Incident  electric  field 
is  obtained,  l.e.,  the  electric  field  tangent  to  the  ground  plane  (TYPE*2) 
and  the  interior  of  the  cylindrical  shield  (TYPE»3)  is  zero.)  In  Group  #2, 
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(N|j) 


we  communicate  the  number  and  the  magnitude  and  phase  of  the  field  at 
the  first  specification  point  at  x»0.  In  Group  #3,  we  compute  the  locations 
of  the  remaining  N^q  specification  points  and  the  magnitude  and  phase  of  the 
field  at  each  of  these  specification  points.  The  cards  in  this  group  must  be 
arranged  sequentially  so  that  each  specification  point  Ls  located  to  the  right 
of  the  previous  point.  In  addition,  the  last  (N^q+1)  specification  point 
must  be  equal  to  the  line  length,/,  i.e.,  located  at  • 

The  remaining  card  groups  (jl>4-//9)  use  the  same  philosophy  as  Group 
#2  and  //3  and  describe,  for  each  wire  from  1 to  n,  the  quantities  (in  this 
order) : 

(1)  longitudinal  field  on  i*-th  wire, 

(2)  transverse  field  at  x«0  between  the  reference  conductor  and  1-th 
wire,  and 

(3)  transverse  field  at  X"X  between  the  reference  conductor  and  1-th 
wire. 

For  example,  after  Group  #3  we  must  have  Groups  #4-#9  for  wire  1,  Groups 
#4"#9  for  wire  2,  — , Groups  #4-//9  for  wire  n.  This  is  Illustrated  in 
Figure  5-6, 

It  should  be  noted  that  the  incident  electric  field  which  one  specifies 
in  Card  Groups  #2-#9  la  the  incident  field  with  the  n wires  (and  the 
reference  wire  for  TYPE  1 structures)  removed.  This  is  Inherent  in  the 
derivations  of  Chapter  II.  Thus  cne  specifies  the  longitudinal  electric 
field  at  points  along  the  positions  of  each  wire. 
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TABLE  4 (continued) 


Format  of  the  Field  Specification  Cards.  Group  HI,  for  Nononlform  Fields, 

FSO  ■ 2 (see  Figure  5-5) 


Card  Group  //I  (total  » 1); 


card  column 


format 


(a)  Frequency  of  incident  field  in 
Hertz 


1-10 


n A n /total  «■  1 if  TYPE  - 1 ^ 

C«d  Group  »2  . 2 or  3> 


(a)  (number  of  field  specification 
points  along  reference  wire 


1-10 


(b) 


(0) 


"£0' 


(magnitude  of  electric  field 
along  reference  wire  in  4>x 
direction  at  x«0  in  volts/meter) 


21  - 30 


(c) 


■rO) 


(phase  of  electric  field  along 
reference  wire  at  x«0  in 


41  - 50 


degrees) 


0—4  r, . at  /total  - if  TYPE-1  x 
<.b.ont  lf%PE  . 2 or  3> 


(a)  X 


(»«) 


(b)  |E 


(ejec'.rlc  field  specification 
point:  along  reference  wire  in 
maters) 

(m) 

,qI  (mc'vnltude  of  electric  field 
at  J.n  •■‘•X  direction  in 
voltfe)!"<»ter) 

(m) 


1-10 


21  - 30 


(c)  /E 


(phase  of  electric  field  at 


’'lo  degrees) 


41  - 50 


Note!  ra  * 1,2,— snd  x^^q  (the  last  specification  point)  must  equal 


the  line  length,  / . The  c trds  in  Group  #3  must  be  arranged  such  that 

( !i)  ^(o+-l) 


lu 


:io 
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TABLE  4 (continued) 
Card  Group  #6  (total  » 1)  continued 


format 


(c)  /e^q^  (phase  of  electric  field  on 
* contour  at  x»0  in  degrees) 


Card  Group  #7  (total  » N^qj^) 


(a)  (electric  field  specification 

point  on  contour  between  re- 
ference conductor  and  i-th  wire 
at  x«0  in  meters) 

(b)  l^foll  (magnitude  of  electric  field  on 

contour  at  x"0  at  p ^5'  in 
volts/meter) 

(c)  (phase  of  electric  field  on 

fc  ■ " - contour  at  in  degrees) 


41  - 50 


1-10 


21  - 30 


41  - 50 


Note;  m ■ 1,2,— and  (the  last  specification  point)  must  be 

located  at  the  center  of  the  1-th  wire  at  x»0.  The  cards  in  Group  #7  must 


be  arranged  such  that 


(m)  (m+1) 

^tOi  ^tOi 


Card  Group  #8  (total  ••  1) 


(a)  N^j^  (number  of  field  specification 

points  on  straight  line  contour 
at  x"  i.  between  reference  conductor 
and  i-th  wire  • N + 1) 


Ie^iI  (magnitude  of  electric  field 

on  contour  at  x«  'i  in  volts/meter) 


1-10 


21  - 30 


(c)  /e^^^  (phase  of  electric  field  on  contour  41  - 50 
4 — - — at  x»i  in  degrees) 
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1 I 


TABLE  4 (continued) 

Card  Group  #9  (total  ■ card  column 


(a)  (electric  field  specification  point 

* on  contour  between  reference  con- 


ductor and  1-th  wire  at  vX  in  meters) 


(b) 


,(m) 


'y'l  (magnitude  of  electric  fie 
contour  at  x«)^at  P^)  in 
volts /meter) 


field  on 


1-10 


21  - 30 


41  - 50 


format 


E 


E 


(c)  /e^^]  (phase  of  electric  field  on  contour 
' ' ■ at  x*;iat  in  degrees) 

Note;  and  p^^^  (the  last  specification  point)  must  be 
located  at  the  center  of  the  i-th  wire  at  x*X • The  cards  in  Group  iP9 
must  be  arranged  such  that 


(m) 

'til 


(m+1) 

'tii 


Note;  Card  Groups  )!*4  - 9 must  be  repeated  for  wires  1 to  n and  arranged 
sequentially  for  wire  1,  wire  2,  - — » wire  n« 
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VI.  EXAMPLES  OF  PROGRAM  USAGE 


Several  examples  of  program  usage  will  be  described  in  this  Chapter. 
These  examples  will  serve  to  Illustrate  preparation  of  the  data  input  cards 
as  well  as  provide  partial  checks  on  the  proper  functioning  of  the  program. 
The  data  input  cards  as  well  as  the  computed  results  will  be  shown  for  each 
example. 

6.1  Example  I 

In  this  Chapter  we  will  show  an  example  of  a two  wire  line  above  a 
ground  plane  (TYPE>2)  illuminated  by  a uniform  plane  wave  (FSO^l).  The 
solution  for  the  terminal  currents  for  a 1 volt/meter  field  with  several 
angles  of  incidence  will  be  shown.  The  image  problem  will  also  be  con- 
sidered by  replacing  the  ground  plane  with  the  Images  of  the  wires  resulting 
in  a four  wire  line  (N*3,TYPE«1).  The  corresponding  currents  in  the  wires 
for  the  problem  of  two  wires  above  a ground  plane  should  be  twice  those  for 
the  image  problem. 

Two  Wires  Above  a Ground  Plane 

The  problem  considered  here  is  shown  in  Figure  6-1.  Wire  ill  has  a 
radius  of  30  mils  and  is  5 cm  above  the  ground  plane.  Wire  H has  a 
radius  of  10  mils  and  is  2 cm  above  the  ground  plane.  The  two  wires  are 
separated  horizontally  by  4 cm.  The  cross-section  of  wire  #1  is  located 
at  y»5  cm,  z»0.  The  cross-section  of  wire  #2  is  located  at  y«2  cm,  z»4  cm. 
The  line  length  is  5m  and  e^»l  (a  logical  choice  although  any 

and  may  be  used  in  the  program).  Each  wire  is  terminated  with  a single 
impedance  (in  this  case  purely  resistive)  between  the  wire  and  the  ground 
plane.  Clearly  one  may  chose  the  load  structure  option  of  LS0«11  with  the 
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terminal  impedance  matrices 


‘lOO 

■■ 

0 

7 BS 

"500 

0 ' 

_ 0 

500. 

_ 0 

1000. 

Three  orientations  of  the  incident  field  will  be  considered: 

(a)  - 1 V/m,  0_  » 30%  - 150%  4 - 40“ 

m . Cl  p p 

(b)  E^  - 1 V/m,  0g  - 0%  0p  » 90%  - 90“ 

(c)  E - 1 V/m,  0„  « 0“,  0 • 180“,  ^ « 90“ 

m £>  p p 

Notice  that  case  (b)  has  the  wave  propagating  in  the  +x  direction  along  the 
line  axis  with  ll  in  the  +y  direction,  i.e., 


E ■ e 


"Jkx 


k 

‘I 

a- 

. 


Case  (c)  has  the  wave  propagating  broadside  to  the  line  (in  the  -y  direction) 
with  E in  the  +x  direction,  i.e., 

f ■ t 

Four  frequencies  of  excitation  will  be  investigated: 

1 MHz,  10  MHz,  100  MHz,  iGHz 

(1E6),  (1E7),  (1E8),  (1E9) 

and  since  the  loads  are  resistive,  the  frequency  iteration  feature  of  the 
program  can  be  used  by  simply  placing  all  four  frequency  cards  as  a group 
at  the  end  of  the  program. 

The  reason  for  using  these  frequencies  is  that  for  1 MHz,  the  cross- 
sectional  dimensions  of  the  line  are  electrically  small.  For  1 GHz,  they 
are  not.  This  will  serve  to  further  illustrate  why  we  require  that  the  cross- 
sectional  dimensions  of  the  line  be  electrically  small.  To  illustrate  this 
let  us  arbitrarily  select  the  distance  between  wire  #1  and  the  image  of  wire 
#1  to  be  the  ’’largest"  cross-sectional  dimension  of  the  line.  This  distance 
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is  given  cy 


The  quantities  kd 
^ max 

d ® 10  cm 

max 

(in  degrees)  at  the  above  four  frequencies 

frequency 

kd  (degrees) 

max  ** 

d /X 
max 

1 E6 

.12 

.000334 

1 E7 

1*2 

.003336 

1 E8 

12.0 

.033356 

1 E9 

120.0 

.333564 

Notice  that  for  the  frequency  of  1 E9,  the  cross-sectional  dimensions  of  the 
line  are  certainly  not  electrically  small.  For  the  other  frequencies,  they 
probably  are. 

The  input  data  cards  for  the  angles  of  Incidence  in  (a),  (b)  and  (c) 
are  shown  in  Figure  6-2 (a),  (b)  and  (c),  respectively.  The  results  are 
shown  in  Figure  6-3. 

6.1.2  Two  Wires  Above  a Ground  Plane  by  the  Method  of  Images 

Here  we  solve  the  problem  considered  in  the  previous  section  by  the 
method  of  Images.  The  image  problem  becomes  a four  wire  problem  (N«3, 
TYPExl}  as  shown  in  Figure  6-4.  Here  we  choose  (arbitrarily)  th?  image  wire 
of  wire  ifl  in  the  previous  problem  as  the  reference  wire.  The  various 
wire  radii,  and  coordinates  are: 


...o.Jig, 

ttaii«iin««iiiia»itiiiifi)«ffii»nBMiiiip»iii*«i«««i«i«fl«iiMiiitijti»iiuiiii«iiir«iM»iiii««fiii«iiMiiiiniiiiii 

I I 1 « t i|l  j’i*<>H«tti}aii»nttnNa»PM»Mitn»M»«HM]»4»fln4iMA«io««i»Mu»t«>BV)i»itRnaHNfiflNNiiiia))XffaiiNaH 
t I I 4 « *1}  «n 

iiiii.hWt. 

mU 


I 1 1 4 1 ti 


I » •««»»»  ••Him  aw  niiinimii,it)i»im»i,,*,oo„„„,,^^ 

soo  eo  J e- Wif 

^^ifl^pS*”^*"******^  *************** **************  «8||8 
a»»»axtfM»4t  hU 


HBBH1B»n»«BBMBNnB«B»»NIIB«BBa« 
ll»BHB»"6l«l’«l«MII«l'ltN»»«IIBIIIinBB« 
««««»IIBBHB»»BBMI<BB«aB>iaBB»*aNI|K»BIIN 


I I I < I I I • ••1*131»0»»»*»»BI(»Hb*«B8II»X»»BB'  _ 

I > I • I • I I ••  ***??**"’'"''"*"»»»'>»'»»»»»»9^^«8«»«»»»»»MB»(IB»»»»BaBBMBB»»»IlBII»)l»Blt 

I I I < I • I ( •l"»»»«»»«4l»n»».»»IIBBl»BaB»B»BBl««B««B»Mf»8BB»Blifc58 


Mt 


""""""'""JIM  t Mil  I III  Miin  . Ill  I iiiii:iiiiiiiiii 

nnnunnnnn-unn^urnnilt;?^^ 

nnjinjnninnnnnnijnjftfJ^M-j^li'^nnjjn.nnjjininjnninn 

4 4 4 4 4 4 44  4 4 mi  4 4 4 44  4 4 4 M 4 44  M«l^l4^H,ra0M44MII4Mm«44l44|«||4,««,4,,«44 

MMMMMIIIIIMIIIIlMMIMMIlfjMjliilllViil^yillllllllllllllllllllllllllMIII 

innnii!nir>nnnnnnnniii^iUM,j{;n>mnnjimnMnnnmnnnn 


• " " " " " " y w ij  m n j 7 j T n n n I n n n n I J n n I n 

••miM.iiiiiiiiiiii 

?ii!l!!l!!******************************’’***”****^***>*>niniiiiiiiiiiiiiiiiii 


■Utn.TIWOW 


Figure  6-2 (a).  Data  cards  for  the  problem  in  Figure  6-1 


with  E - IV/m,  0 •30®,  0 -ISO®,  4 -40®. 
“•  £>  P P 
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} I * » tl 

I 1 « t 

I 1 4 1 • ^ « n 


iiaiiHtti«i)«tinfiitiiHnMffii»M»ttii)fifi»unii4«««**N««N««i«itifuwifMiJ»MHiitfu««Mii«H»iiiinMiiiini«i«ii 

• amiantitaiianHnitiiiiaxiiiiaKiiiiwiiitwuiiaMiiviiaHNiiuuHiiHVHiiMiiaiiKaHiiHWitnnnNiiiiiiiiiiN 

• ••«aannNaatia»a«aaMaK»a»<««au(i««««<iN»HMM*xiiMa«ii«aMaHiiMM*iiaaN>i«»«nH 


iii»*Maani4n»naaaaaaM»»»aN««naH««««<«M«HUHKi>uHNN«««HaHva«»ii  niMaar  . 

1 1 < « •^^^)"*ax*»»«<aa(4Ctaaaaa«aaaaKlk$tfaaaaaaaaatlbCtaaaaaaaeaaaaaaaaana  . t • • "v 

1 I •^^^'^^"•''"''»^4^8<'»»»naBaBa»anaa»»a»a«Bo«5lls?JiiaaMaal!4iM»anaaaBaaaiia»ai,«^»aa 
11.1.1  "•^'?*M«imaaaiii.»i»l^,«aaaa»aaal»Waaaa»»8UBia.iaaa«aaaaaBaa»i.»aa»»n»aa 
■•^^»^^*»*'*»>'»»»»:«J»?t««»B»«»»»«^l5:4a»aaaB-»»BaB« 


!!!*!!!*********************************HIMIIII'||IM||||  liiiiitiflittttxi 

1 1 m 1 1 1 1 1 n n M n‘<!' 

innjnj'nnnnjnnjnMMnuXc^rj^miffMjiiMMnMMMjnninnimi 

»niiiJnmnjnnjn>njnn>6*r>ffivV|»ii.ijVfvnijjjj'j}nnn)’mjijn>n 

H44<4l44444M444n  4 4 4 4 M4«MuhiM«4U?i4444lf*»ll4«4IMm4ll44n4l44IHIUI4 

iiMMiHMi4in4ii4umiiMiu4nMiUunM'ffiiiinMiinin,,,,,,,,,,,,,,,,| 

'''.nMMnMnnnnninnnnV.v'yinnrTnjTnnnMMnnnMinn.nni 


Figure  6-2 (b). 


Data  cards  for  the  problem  in  Figure  6-1 

with  E » 1 V/m,  G -0  , 0 -90  , (p  -90  . 
m E p p 
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m 


I 1 1 a I tS\^^v\m»a»ii**tntnp»ntinntni*nitita»tiumimitam»ini9m»taaunua»»m»timmmm»nmnnnnnn»n‘»nm 
I I I a I >\\Sl^vaaH»ajin*n*n«nHa«mnnpaaaaaaaaaaaaaamaaanaaaaaapaaaaaaaamaaannanaa»aann 
I 1 1 a I )^^(nta««<ii(«ax9tCllaxaHaxa»xUfk{Sa)a«ai«aiaiaa«t{lk((a»«aa«aaMaa««»»»«»N«»»»*a 

o.a^9M5„.,«,.t»W  D It » i«  n If  ti  N M N II  It » It » II  ir » Ji*i « « « iltfltl tlfl....aaaa. 

,,,a^?9^^S....,.nCta..a..a.a....a....aa....9lf|{l 

«ii4«Nii*ati 

vnaaMtiaaMna  m,  a aaaaaaaa^&Cataaaaaaaaaaaaaaaaaaaaaaaaa  a V ’***''*<"' 

...  a 1 1 . I ..•aa«»).«tliaaaaaaaaa(«aaaaaaaa(l«aaaaalkfilaaaHaattCiaaaaaa(i{t>aaaaaaaaa 


t Mill IIMII III MIIIIMI IIMMII M tlMMIMMIII  'MIMIMriHIMIII  'MIIMIMI 
1 1 a o . I ••n«aaa«aa«aaaaaaaaaaa«aaaaa«hati*|rtaaaattiiaaa8aMBaaaaaaaaaaaataaRaaaaaaaaa 

I III  I III  I III  1 1 111  “ n 1 1 1 1 1 1 1 ‘ III  iH(^iV  iW M 1 1 1 1 1 1 1 It  ‘ 1 1 1 n 11  It  ‘ 1 1 1 n n ti  1 1 n 

I n } 1 1 n ' > n i n M n n M < ] > 1 2 i M t ff)  t M n 1 1 n M n } M M M 1 1 1 i M t M It  i 

n n ) I n n n I n n 2 1 n 1 2 2 n M n i ] 7?)  1 1 M i n > ' I n 1 1 1 n I * n 1 1 1 1 1 M 1 1 1 

nj 

4 4 4 M4«44MM  4 4 4 Mt<<M4MI  4 4 4 ljn4m4y|lM4^4M44ltM  4 4 4 M4M4MMMIM<<M 


M 1 1 } M M M M M M i M i M M 1 1 M I U2 1 M i 1 1 1 1 1 *1 1 M i 1 1 1 2 ’l  I II I M ill  I 


Miiiiiiiiiiiiiiiiiiiiiiiiiiiii^i  i|(^l|k^VMt^niiiiiiiiiiiiiiiii(itiiiiiiti<i 

n M I MI  n I M n 1 2 1 M 1 n 2 2 2 2 2 2 2 2 2 2^1  ijy  2J  M rtTl  2 2 2 2 2 2 2 2 2 M M M 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 

1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ;:i  1 1 1 1 1 1 1 iZi  1 1 1 1 1 1 1 1 1 1 1 

IIMIIIillilMllllllillllllllltMMIIIIlTIMririllllllllllllllllllllllllllltfn 

I ) 4 1 1 > 1 1 n 2i9um)it  ii  «iiNtinJM»Nita«»ii»>iM»«}i«»««4f«fi«tftf««ifiittui«»NNi}«ii«v«<i««v««Niiitaaiiaif auM 
IMMCTT  m999 


I I 


Figure  6-2(c).  Data  cards  fw  the 

with  - 1 V/m,  0-0  , 0 -180  , -90  . 
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(a)  E » 1 V/m,  - 30,  0 » 150,  (}.  = 40 
n E>  p p 


f ■ 1 MHz 


1^(0)  - 3.298E-6  /89,41  I^«)  - 2.837E.7  /86.22 
12(0)  - 7.336E-7/  88.68  l2(';0  - 1.782E-7/-91.58 


f « 10  MHz 


Ij^(O)  » 3.315E-5  / 84.07  - 3.116E-6  / 53,64 

12(0)  - 7.191E-6/  76,96  IjU)  - 1.732E-6/-105.55 


f « 100  MHz 


Ij^(O)  • 2,495E-4  /-1, 650  ■ 1.024E--4/-142.78 

12(0)  - 3.450E-5  / 4,802  l2(x)  • 1.101E-5/-177.51 


f « 1 GHz 


Ij^(O)  • 2,089E-4  / 3,521  Ij^(a0  - 9.315E-5/<-139.76 
12(0)  - 3.317E-5/-IO.474  l2(3C)  - 1.089E-5  /l72.48 


Figure  6-3.  The  problem  in  Figure  6-1. 


I 

I 


’'<>  ’^■3- 


(b)  E • 1 V/m,  0_  » 0,  0^  - 90, 
in  £■  p 

f « 1 MHz 

Ij^(O)  ■ 9.294E-6  /89.09 
I^W)  - 1.963E-6  / 88.44 

f " 10  MHz 

1^(0)  - 9.316E-5  /80.85 
12(0)  - 1.920E~5  / 74.56 

f ■ 100  MHz 

Ij^(O)  - 4.638E«4/-37.08 
12(0)  - 6.602E-5/-24.28 

f - 1 GHz 

Ij^(O)  - 4.S87E-4/-37.91 
12(0)  - 6.567E-5/-24.92 


(|)  = 90 
P 

- 2.333E-6  /87.87 
l2«)  - 1.432E-7/-93.46 

- 2.336E-5  /68.63 
l2(r)  - 1.3838-6^-124.51 

• 1.150E-4/-156.86 
l2«)  - 3.021E-6  / 70.15 

I^(/)  - 1.138E-4/-158.43 
I2OO  - 3.054E-6  /68.47 


I 


f 


I 


Figure  6-3.  The  problem  in  Figure  6-1. 


(c) 


! - 1 V/m,  0_  - 0,  e - 180,  ♦ - 90 

Ml  D U 


f - 1 MHz 


lj(0)  - 3.A94E-6  / 90.08 
12(0)  - 5.590E-7  / 89.95 


lAt)  - 3.493E-6  /89.27 

X Lrnrnmmmmmm 

- 5.589E.7/ 89.44 


f - 10  MHz 


1^(0)  - 3.553E-5  /9O.71 
12(0)  - 5.656E-6  / 89.41 


Ij(f)  - 3.500E-5  782.65 
l2«)  • 5.581E-6/84.45 


f ■ 100  MHz 


1^(0)  • 5.316E-4  /33.83 
12(0)  - 8,392E-5  /52.8O 


• 1.988E-4  /■.6.817 
12(0  • 4.634E-5/  35,77 


f " 1 GHz 


IjCO)  - 4,402E-4  / 33.09 
12(0)  - 8.585E-5  /52.98 


Ijvt)  - 1.632E-4  / -7.429 
l2(X)  • 4.664B-5  /37.48 


Figure  6-3.  The  problem  in  Figure  6-1. 
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Wire 

Radius 

h 

0 

30  mils 

0 

0 

1 

30  mils 

0 

10  cm 

2 

10  mils 

4 cm 

7 cm 

3 

10  mils 

4 cm 

3 cm 

I 


The  four  frequencies  of  excitation  for  the  ground  plane  problem  In 
the  previous  section  (IMIlz,  lOMHz,  lOOMHz,  IGHz)  as  well  as  the  three 
orientations  of  the  plane  wave  will  be  considered  here.  Note  here  that  the 
zero  phase  reference  for  the  plane  wave  Is  not  the  same  as  for  the  ground 
plane  problem.  Here  the  zero  phase  reference  Is  displaced  downward  (in 
the  -y  direction)  from  the  zero  phase  reference  for  the  ground  plane  problem 
in  the  previous  section  by  3 cm.  This  means  that  the  phase  angles  of  the 
currents  In  this  problem  will  differ  from  the  phase  angles  of  the  corresponding 
currents  In  the  ground  plane  example  by  k(5  cm)  degrees  or 


cm}  (degrees j 


1 MHz 
10  MHz 
100  MHz 
1 GHz 


.0600 

.6004 

6.0042 

60.0415 


The  next  problem  remaining  is  to  determine  the  appropriate  representation 
of  the  terminal  networks.  This  type  of  situation  was  considered  In  Section 
2.6  of  Chapter  II,  From  Figure  6-4  we  may  write  (note  that  the  line 
voltages  are  with  respect  to  the  reference  wire  here) 
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I 


1^(0)  - -(1/200)  V^(0) 

I^CO)  - (1/1K)(V3(0)  - V,(0)) 

1^(0)  - (1/1K)(V2(0)  - VjCO)) 

I^(t)  • (1/lK)  V^et) 

- (1/2K)(V2(:^)  - V^W)) 

I^Cgf)  - (1/210  - V^CX)) 

Thus  we  select  the  load  structure  option  LSO  ■ 22  and  the  terminal  admittance 
matrices  become 


5E-3 

0 0 

’lE-3 

0 0 

0 

lE-3  -lE-3 

" 

0 

5E-4  -5E-4 

0 

-lE-3  lE-3 

-t 

. 0 

-5E-4  5E-4, 

The  input  data  cards  are  shown  in  Figure  6-5.  The  results  are  shown  in 
Figure  6-6. 

Note  that  for  all  angles  of  incidence  the  magnitudes  of  I^CO)  and  1^(0) 
for  each  frequency  are  equal  as  are  the  magnitudes  of  and 

Further  note  that  12(0)  and  I^CO)  are  precisely  180*  out  of  phase  as  are 
1 (t)  and  I (D.  Therefore  we  have 

m W 

12(0)  + 13(0)  - 0 

l2<^  + 13(0  - 0 

for  all  frequencies  as  they  should  be. 
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t • ) « 1 .^•*^?i*!inMBiiii'*flMit»nM»)iiiii«iinBa>*n»»»»«««««<'*i««***«««Httii««»»iiiiu»i»«»«»m>'n»**'n»n'<»» 

1.C8 

I I 1 4 j 4^*«^?il«OHB«il«B»n»n}«nMtlNI*a»ttaM»N)ia»NUtt4}4li4NII«llttll|;UHttMPH»HII«liMHIIII«MannnNnMII«aN 

t t ) 4 1 4^ It « OH « ii n hbh D ;in:t w»i]aun»}t»»MN4j«}i4ii«ttaitiitittt]Maa»r»M««««««a«iM«an anNaaitsait 
I ; j 4 i 0 0 0 H « « I)  0 jitTMaap««M«Hif«aaiia»Naaaaaii 

I I 1 o^^^*^ntQH»o9%^ipitiitinMftatt»SNa»aairaa««a«Nl^%f^«aua»»  f»ll 

<i4« m Mi  HUMMatt  m 

m m itauaaa  m 

1 4 Hunwaiif^KvaaMaNNMMaR aaM«Nff aaii 

I j 1 « if >1  •oHoof^fijl It ■»«»■!!««« »»«»»» »»a4iHtt«M«9l^(^aNNaa»  Ml 

..,..M-;i. Ml  UU  auuMsa  Ml 

I I 1 < 1 « I I I *lllf  l.ll  .....»M.lll7l..nM......ll(Nl...... .................... ...... (.g 

1 1 1 1 > 1 1 1 ) « I'lill. Ii5»l..........?i5^l... ............................... u 
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Figure  6-5(a).  Data  cards  for  the  problem  in  Figure  6-4 

with  E -1  V/m,  0^-30  , 0 «150  , <|i  -40  . 
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Figure  6-5 (b).  Data  cards  for  the  groblem  in  Figure  6-4 
with  Ej^»l  V/m,  ©g-0  , 0p-9O  , (}>p-90  . 
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Figure  6-5(0.  Data  cards  for  the  problem  Ig  Figure^fi"^ 
with  E^*=l  V/m,  0j.-0  , 0p-l8O  , <j>p»90  . 
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= 1 V/m,  0_  = 30,  0 = 150,  <|)  = 40 
m E p p 


f = 1 MHz 

1,(0)  » 1.649E-6  /89.46 

(f)  = 1.419E-7 

/86.28 

1,(0)  ='  3.668E-7^  88.79 

l2<t> 

» 8.912E-8  ^ 

"91.43 

1,(0)  » 3.668E-7/-91.21 

l3(^) 

= 8.912E-8 

/ 88.57 

f - 10  MHz 

1,(0)  - 1.658E-5  / 84.57 

13(2’) 

- 1.558E-6 

/ 54.19 

1,(0)  - 3.597E-6  1 78.11 

l2(/) 

- 8.661E-7^ 

■104.07 

1,(0)  - 3.597E-6/-101.89 

i3(;r) 

8.661E-7  ^ 

/ 75.93 

f ■ 100  MHz 

1,(0)  - 1.247E-4  / 3.519 

\(X) 

- 5.II8E-5/-I37.55 

lj(0)  - 1.763E-5  / 15.05 

IM) 

- 5.657E-6/-165.87 

1,(0)  • 1.76eE-5/-164.95 

- 5.657E-6  ^ 

1 14.13 

f - 1 GHz 

1,(0)  - 1.037E-4  /55.36 

» 4.631E-5  [ 

-87.41 

1,(0)  - 2.591E-5  ^74.23 

l^iX) 

* l,024E-5/^ 

-93.i;o 

1,(0)  - 2.591E-5/-105.77 

l^it) 

- 1.024E-5  ^ 

' 86.50 

Figure  6-6.  The  problem  in  Figure  6-4, 
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6.1.3  Comparison  of  the  Two  Solutions 

The  terminal  currents  (1^^,  I^)  in  Figure  6-1  (the  ground,  plane  problem) 
should  be  exactly  twice  the  magnitude  of  the  corresponding  currents  (Ij^,l2) 
in  Figure  6-4  (the  image  problem).  For  the  case  of  propagation  in  the  +x 
direction  with  E in  the  +y  direction: 


(b)  0g  = 0,  0p  » 90,  = 90 


-jkx 


E = e 


and  the  case  of  propagation  in  the  -y  direction  with  E in  the  +x  direction: 


(c)  0 “ 0,  0 ■ 180,  <|)  « 90 

E p P 

t - t 


this  is  precisely  the  case  for  all  frequencies.  (Although  not  shown  here, 
thn  results  were  printed  out  to  16  digits  and  agreed  to  15  digits.) 
However,  consider  the  case  of 


(a)  0„  * 30,  0 - 150,  (J>  » 40 

E p p 


Note  that  the  corresponding  currents  for  the  ground  plane  problem  in 
Figure  6-1  are  precisely  twice  those  for  the  image  proolem  in  Figure  6-4 
for  1 MHz  and  10  MHz.  For  100  MHz,  and  Ij^(X)  correspond  exactly  and 

I„(0)  and  I„(J’)  correspond  very  closely.  However,  for  1 GHz,  only  currents 

4 4. 

Ij(0)  and  correspond  whereas  12(0)  and  l2(^»  although  within  a factor 

of  two,  do  not  correspond.  The  reason  for  this  becomes  clear  when  we  con- 
sider the  definition  of  line  voltages  used  for  the  two  problems. 


•a 
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Consider  Figure  6-7.  We  have  shown  the  cross-section  of  particular  wire 
above  a ground  plane  in  Figure  6-7 (a),  the  wire  voltage,  V,  is  shown  and  the 
potential  difference  between  the  wire  and  its  image  is  2V.  In  Figure  6-7 (b), 
we  have  shown  the  corresponding  image  problem  and  the  voltage  of  each  wire  is 
defined  with  respect  to  the  reference  wire,  i.e.,  and  V2.  For  the  tv;o 
representations  to  yield  corresponding  results,  one  would  expect  that 


2V  i (V^  - V^) 

These  voltages  are  related  to  the  integral  of  in  the  transverse  (y,z) 

plane  along  the  contours  shown  in  Figure  6-7  and  are  included  in  the  vectors 

E^inc)(0)  and  Clearly  these  will  correspond  only  if  is 

curl  free  in  the  transverse  (y,z)  plane,  i.e,,  only  if  there  is  no  component 

of  H in  the  x direction  which  penetrates  a transverse  contour.  For 

angles  of  Incidence  0_  ■ 0,  6 • 90,  ^ ■ 90  and  0_  « 0,  0 » 180,  <))  » 90, 

® E p p E p ’p’ 

this  is  clearly  the  case  and  the  results  show  this.  However,  for  0^  ■ 30, 

0 ■ 150,  ^ • AO,  there  is  a component  of  in  the  x direction.  How- 

P P 

ever,  for  f ■ IMHz,  10  MHz  100  MHz,  the  cross-sectional  dimensions  of  the 
line  are  electrically  small  and  the  fact  that  is  not  curl  free  in  the 

y,z  plane  does  not  matter.  For  1 GHz,  the  cross-sectional  dimensions  of  the 
line  are  not  electrically  small  and  it  does  matter  as  is  evidenced  in  the 
computed  results. 

6.2  Example  II 

In  this  Section  we  will  consider  a problem  which  was  investigated  by 

Harrison  using  an  alternate  formulation  [5,8],  The  problem  consists  of 

-3 

three  wires  in  free  space  all  of  radius  10  » which  lie  in  the  x,y  plane  with 

-2 

adjacent  wire  separations  of  10  mas  shown  in  Figure  6-8,  The  line  is  10  m 
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long  and  various  loads  connect  each  wire  to  a central  point. 
The  termination  impedances  are 


= 50  - j 25  Zq2  = 100  + j 100  » 25  + j 25 


= 50  + j 25  ^2  “ ^00  ■ i 50  2^3  * 150  - j 


50 


Obviously  we  should  select  LSO  = 12  and  the  termination  impedance  matrices 
become 


\l  * ^03> 

^03 

B 

’ 75  + jO 

25  + j 25  ' 

- ^03 

<^02  ^ ^03>- 

25  + j 25 

125  + j 125_ 

2^  • 

* w 

■ 

m 

“ 200  - j 25 

150  - j 50  ' 

1 

««  + ^X3  1 

_ 150  - j 50 

250  - j 100_ 

AlV/m  uniform  plane  wave  is  propagating  in  the  +y  direction  with  E in 

the  +x  direction,  i.e.,  « e x . Therefore  £»  ■ 1,  &_  ■ 180, 

w*  E 

0 ■ 0,  (|i  =90.  Harrison  showed  the  result  for  the  terminal  currents  com- 
P P 

puted  by  his  method  for  k;;(  » 1.5,  The  line  is  10  m long.  Therefore  the 
frequency  is  7157018,74  Hz.  The  input  data  cards  are  shown  in  Figure  6-9. 

The  computed  results  are  shown  In  Figure  6-10  and  compared  with  those 
obtained  by  Harrison.  Note  that  the  results  computed  by  this  method  agree 
with  those  computed  by  Harrison  to  within  three  digits.  The  main  reason  that 
the  results  do  not  agree  precisely  is  that  the  ratio  of  line  length  to  wave- 
length at  this  frequency  is  .239.  Thus  we  are  entering  a frequency  range 
where  the  variation  of  line  responses  with  frequency  is  generally  quite  rapid 
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Figure  6-9.  Data  cards  for  the  problem  In  figure  6-8 

with  E^-1  V/m,0£»180  , 0p-O  , <t»p-90  , FSO-1. 
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-f 

I 

i 

I 


E = 1 V/ra,  0_  » 180,  0^  - 0,  ♦ » 90 

m a p _ p 


f n 7357018.74  IIz  (kX  » 1.5) 

Computed  with  WIRE; 

1^(0)  - 1.066E>5  ^->99,83  , I^(af)  • 

I^CO)  - 5,647E-5/-159.07  l2«)  - 

Computed  by  Harrison; 

|I^(0)|  - 1.065E-5  ll^Wl 

|l2(0)|  - 5.64AE-5  |l2(i;)| 


t 


! 


( 

1 

! 


1.221E--5  / 158.65 
2.784E"5Ai48.26 


1.220E-5 

2.784E-5 


Figure  6-10.  The  problem  in  Figure  6-8  with  FSO  ■ 1 


and  any  seemingly  insignificant  approximations  in  the  input  data 
(k3?:  1.5)  can  cause  significant  chang::s  in  the  result.  In  fact  the  program 
was  internally  modified  (temporarily)  such  that  k;fwas  precisely  1.5  and 
the  results  agreed  exactly. 


One  additional  case  will  be  computed  in  which  E = 1,  0_  ® 0,  0 « 90, 

map 


<l>  “ 90,  i.e.,  the  wave  is  propagating  in  the  +x  direction  with  E in  the  plane 
P 


of  the  wires,  i.e.,  the  +y ’direction. 


■+(lnc)  -jkx  -+• 
E - e y 


The  input  data  cards  are  shown  (for  k;;(  ■ 1.5)  in  Figure  6-11  and  the  computed 
results  are  shown  in  Figure  6-12. 


6,2,1  Use  of  the  Wonunlform  Field  Specification  Option.  FSO  » 2 

In  this  Section  we  will  solve  the  two  problems  considered  in  the 
previous  Section  by  using  the  nonuniform  field  specification  option. 
For  the  first  example  we  consider  the  problem  in  Figure  6-8  with 


E ■ 1 V/ra,  0^  ■ 180,  0 ■ 0,  ■ 90,  To  use  the  nonuniform  field 

m E p p 


specification  option,  we  must  describe  the  magnitude  and  phase  of  the- 
incident  elactric  field  along  the  three  wires  and  along  straight  line 
contours  (the  y axis  in  this  case)  between  each  of  the  two  wires  and  the 
reference  wire  at  x 0 and  x =jf.  Because  of  this  particular  field  orienta- 
tion, the  specification  of  these  quantities  is  quite  simple.  Clearly,  the 


transverse  fields  are  zero.  The  longitudinal  field  at  all  points  along  the 

reference  wire  are  l/o;  along  wire  1 are  x 10  ^ * 1 /-  8.5944E-2 

and  along  wire  2 are  l/v.  x 2 x 10~^  = 1 /-  1.7189F.--1.  Although  redundant,  11 


specification  points  for  the  longitudinal  fields  and  6 specification  points 
for  the  transverse  fields  will  be  used.  The  data  cards  are  shown  in  Figure 
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i % 
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Figure  6-12.  The  problem  of  Figure  6-5  with  FSO  » 1, 
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6-13  and  the  computed  results  are  shown  In  Figure  6-14.  These  results  com- 
pare (and  should  compare):  c::actly  to  those  in  Figure  6-10  where  the  uniform 
plane  wave  option  was  used. 


The  next  problem  is  to  use  the  field  orientation  of  9^  *■  0»  9p  ■ 90, 

* 90.  In  this  case,  the  longitudinal  fields  along  all  wires  will  be 
zero  whereas  the  transverse  fields  at  all  points  along  the  contours  at  x «•  0 
will  be  l/o  whereas  those  of  x ■ ^ will  be  1 85.943669.  The  data 

cards  are  shown  in  Figure  6-15  and  the  computed  results  are  shown  in  Figure 

t 

6-16.  These  results  compare  (and  should  compare)  exactly  to  those  In 
Figure  6-12  where  the  uniform  plane  wave  option  was  used. 
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t(l)K(aaaaaaaaaaaaaaaavaaaiiaaananaaa 
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1 n a aim  a.ifinliftiKliBaum 
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1 1 1 « 
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BS4Kj|i4uaMB8iiaBBBe8Baai)BBai)aaaa.'4naaB 

1 ) ) 4 t^lt^Untottagtaataaaal^KSaaaa: 

449^KS)iaaBaBaa88BasBaa4)aait)4«aaaaaaaa 



a84Ki)iaaMBaBaa8BB8Baa«a88iianaaa)iaaa 

.).4.?1«S 

1114  tl^f^{na«aaaaaa}ianni4a(frjK8nN»M 

««)^KlMauH»»Ptt»aa««MnMii«N»a  aaaaa  n aan 
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. . . 

mfVMimiiuMMiiwiisiiHrDMiiiiiiiiiiiiimvMiiiimma 

MU  a It  mm  mi*  Ul;$  Ntt«imaumttiaPii»fia«»Mii«i)ttiiai:aaNaaaaaii 

anaaiiaanfta«ttaaiiaifaa»imaii«i««mM««a«i!»NwuHaawNMiiatf«MHiiimMi»RaaaiianaaN 
1 ) 1 4 4 4 1 4 1 In  to  ••i)tlfi)taaa«»«a»$))a«»»«»»»i  444>4)««l(KtiiuuH«it>KtK««H«ii))Klli)  »»«»«))»»■ 


MMIIMItMUIMMIIIliUllllMItlDHIIIHIM  'I  'lltIurrillllH  'IIIIIKIIMI 

i))4iii)4»iiti)atat  t •N)i)))))4»a)'aa»«aa»«>^)i(M^)fA)44M44iii«w»HUM)i»4>a»»««))Naa)i)iH«>i « i)i4at  n tat 

1 1 1 1 1 1 n r 1 1 1 1 111  I * 1 1 1 111  1 1 1 1 11 1 1 i-^^iV  i>f  ni/*  1 1 1 1 1 1 1 1 1 ‘ 1 III  1 1 1 1 1 ‘ 1 1 1 ‘ 1 1 1 1 1 1 1 1 1 1 

iiininimtmm’nmnn'n  nmnnnittmiminmttn 

in  n ) n 1 n I ) 1 1 ) ] ] 1 ] 1 1 ) ) ni  n n i i n n n ' i n n n n ' n n n n n i n 


ittnnuutnnuMniiuiMnl 


114  •'  - Af 

I tn  < y « 4 (|<l  n 4 j ^ « U < < 4 < < t 4 4 4 4 4 <4  I « M 4 4 ( 4 « 4 < « 4 < 

i n i n i n i n n n n n n n n n f n U4  nl)  i'l  i } inO  4 1 f n n n I n * n n n n i ’ n n n i n n 

1 1 • 1 1 1 M 1 1 i 1 1 i i 1 1 1 1 1 1 1 i i 1 4 1 1 1 1 1 1 <'-11 4 1 M 1 1 1 '«  i I!  a<n  M i 1 1 M i I i i I S 1 1 1 1 1 1 1 1 1 1 1 i i I i m 

" V 

1 1 11  n 1 1 1 n n M 1 1 n n n 1 1 M I n I mf  I i,nn.;  n 1 M 1 1 11  n 1 11 M n n I n 1 1?  M I n 1 1 II 1 1 1 

/ U ■* 


1 1 1 1 1 1 r.  1 1 1 1 1 f 1 1 1 1 1 1 1 1 1 1 1 1 1 


V B U 1 1.1 1 II  nj  1 1 1 !;  1 1 1 1 1 1 1 1 1 ' 1 1 1 1 1 M 1 1 1 1 1 1 


PUTING  CENT  r.R 
I'nniiiiiiiFsnnis! 


1 1 n M I IS  M n I n n IS  n n 1 1 n n Hi ; 1 1 li  f 1 1 1 1 n i s n i s i s 1 1 1 s 1 1 1 1 1 1 s 1 1 1 1 1 1 1 1 1 1 1 1 1 s s s 

I I ) 4 1 I I I «iaii9ii'm<'>««Nrinn}4»Nra»a«an.«t)«a»aa4it:«4r««t««i«4U*tMMus>4»aiiM:ii»»><'»t».Nm;MM  }ii7.*)waaT>ii« 

mtmrf  mtom 


Figure  6-13.  Data  cards  for  the  problem  in  Figure  6-8 
with  E - 1 V/m,  0£»18O  , 0-0  , 4-90  , 
FSO-2  ^continued) . 
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) . 


I 


SI » n H n fi  sr » » • • • ■ « « « nil  a I*  n « n N « « fl  II N N 

I 7 1 « w*^i«ii<roHtf«n«itN7ti;nMn  M » » » » HU  M M a fl  n a « « mt 

I  7 j 4 l^tl•aNfts^l■ll7S7lnIIMn  hm  »»»  HitflisaaiiaaMa  l»KB  «»9  •smaaaftaMMaNaaaanvMNaaflaNN 

I 7 ) 4 j^^^twiivoMnin « 7i»7)ttaanil^li^i2ftanii»NVMaa««a«  

11*4  f^TiirsutiiintinaTianMM  8»8S9  K nauaNMiieaM#  

7 1 4 1 I I I i?ii0aMn«niiaan»n»n  aauaHaitaaN*  »KS  aaitaaHaaaHaatiaaaaaaaaaaaaa 

q M 17  « na 71  u TiNa  aavaa  hKi  »nv»umvmmmmm«mmnnnnnnniinnm 

aanaqaannaa  aapMaaaaaMaavwaaRaaaaanaaa 

...M7?  d « fi « q a 11  tt  11  }4  a aa»aaa4i«i4}M«  sK8  atMtaapaavaaaaaaaaaafiaaaaaaaaa 

d d n ddaTnaa  m a a i>  a a ...-.tiW 

dddddananaa  tsKI 

1 1 4 1 « I I i)itdOHddiiddafinnNilt%fC}«aaaaapaa«4iaa««I^K(aaaaaaPaaa«aaMaaa««aiiaaaaaaaaa 


aaaaa: 

1 

aaaa 

maaai 

lanaa* 

lavsaa 

avaa 

laaaai 

ImtaMa 

laaa* 

laaaai 

*aaaa«: 

taaaMa 

aiiad 

saapam 

tnunas 

('aval 

lanaai 

waaaaa 

(attaaii 

taaai 

laaaai 

daaaaa 

aaa aa» 

.6.«l .UKf..... 

• v a a MalsItlK  •{) » 
» u a a a altlttHHal  a 

aaaaaltsitlllat. 

t4.t» 





ftddariana  alliKt  Maaaaaaaa  a dtl  t ati  n 

UKi 



i1^^?9iid0ddi 

t If  d d a a a n a al  t a a a a a a P a a al*^  t 

,^,9^9,, 



2^^  11  V 0 H d I 

1 n d dan  ananl^Ktaaaan 

l.tiWSn,.... 

1 If  d ddfinnanlif  fClaaaaaaaaaanf 

1 1 1 3 il^itnaddd'iddanii nanli^^Siiaaai 

1 a P a a 

MMUMBawaaMaaaaaaaaaaiiaaaaaaaaa 


I8IM  till  HIM  IIIIIIIMliri 'MliniMIJIIItllMMIMIIMIIMIMIMIMMIIIIIM 

iijiii’iidiidQMddn  d da  nanafiaaaaaad  a -sAmIh^  iHiwaaaiiMaaaBPaaMaaaaeaaaaaii  Baa  a*  a aav 

1 n It  I M ‘ 1 1 1 1 1 1 1 n 1 1 1 H n M 1 1 1 1 1 1 iVf  n i^i  n n 1 1 n it  1 1 n 1 1 n ii  1 1 1 j n i ii  1 1 1 1 
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]]  1 3 3 3 1 3 1 1 : 1 3 3 3 3 ] 2 2 2 2 2 2 3 3 2 ' 1 1 2 2 1 ! 2 ^ 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 2 3 2 1 2 2 1 1 1 1 1 1 1 2 

4 4 4 < 4 4 4 4 4 4 4 4 « 4 4 < 4 4 4 4 4 4 4 4 4 4 4 4 4 4 4 b'l  < l(*  f ^ ^4  < M 4 4 4 4 4 4 I « 4 U I 4 4 4 4 4 4 4 1 1 1 4 4 M 4 4 

2 3 S 3 2 3 2 2 3 3 3 3 3 2 3 S 3 2 2 2 2 1 2 2 2 2 2 2 * 2 2 2 i(2|{*i  < 2 2 2 2 2 2 2 2 2 2 2 2 2 2 1 2 2 1 2 2 2 2 2 1 2 2 2 2 2 1 2 

--y  « ! -4^ 
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3 I 1 ?)  M 3 ) m I 3 I ) 3 3 ? 1 3 I 3 3 3 3 1 3 3 2 3 3 3 Ml  3 2 3,3  37  H 3 3'?3  3 3 3 3 3 2 3 1 n 3 3 3 3 3 3 3 2 1 2 1 2 2 3 3 3 1 3 1 3 2 3 

I i S 1 1 1 1 1 3 1 1 1 1 1 4 1 1 1 4 2 4 4 I 4 4 4 ^ il  j^l  I m 1 1^1 1 4 4 1 1 1 1 1 1 1 1 1 1 2 1 1 1 1 1 1 1 1 1 1 • 1 1 

3HS1JMS!1!JM3H24l!n!24‘2r4m:S  4 4 4 4 4 4'l24'  M 0 4 24*4  4 2 4 4 4 4 4 4 4 4 4 4 1 4 1 1 1 1 4 4 4 2 4 4 4 4 


4 4 4 4 4 4 4 1 4^4  2 4 4 4 4 4 4 4 4 4 1 4 4 4 4 1 4 4 1 4 1 4 4 4 1 4 4 
4^  24244  n'  4 4 4 4 4 4 4 4 4 4 4 4 4 1 4 1 1 1 1 4 4 4 2 4 4 4 4 4 


0)919  M19)nM3Hin!ni!4434SS  4:  9 4 4 4 4 4 4 4 2 4 4 4 4 4 4 1 4 4 4 1 4 4 4 4 4 4 4 4 4 4 1 4 4 1 1 1 4 4 4 2 4 4 4 4 4 

I 7 4 4 2 • I I 2 ' }»  » 4-*a.«'7<  u/'.rfiiaaailUo.daaiiadaiitiUMddiiaaaKuuaaavaaMiattttwaaaaMaiianaaddaaM 

N.;uir  20K4» 


Figure  6-13.  Data  cards  for  the  problem  in  Figure^6-8 
with  E =1  V/m,  0-180°,  0-0.  4-90  , F£ 

m t P I* 


*“»  ''E 
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0-2  (cent.) 


I 


2.E-2 

t I I 4 I • » I •« 


noaM»Mir«»ii}tnDMn 


i I *11  *T«l|ntat«»»ii«iii»iionMli 


8.E-<3 

t } ) « } • l7l  IH 

4.E>3 

I } ) « } « >7  « « n 


1 0 H « « R « 4 M nn  » I*  n 


hm 

hm 


k s a M n » u » » «•  « « « 41  «i 


2i,K2 


«•«««««  SI  MnancaMniiffaifM 


• OH««o««»}innisiii 


aMMUDNuan  nnMAMA  anit 
{trtSSkunisasvfiMitnvaMiiMHiiMafi  vnMff  air  Nftn 

UsSSii......... 
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naan  a aa  9ti»nMaap»»0lHi%?MKakt 
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1 I 1 4 1 1 1 «^9i4faH»aita«nattnHatf9Cti 



MttHts«a>iK»t»Mn0iiaMiiN»aMaaaaanaaM 

tittuisaaPM»0MttO0«M00«aii«aaaaffaaN 
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"•w  I'i 
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4 4 4 4 4 n n 4 n 4 4 I ( « 4 4 4 4 4 O 4 4 4 4 4 4 4 I 4fl  ^ « 4^  4 444  < 4 M 4 C 4 4 4 I 4 44  4 4 ( 4 4 4 M 4 4 4 4 ( « 4 4 4 4 
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IISIIIIII'  4ISII4IIIIIIIISIS'  I ' IIIISIIIMJIIIKMIlllllllllllllllllSSIIMIIISlISS 
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sn-OA^PuTiNG  CENitR 

tSSSSSSISIISSf  SS4IS4IS44ilSI4ISII!SilSSIIII44r  4 4 SSI  4 1 4 1 4 4 4 4 4 4 4 4 1'S  4 1 4 4 4 1 1 4 4 S 1 4 4 4 

4 I 4 «s  t I t t iiiatju3<*u«aai'naa»ara«0iia  .4»aH»a««tf«4s««va4i’s  iv«ivsaa)iiai.ti"utiMOsiHM«iii}aMn'i>i'ia« 

MMtTf 

Figure  6-13.  Data  cards  for  the  prgbl^  ig  6-8 

with  E^-1  V/m,  0J.-18O  , 0^-0  , (J»p*'90  , FS0»2. 
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E = 1,  0„  - 180,  0 «=  0,  (fr  =90 
m * E * p * 

f = 7157018.74  Hz  (k;t  » 1.5) 

Ij^(O)  = 1.066E-5  /-99.83  I^(^)  » 1.221E~5 /l58.65 

IjCO)  - 5.6A7E-5/-159.07  - 2.784E-5/~148.26 


Figure  6-14.  The  problem  of  Figure  6-8  using 

the  nonuniform  field  specification  option. 
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,aaa9l.959,. 

.2iK9 





AM 

VW9.., 

AM 

5-.W9., 

bUKS 

.ti^Kt 
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MSMVSVSHSUSSSftSSIt 

\9. 



..9.;5e8 
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2S.C0 

4 1 1 4 » 4 » »* 
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n 1 > > n ) ) n n j n ] M n ] 1 n 1 M n } iii  i M'lV-i'/i  nii  j'tfi  i n i n i n ) ) n n > ' i m n n i n 1 1 1 

,1.;  ; 

M 4 M 1 4 I M 4 4 4 4 4 4 M 4 ( M M M M M 4 M M ( 4 M (,(  I 4 M I 444  1 1 1 4 M 4 I « i M ( M < 4 4 M M 4 « 4 M 4 4 4 4 

> n / *-- 

M ) s M s } M n n ) ) M M M M n n n ) 44.4 ) n vi  i i i n > t ) i > s } m ! s s t > > > i i s } i s > i i j 1 1 4 1 

• r ‘ ^ 

f i I 4 4 4 i I 1 1 1 ( I 4 1 4 • 4 4 4 4 4 4 4 ! ! ] ( s t S ( 4 1 1 ( 1 1 1 1 4.4  4 4 1 t:iF|-4  4 4 I 4 4 4 4 4 4 I 4 4 I 4 4 I 4 I 4 I 4 4 4 4 4 4 4 4 44  4 

•*^1'^  ” -'Am 

I t ( I M M I M M M M M M M M M n M M M I y y I J M M\|  M » M M M M I M M M M n M n M J n n 

liMH4ll4lllll4lll4ll4IIII^^MI^J4l 
M : 1 4 4 M 4 14  M 4 4 1 4 4 1 1 4 4 4 4 1 4 4 4 4 4 4 4'4 


M ^4  4 4'4  4 4 ! 4 4 4 4 4 4 4 4 e ! 4 4 4 4 1 1 4 4 4 ! 4 1 4 4 4 4 4 4 4 4 4 4 4 


I444l44444t4t4t44444:444i44l 
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Figure  6-15. 


Data  cards  for  the  oroblem  in  Figure  6-8 

with  E -1  V/m,  0^-0°,  0-90°,  <j)  -90  , FSO-2  (cont.) 
in  D P p 
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l.E-a  l.OEO  -61 

I 4 i i I I nwnnniiuuuM  nMiYMii4t4itt 


8.E-3  l.OEO  -&S.  943669 

I 2 1 4 } ( > I • « ti  0 0 M B 9 H » II It 21  s nMnMiliinM H nu)(i>)iuMJi«o«4»ow««««ioM n n nN n n n MWM 


I ; ) 4 l«tiDQM»9  II  BHItTi  II  nHn  b S K n D 1.  > ■ » 1. 1. 


4.t-3  i.0£0  *85 

ii)4>(>ii9iiiioMB9ii9  Bit  n n»}«n2iniiait«»)iM>i»)*jt»«t4i«itt 


2.£««3  U9E0 

f I I 4 1 I I $ i«iiouHBiiii«iiititii DunMliiilfM r Da>iiiit»)tjt«t«i 


I I 4 } t I I |9il0OHA9ll«T}NIinn}in9flMllN»ffUM»)t>IM  »4t « 


4.e-3  UKS 

I ) « 1 « I I *«HoaHB9naiiMii a a}« niiliiiWiiiia nnaaiiii 


bSSS» 

ibS^Sr. 

U9S9.. 

>b9^9ii 

>b9l9> 

i>b9$9> 

nb9^9]i 


BVMuwiiKaMititiiatttiaiiiiMitananMiiaaanii 


iia»M»»}>M»itiianMnwuititanaDNanaaNa 

MuuliNMMMStMiivaMBMnuiiatiaaaaMaanai 

vuuHimuHitittivuwBiiiTMNananNitaaNan 
9r9S9  iiu»»4aiiuMiiMii*i]iiNMn«iiMiinnNaaaRnit 
«4t4I«MW  9i,9$9  fluuHMiiuHitiiMaiiMaMUMMar.  anaxaffann 

9<,9I9  ttuuMaapMaaMnaMaMvai^aaai>Maaffaa« 

......9I.9S9  vunMaapaaMiieaHManaaaNaaaaanaaM 


I I 1 4 hm  n u n M IS  M p a PM  41  u tt  M s 9r9S9  PUBpaapaaaaapiiaMnaaaiinnaaaaaait 

8.C-*3  l.OCO  O.OEO 

I I ) 4 s I > I I « II 0 q « a « >1 « naa  ana  aan.TiiaiiuaMaMPaMa4ippttM8ii««9apuMaaPMaaPsaMB«paMar.  anaaaaaaa 

8.C*-3  ItOCO  e.OCD 

I I ) 4 1 I I t «itiioaMa«naiiaa a aa »anilaa»aaaaapaaa4ipp«4 ««uaa9ttyuMaaPMaaapaMaapai^aa anaaaaaaai 


iPUPMaapaaMppaMMapaaaaaaaaanaaii 


I 1 I w*wi9iiouMtiv  'aitaaanHnlliSSMapMMnapaaitiiQuMcSTiKSMuuMBaPMMnapaHBHpaBaNaDaaaaaaM 
iiiii’ii^iivnMBB  artananna  b9^9  )ianaaaPMaw4iao44M  9^9«9  MPUwaapttaaaeaMaapaaattitaNaaaaaM 
iidaMBNiianaiiaann  hm  iinuaaaPMaa4tttpw« 

«upwattPaaa«p8Haapa«aaaaaaanaaM 

9«.  (niMtianartananwa  9(9«9  naBMtiaiiii.'ap4iu4«.«  9»9» 

ttupMMapaaaHDBMaapaaaaaaaaaaaaM 

I I I 1 ^ * 9^  aiidDHBaiiaiiaiinnaa  9i,9S9  nuuaaMpaMa««aM«  9^9$9 

auuMaapaaanppMaapaaan  anaaanaaa 
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Figure  6-15.  Data  cards  for  the  pr6blem  in  Figure  6-8  with 

E =1  V/m,  0=0°,  0 =90°,  ({)  =90  , FSO=2  (cont.) 
m E p p 

-147- 


n H 


\ 
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Figure  6-15.  Data  cards  for  the  problem  in  Figure  6-8  with 
E =1  V/m,  0^=0  , 0 =90  , <t>  =90  , FS0=2. 
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VII.  Sum-IARY 


A digital  computer  program,  WIRE,  which  is  designed  to  compute  tht 
terminal  currents  induced  in  a multiconductor  transmission  line  by  a 
single  frequency,  incident  electromagnetic  field  has  been  described.  The 
transmission  line  consists  of  n wires  (cylindrical  conductors)  and  a 
reference  conductor.  The  reference  conductor  may  be  a wire  (TYPE=1) , an 
infinite  ground  plane  (TYPE=2)  or  an  overall,  cylindrical  shield  (TYPE=3) . 

All  (n+1)  conductors  are  assumed  to  be  perfect  conductors  and  the  surrounding 
medium  is  assumed  to  be  linear,  isotropic,  homogeneous  and  lossless.  The 
line  is  assumed  to  be  uniform  in  that  all  (n+1)  conductors  have  no  variation 
in  their  cross-sections  along  the  line  length  and  are  parallel  to  each 
other. 

Two  types  of  Incident  field  specification  are  provided  for.  Uniform 
plane  wave  excitation  can  be  specified  for  TYPE  1 and  TYPE  2 structures 
whereas  nonuniform  field  excitation  can  be  specified  for  all  structure 
types. 

The  primary  restrictions  on  the  program  validity  is  that  the  cross- 
sectional  dimensions  of  the  line,  e.g.,  wire  spacings,  must  be  electrically 
small  and  the  smallest  ratio  of  wire  separation  to  wire  radius  must  be 
larger  than  approximately  5. 

General  linear  termination  networks  are  provided  for  at  the  two  ends 
of  the  line. 
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APPEi^DIX  A 


A Note  on  Common  Mode  and  Differential  Mode 
Currents 


An  important  assumption  in  this  method  is  that  the  sum  of  the  currents 

in  all  (n+1)  conductors  at  a particular  x along  the  line  is  equal  to  zero. 

This  is  the  conventional  notion  of  transmission  line  currents.  The  purpose 

of  this  Appendix  is  to  provide  some  justification  for  the  assumption. 

As  a prelude,  consider  the  two  conductor  line  (n*l)  shown  in  Figure  A-l(a), 

At  a particular  longitudinal  coordinate,  x,  we  have  separated  the  total 

current  into  a common  mode  component,  I„,  and  a differential  mode  component, 

u 

Ip.  This  is  purely  a mathematical  operation  and  given  the  currents  Ij^(x) 
and  Iq(^)»  one  can  always  resolve  them  into  these  components  as  shown  by 
the  following.  We  are  simply  looking  for  a unique  transformation  which 
performs  this  separation.  If  we  write 


Ij^(x)  » lp(x)  + Ip(x) 

(A-la) 

l2(^)  “ - Id(x) 

(A-lb) 

then  in  matrix  form  the  equations  become 


r 

1 

1 

’ipC*)' 

X 

CM 

,J 

-1 

« 

1 , 

T 


(A-2) 


The  essential  question  here  is  whether  T is  nonsingular  which  would  represent 
a unique  transformation  between  the  two  sets  of  currents.  Clearly  T is 
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nonsingular  and  we  may  write  (multiplying  (A-2)  on  the  left  by  T 


Therefore,  given  1^  (x)  and  l2(x)  for  a particular  x,  one  can  uniquely 
determine  Ip(^)  and  I^(x). 

The  question  of  physical  significance  of  I and  I is  essentially  ir- 

D O 

relevant  since  this  is  merely  a transformation  of  variables.  The  essential 
point  is  that  as  far  as  the  terminal  responses  are  concerned,  we  need  only 
consider  the  differential  mode  (transmission  line  type)  current,  Ij^, 
since  the  common  mode  current  (commonly  called  antenna  type  current)  has 
essentially  no  effect  on  the  terminal  responses.  The  justification  for  this 
statement  lies  in  our  fundamental  assumption  that  the  cross-sectional 
dimensions  (wire  separation)  of  the  line  are  much  less  than  a wavelength. 
Tlierefore  we  may  consider  the  terminal  impedances  (Z^  and  Z^)  as  lumped  and 
if  we  apply  Kirclioff's  current  law  to  the  "nodes"  containing  the  impedance 
we  can  only  conclude  that  the  common  mode  current  is  zero  £t  the  endpoints 
of  the  line,  i.e,,  1 (0)  = I (t.)  = 0,  At  points  along  the  line,  this  is 
not  generally  true  and  the  line  currents  will  not  be  simply  due  to  the 
differential  mode  current  but  will  be  a combination  of  Ip^(>t)  and  Ip(x)  as 
shown  in  Figure  A-l(a).  The  essential  point  here  is  that  if  we  are  only 
interested  in  computing  the  terminal  response  of  the  line  (as  we  are  in  this 
report),  we  may  disregard  or  omit  consideration  of  the  common  mode  current. 

The  extension  of  this  result  to  multiconductor  lines  is  quite  similar. 
Consider  Figure  A--l(b)  where  we  have  decomposed  each  line  current  into  a 
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differential  mode  current,  I . , and  a common  mode  current,  I . Note  that 

Dl  vf 

we  have  taken  the  common  mode  currents  to  be  the  same  at  corresponding  points 
(values  of  the  x coordinate)  along  the  line  In  all  line  conductors.  The 
justification  for  this  Is  our  primary  assumption  that  the  maximum  cross- 
Sectional  dimension  of  the  line  is  "electrically  small",  i.e.,  much  less 
than  a wavelength.  The  essential  question  here  is  whether  we  can  define 
a unique  (nonsingular)  transformation  between  the  actual  line  currents, 

Iq,  , I^,  and  the  decomposition  components,  1^^^,  * 

This  is  easily  found  from  Figure  A-l(b)  from 


I (x)  - I-(x)  + I_  (x) 
n C Dn 

Ij(x)  » IqCx)  + 
lj(x)  ■ Iq^x)  + 

n 

Iq(x)  - Ij,(x)  - E 


(A-4) 


which  becomes  in  matrix  notation 


• «« 

I (x) 
n 

1 0 0 1 

• 

• 

0 10  0 1 

1 \ V , 

1 V.  V ‘ 1 

• 

SS 

I X \ ( t 

I^(x) 

1 's  \ 'v'  ' 

• 

• 

. ^ \ 0 1 
* \ % 

1 N 

Ij^(x) 

* Oil 

-1 -1  1 

"N — — ^ 

1 


LicW 


(A-5) 
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One  can  easily  show  (use  elementary  row  operations  to  reduce  T to  echelon 
or  upper  diagonal  form)  that  T is  nonsingular  and  therefore  represents  a 
unique  transformation.  Thus  for  a particular  x,  given  the  actual  line 
cun  ants,  1q(x),  one  can  obtain  the  components. 


• 

• 

Idi(x) 

• 

I^(x) 

• 

= T 

• 

• 

1di(x) 

• 

Ij^(x) 

.Ic(x)  . 

Again,  assuming  the  line  cross-sectional  dimensions  to  be  electrically  small, 

we  may  conclude  that  the  common  mode  currents  at  the  endpoints  of  the  line, 

I (0)  and  I (0»  are  essentially  zero  and  have  no  effect  on  the  terminal 
C C 

networks.  Therefore  it  suffices  to  consider  only  the  differential  mode 
(transmission  line  mode)  currents  when  computing  only  the  terminal  responses 
of  the  line. 

For  a parallel  discussion  of  this  problem  see  [14]. 
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APPENDIX  B 


WIRE 


Program  Listing 
Flowchart 
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C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


^,^^*^,t**»**  ******:****  ****************************  ********Mlt(ZOOO^ 

■IRE0002 

PROGRin  WIRE  RIRE0003 

(FORTRAN  IV,  DOUBLE  PRECISION)  NIRE0004 

WRITTEN  BT  WIREOOOS 

CLAYTON  R.  PAUL  WrBB0006 

DEPARTHENT  OP  ELECTRICAL  ENGINEERING  WIRE0007 

ONIVEHSITI  OF  KENTUCKY  BIRE0008 

LEXINGTON,  KENTUCKY  40506  WIRE0009 


C 

c 

c 

c**** 

c 

c 

c 

c 

c 


A DIGITAL  CONPUTER  PROGRAfl  TO  CONFUTE  THE  TERBINAL  CURRENTS 
AT  THE  ENDS  OF  A NQLTICONDUCTOR  TRANSNISSION  LINE  HHXCK  IS 
EXCITED  BY  AN  INCIDENT  ELECTRONAGNETIC  FIELD. 

THE  DISTRIBUTED  PARANETER,  NULTICOHDOCTOB  TRANSNISSION  LINE 
EQUATIONS  ARB  SOLVED  FOR  STEADY  STATE,  SINUSOIDAL  EXCITATION 
OF  THE  LINE. 

THE  LINE  CONSISTS  OF  N HIRES  (CYLINDRICAL  CONDUCTORS)  AND  A 
REFERENCE  CONDUCTOR.  THE  REFERENCE  CONDUCTOR  HAY  BE  A HIRE 
(TYPE=1),  AN  INFINITE  GROUND  PLANE  (TYPE=2) , OR  AN  OVERALL 
CYLINDRIdAL  SHIELD  (TYPE=3) . 

THE  INCIDENT  FIELD  HAY  BE  IN  THE  FORN  OF  A UNIFORH  PLANE  HAVE 
FOR  TYPE  1 AND  TYPE  2 STRUCTURES  OR  A NONUNIFORH  FIELD  FOR  ALL 
STRUCTURE  TYPES. 

THE  N HIRES  ARE  ASSUMED  TO  BE  PARALLEL  TO  BACH  OTHER  AND  THE 
REFERENCE  CONDUCTOR. 

THE  N HIRES  AND  THE  REFERENCE  CONDUCTOR  ARE  ASSUMED  TO  BE 
PERFECT  CONDUCTORS. 

THE  LINE  IS  IHHERSED  IN  A LINEAR,  ISOTROPIC,  AND  HOHOGENEOUS 
HEDIUH  WITH  A RELATIVE  PERNBABILITY  OF  NOR  AND  A RELATIVE 
DIELECTRIC  CONSTANT  OF  ER.  THE  HEDIUH  IS  ASSUMED  TO  BE  LOSSLESS. 

LOAD  STRUCTURE  OPTION  (ISO)  DEFINITIONS: 

LS0=11 ,THEVRNIK  EQUIVALENT  LOAD  STRUCTURES  WITH  DIAGONAL 
IMPEDANCE  HATHICES 

LSO= 12,THEVENIN  EQUIVALENT  LOAD  STRUCTURES  WITH  FULL 
INPEDANCE  MATRICES 

LS0=21, NORTON  EQUIVALENT  LOAD  STRUCTURES  WITH  DIAGONAL 
ADMITTANCE  MATRICES 

LSO=22, NORTON  EQUIVALENT  LOAD  STROCTOHES  WITH  FULL 
ADMITTANCE  MATRICES 

FIELD  SPECIFICATION  OPTION  (FSO)  DEFINITIONS; 

FS0=1, UNIFORM  PLANE  HAVE  (TYPB=1,2) 

FS0=2, NONUNIFORM  FIELD  (TYPE=1,2,3) 

FUNCTION  SUBPROGRAMS  USED:  E1,E2 
SUBROUTINES  USED:  LFQT1C 


HIRE0003 

HIBE0004 

HIRE0005 

HIRE0006 

HIRB0007 

HIREOOOS 

HIRE0009 

HIBE0010 

HIBE0011 

HIBE0012 

HIRE0013 

HTRE0014 

HIHE0015 

HIRE0016 

HIRE0017 

HIRE0018 

HIRE0019 

HIRE0020 

HIRB0021 

HIRB0022 

HIRE0023 

HIRB0024 

HIRB0025 

HIRB0026 

HIHE0027 

HIRE0028 

HIRE0029 

HIRE0030 

BIRE0031 

HTRE0032 

HIRE0033 

HIRE0034 

HTRB0035 

HIRB0036 

HrRE0037 

HIRE0038 

HIHE0039 

HIRE0040 

HIREOOHI 

HIRE0042 

HIRB0043 

HIRE0044 

HIHR0045 

HIRB004b 

HIRE0047 

HIRE0048 

HIRE0049 

HIREOOSO 

HIRB0051 

hire6o52 

HIRE0053 


HXRE00S4 

*•*•«*******««**«*•**•*«»***«******»****••»»«***•***«•**•***»*•**«*)(  jOgOOSS 

HIRB0056 

ALI  VECTORS  AND  MATRICES  IN  THE  FOLLOWING  DIMENSION  STATEMENTS  HIRE0057 

SHOULD  BE  OF  SIZE  N WHERE  N IS  THE  NUMBER  OF  WIR ES (EXCLUSIVE  OF  WIRE0058 
THE  REFERENCE  CONDUCTOR),  I. E. , V 1 (N) , V2 (N) , YO (N  ,N) , YL  (H , N) , B (N)  , WIBE0059 


A (N,N)  ,WA  (N)  ,M1  (N,N)  ,M?(N,  N)  , ETL  (N)  , ETO  ( N)  , H 3 (N,  N)  , V 3 (N)  , V4  (N) 


WIBE0060 

WIRE0061 
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mmciT  BBtL*8  (Jt-B,0-Z)  tISE0062 

IITEGKI  TTPB,rS0  IIBB0063 

BBBL*8  L.T3(  2) ,t0B,ILa,IOR,I0X,ItB,lLI,I0t,IL«,B0B,H0O2PI  VIBE006« 

1,nX,1II,HB,BP  BXBB0065 

caaPLBX«16  ZJ,ft(  2),T2(  2),r0(  2.  2i,ni  2,  2),B(  2,  2),B(  2},  BXBB0066 
1SORO,SOnL,XO,XL,ZBBOC«BB(  2)  .HM  2,  2)  ,R2(  2,  2)  ,XTL(  2),BT0(  2),  BXkBOOB? 

2C,i1,B2,OBBC,B3(  2,  2)  ,BBXL,BBL,f 1B«?2R,BP,Bf, BJBZ«BBTPBZ,BPBL,  BIBB0068 

3BHBL, BJCI,S0BC, SOBS, BLOC, BL0S,B1,B2,BJBT  0IBB00B9 

COBBOM  ZJ,tBBO,TBO,OKB,OBBC  BIBB0070 

DATA  CflTB/2.5*D-5/,H002PX/2.O-7/,P5/.5D0/,PO0H/*.D0/,  BIBB0071 

10HE80/t80. DO/, T/2. 99792508/  BXBB0072 

ZBBO^O.DO  BXBB0073 

ORB<=1.00  IXBB007t| 

T00>2. DO  BXRE0075 

ORBC«OC8PLX(1.DO.O.DOl  BXRB0076 

ZBBOC*DCRPlZ (0.00,0.  DO)  VIBE0077 

IJ^DCBPLX (0.00,1. DO)  IIBB0078 

PX>POOB*OBTtR(OBB)  BIBB0079 

8ADBG>0RB80/PX  VXBBOOBO 

VXBZ0081 

•♦PSEQOBlICt  IBDEPEBDBBT  CALC0L*TIOiS**************M***-*************BIRB0082 

BXRB0083 

BEAD  AID  PBXBT  XRPOT  DATA  BIBBOOBH 

■XBB008S 

BBA0(5,1)  TTPB, ISO, P"0,B,ER,B0R,L  VIRB0086 

1 P0RRAT(9Z,X1,8X.X2,9X,X1,8X,X2,3(B10.3))  BXRE00B7 

IF(TTPB.GB.1.AB0.TTPB.LE.3|  GO  TO  3 BIREOOBB 

BRXTB(6,2)  BXBE0089 

2 POBHATC  STROCTORB  TfPB  BBROB'//*  TTPE  ROST  BQOAL  t,2,0R  3*///)  BXBB0090 

GO  TO  121  BIBB0091 

3 XP(LSO.BQ. 11. OB. LSO.BQ. 12)  GO  TO  S BIRB0092 

XP(lSO.BQ. 21.0R.LS0. BQ.22)  GO  TO  5 BIRB0093 

■RITB(B.4)  IXRB009« 


BBA0(5,1)  TrPB,LSO,P"0,B,ER,B0R,L  VIRB0086 

1 P0RRAT(9Z,X1,8X.X2,9X,X1,8X,X2,3(B10.3))  BXRB00B7 

IF(TTPB.GB.1.AB0.TTPB.LB.3)  GO  TO  3 BIRBOOBB 

BRXTB(6,2)  BXBE0089 

2 POBHATC  STROCTORB  TfPB  BBROB'//*  TTPE  ROST  BQOAL  1,2,0R  3'///)  IXBB0090 

GO  TO  121  BIBB0091 

3 XF (LSO.BQ. 11. OR. LSO.BQ. 12)  GO  TO  S BIRB0092 

XF (LSO.BQ. 21. OR. LSO. BQ.22)  GO  TO  5 BIRB0093 

■RITB(B,4)  BXRB009« 

A FORBATC  LOAD  STROCTORB  OPTXOB  BBROB'//'  LSO  HOST  BQOAL  1 1, 12, 28XRB009S 
11, OB  22'///)  BIBB0096 

GO  TO  121  RIRB0097 

5 IF<FSO.BQ. 1,0a.PSQ.BQ.2)  GO  TO  7 HXRB0098 

BRITB(6,6)  HIBB0099 

6 POBHATC  FIBLD  SPBCIFICATIOR  OPTXOB  BBBOR'//'  FSO  HOST  BQOAL  1,VIRB0100 

1 OB  2'///)  flBBOIOI 

GO  TO  121  HIRB0102 

7 IP(TTPE.BQ.3.AHD.FS0.BQ. 1)  60  TO  8 BIBB0103 

GO  TO  10  BIRB0104 

8 HRITE(6,9)  HXRB0105 

9 FORBATC  OHXrOBH  PLABB  lATB  EXCXTATIOB  CABBOT  BE  SPBCXFXBD  FOR  THBBXRB0106 


1 TXPE  J STROCTORB'///) 

GO  TO  121 

10  tBITB(6,11)  B, TIRE, LSO, FSO, L,ER,HOR 

11  P0HHAT(1h1, SIX, 'SIRE'/// 

1«5X,X2,'  PARALLEL  HIRES'/// 

2A3X,'  TTPE  OP  STROCTORB*  ',11/// 

341X,'  LOAD  STROCTORB  OPTION*  ',12/// 

«*0X,'  FIELD  SPECIMCATIOM  OPTIOB*  ',11/// 

539X,'  LIRE  LENGTH*  •,1PE13,6,'  HBTBRS'/// 

633X,'  DIELECTRIC  CONSTANT  OF  THE  HBDIOB*  ',1PE10.3/// 
733X,'HBlATHrE  PERHEABILITI  OF  THB  HEDIOH*  • , 1PE1 0. 3///) 
GO  TO  (12,20,16) ,TrPB 

12  READ(S,13)  HHO 

13  FORBAT(5X,B10.3) 

RRITB(6,1<t)  RHO 


BXRB0107 

HIBE0108 

NIRE0109 

HIRE0110 

HIRB0111 

HIRE0112 

HIBB0113 

VIBE0114 

0IRB0115 

HXRE0116 

HIBB0117 

HIBE0118 

HXRB0119 

HIRB0120 

HIRB0121 


14  FORHATC  REFERENCE  CONDUCTOR  FOR  LINE  FOLIAGES  IS  A BIRE  WITH  HADIHIBB0122 
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10S=  aiLSV///)  III*B0123 

HB0=fiW0*CaiH  tIBB012* 

BRITB(6,15)  WXRB012S 

15  POBB&T(*  VISE  RaHBBR',4X»’BIBB  RRDIOS  (HILS)',18X,  fIIB0126 

1'Z  COOBDIHATE  (HETBRS) * ,241, 'T  COORDXMtTE  (HETERS| • ,//)  RIRE0127 

GO  TO  23  RIRB0128 

16  RERD(5,17)  RS  RIRB0129 

17  FORnkT(SX,B10.3)  RIRB0130 

fRITE(6,18)  RS  tIHE0131 

18  FORHATC*  REFERERCE  CONDOCTOR  FOR  LIRE  VOLTAGES  IS  A CTLIRDRZCAL  OTRXRE0132 

1ERALL  SHIELD  HITH  IRTEBIOR  RADIOS^  ',18810.3,'  HBTBRS'////)  HXRE0133 

RS2>RS*RS  IXRB0134 

8RXTX(6.19)  HXRB0135 

19  FORHATC  HIRE  RDHBEB' ,2X,'8XRB  RADIOS  (HXLS)',  2X, • SBPARATIOR  BBTHIXRB0136 

IBER  HIRE  ARO  CEFTER  OF  SHXBLB  <RETBtS) ',61, • ARGOLAR  COOROIRATB  (OBRXRB0137 
2GREBS) *//)  RXIB0138 

GO  TO  23  VXaB0139 

20  RRXTB(6,21)  RIRBOUO 

21  FORHATC  REFERENCE  CONDOCTOR  FOR  LXHB  VOLTAGES  IS  AR  IRFIRXTB  GRO0RI1B0141 

1RD  PLANE'////)  RI1B0142 

RRITB(6,22)  RIBB0143 

22  FORHATC  HIRE  NOnBER • , «X,' HIRE  RADIOS  (HUS) ',181,  HIBB01M 

I'HORTZORTAL  COOROIRATB  (HETBRS) • ,161, 'HIRE  HEIGHT  (HBTERS) '//)  HIRE01A5 

HIRBQ146 

READ  ARD  PRIRT  LINE  OIHBRSIORS  ARD  COHPDTB  THB  CHARACTERISTIC  HIRB0147 

IHPEDARCE  HATRIX,ZC  (STORE  ZC  IN  REAL  PART  OF  ARRAY  HI)  HIRBOUS 

c HIRB0149 

23  C>H002PI«QHEC*V*0SQRT(H0R/BR)  HIRB0150 

DO  29  1^1, R HIRB0151 

RBA0(5,24)  RH,Z,T  HZRB0152 

24  FORHAT(3  (SX,B10.3) ) BIIB0153 

HRITB(6,25)  I,Rf,Z,Y  NIRB0154 

25  FORHAT(2X,I2,13X,1PB10.3,2TX,1PR10.3,35X,1PB10.3/)  •IRB0155 

T3(I)-t  HZIB0156 

V4(I)-t  HItB0157 

RH«RR*CHTH  VIRB01S8 

GO  TO  (26,27,28)  .TYPE  RZRB0159 

26  DI2«Z«Z«Y*Y  flRBOHO 

HI (I,I)»C40LOG(0I2/(RH«BH0))  RIRB0161 

GO  TO  29  BIEB0162 

27  HI  (I,I)=*C*OLOG(THO*Y/RH)  HIRBOieS 

GO  TO  29  HIRB0164 

28  HI (I,I)>C*OLOG((RS2-Z*Z)/(RS*RB))  RIRB0165 

29  CORTIROE  RIBE0166 

IF(H.BQ. 1)  GO  TO  34  HIRB0167 

K1>R-1  HIRB0168 

DO  33  I«1,K1  HZBB0169 

K?>I»1  BIRB0170 

DO  33  J>K2,R  HIRB0171 

ZI«V3(i;  fIRB0172 

7iJ>V3(J)  Ri'RB0173 

YI>T«(I)  HIRB4174 

YJ>T4(J)  NIRB0175 

GO  TO  (30, 31,32) , TYPE  HIRE0176 

30  DI2’=ZI4ZI*II*YI  HIRE0177 

DJ2>ZJ*ZJ'IJ«YJ  HIREOIYS 

ZD=ZI-ZJ  HIHE0179 

YO=YI-YJ  HIRE0180 

DXJ2^Z0«ZD«Y0*fD  HIBE0181 

Hi  (I,J)=PO''C’»DLOG(DI2*OJ2/(BiO*RaO*BIJ2) ) HIRB0182 

HI  (J,I)»H1 (I,J)  HIRR0183 
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c 

c 

c 

c 


c 

c 

c 

c 

c 


GO  TO  33 

31  ZD=ZI-ZJ 
yD=YI-YJ 

DIJ2=ZD»ZD*YD»YD 

HI (I,J)=P5*C*DLOG(ONE+POUR*YI*YJ/DrJ2) 

HI  (J,I)  = H1  (I.J) 

GO  TO  33 

32  THET*=(YI-YJ)/RADEG 
RI2=ZI*ZI 
RJ2=Za*ZJ 

HI  (I,  J)=P5*C*DL0G(  (RJ2/RS2)  ♦ (Rr2*RJ2*RS2*RS2-THO*ZI*Z.T*HS2* 
IDCOS(THETA) ) / (RI 2* RJ 2»RJ2*aj 2-T«0«ZI*Z J*RJ2* DCOS  (THETA) ) ) 

HI  (J,I)=H1  (I,J) 

33  CONTIHGE 

COHPOTE  THE  INVERSE  OF  THE  CHARACTERISTIC  IflPEDANCE  HATRIX,ZCINV 
(STORE  ZCINT  IN  ARRAYS  n2  AND  N3) 

34  DO  36  1=1, N 
DO  35  J=1,N 
A(I,J)=H1{I,J) 

3-:  H2  (I,J)=ZEROC 

36  12 (I,I)=ONEC 

CALL  LBQT1C{A,N,N,H2,N,N,0,1(A,KER) 

KEB=RER-128 
DO  37  1=1,  N 
DO  37  J=1,N 

37  H3  (I.J)  =H2  (I,J) 

IF (KER.NE.  1)  GO  TO  39 
HRITE(6, 38) 

.3  FORHATC//, ' ♦♦♦♦♦CHARACTERISTIC  IHPEDANCB  HATRIX  INVERSION  ERROR 
GO  TO  121 

READ  AND  PRINT  ENTRIES  IN  LOAD  A DHITTAKCE (IH PEDANCB)  HATRICES 
(STORE  ADHITTANCB(IHPEDANCB)  HATRICES  AT  X=0  IN  ARRAY  YO  AND 
THOSE  AT  r=L  IN  ARRAY  YL) 

39  IF  (LSO.EQ,  II.OR,  Li.0.  EO.  12)  GO  TO  42 
»PITE(6,40) 

40  PORHAT(//, 18X,'ADHITTANCE  AT  X=0' ,43X, • ADHITTANCE  AT  X=L'/) 
WRITE(6,41) 

41  FORHAT(21I,' (SIEHENS) ' ,51X,*  (SIEHENS) */) 

GO  TO  45 

42  HRITE(6,43) 

43  FORHAT(//, lex,' IHPEDANCE  AT  i=0* ,44X,* IHPEDAHCE  AT  X=L'/) 
*RITE(6,44) 

4 4 FORHAT  (2  3X,'  (OHHS)  •,  341,  • (OHfiS)  • /) 

45  R8ITE(6,46) 

46  FORHAT  (•  ENTRY  • , 10X,  • REAL*  , HX  . ' IN  AG  • , 4 1 X,  ‘ H EAL' , H X,  ' IH  AG  • //) 

DO  49  1=1, N 

READ(5,47)  TOR, YOI , Y LR , YLI 

47  FOHHAT(2(E10,3)  ,20X,2(E10.  3)) 

Y0(I,I)=Y0R»XJ4Y0I 

YL  (I,I)  = YLB*XJ^YLI 
HRITE(6,48|  I,I,I0  (1,1)  ,YL(I,I) 

48  FORHAT(1I,I2,2X,I2,2 (5X, 1PE10. 3)  ,30X,2 (SX, 1PB10.  3)  /) 

49  CONTINUE 

IP  (LSO.EQ. II.OR. LSO. 80.21)  CO  TO  52 
IF(N.BQ. 1)  GO  TO  52 
DO  51  1=1,  K1 


HIREOIBN 
HIREOiaS 
W1RE0186 
WIRE0187 
HIRE0188 
HIRE0189 
HIRF,0190 
HIRE0191 
WIRE0192 
HIRE0193 
HIRE0194 
HIRE0195 
WIREOISft 
HIRE0197 
HIRB0198 
HIRE0199 
HIRE0200 
HIRB0201 
RIRE0202 
8IRP.0203 
HIRE0204 
8IRE020 
NIRR0206 
HrRE0207 
»inR020H 
WIRE0209 
WrRE02  10 
WrRE021 1 
WTRP02  12 
WIRE0213 
♦♦HIRR02  14 
WIRE0215 
WIRE0216 
HIRE0217 
HIRE0218 
HIBB0219 
WIRE0220 
HIRE0221 
WIRE0222 
WIRE0223 
HIRE0224 
WIRE0225 
HIBE0226 
VIRE0227 
HIRR0228 
WIRE0229 
WIRE0230 
WIRE0231 
HIRE0232 
UIRE0233 
HIRE0234 
HIRRQ235 
WIRE0236 
«IRE0237 
WIRE0236 
HIRR0239 
HIRE0240 
WIRE0241 
HrREn?42 
WI  RR024  1 
WrRR0244 
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c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


K2=I+1 

DO  51  J=K2,N 

BBAD(5,50)  y0P,y01,YL«,JtI 

50  FOBNAT(2(B10.3»,20X,2(E10.3)) 

yo  (i,jj  =roBtxj*yoi 
TL  (i,j)  = yLfi*xj*rLi 
yo  (j,i)=YO (i,j) 

yL(j,ij=yL  (i,j) 

HRITE(6,48)  I,  J,  TO  (I , J)  , YL  (I , JJ 

51  CONTINUE 

IF  THEVENIN  EQUIVALENT  SPECIFIC,  SNAP  ENTRIES  IN  ARRAYS  Hi  AND 
HI  WILL  CONTAIN  THE  INVERSE  OF^ZC  AND  H2  BILL  CONTAIN  ZC 

52  IF(LS0.EQ.21.0R.LS0.EQ.22)  GO  TO  54 
DO  53  1=1, N 

DO  53  J=I,H 
A1=H1  (I,  J) 

A2=M2 (I, J) 

HI  (I,J)*A2 
HI  (J,I)=A2 
H2(I.Jl=A1 

53  H2  (J,I)=A1 

COHPOTE  THE  HATBIX  2C»2L*ZCIHV*Z0  FOB  THE  THBVENIN  EQUIVALENT 
OR  ZCIHV»YL*2C*Y0  FOB  NORTON  EQUIVALENT 
STORE  IN  ARRAY  H2 

54  IF (LSO.EQ. 12.0R.LS0. EQ.22)  GO  TO  57 
DO  55  1=1, N 

DO  55  J=1,H 

55  A(I,J)=*H1(1,J|*I0(J,J) 

DO  56  1=1, N 

DO  56  J=1,N 

56  H2  (1,J)  =YL  (1,1)  •Ad,  J)  ♦fl2<I,  J) 

GO  TO  62 

57  DO  54  1=1, H 
DO  59  J=1,N 
SOnL=ZBROC 
DO  58  K=1,N 

58  S0HL=S0HL«H1 (I,K)*T0  (A,J) 

54  A(I,J)=SUHL 

DO  61  1=1, N 
DO  61  J=1,N 
SUnL=ZBROC 
DO  60  K=1,N 

60  SUHL=SUHL*YL  (I,K) •A (K, J) 

61  H2  (I,J| =SUHL*H2  (I.J) 

62  B8=THO*PI*DSQ8T(ER*HDB)/V 
B8L=BB*L 

IF  FIELD  SPECIFICATION  IS  A ONIFOHH  PLANE  WAVE,  BEAD  DATA  AND 
COHPUTE  THE  COHPONENTS  OF  THE  ELECTRIC  FIELD  INTENSITI  AND  THE 
PHASE  CONSTANT(FOB  ONE  HERTZ)  IN  THE  X,y,  AMD  2 DIRECTIONS 

IP(FSO.EQ.2)  GO  TO  66 
READ(5,63)  EH,THE,THP, PHP 

63  FORHAT(4(B10.3,5X)) 

HRITE(6,64) 

64  FORMAT (///'  EXCITATION  SOURCE  IS  A UNIFORM  PLANE  WAVE'//) 
HRrTE(6,65)  BH,THE,THP, PHP 


HIBE024S 

WIRE0246 

WIRE0247 

NIRE0248 

BXRE0249 

HIRE0250 

WIRE0251 

HIRE0252 

HIBE0253 

WIBE0254 

HIRE0255 

H2HIRE0256 

WIHE0257 

HtRE0258 

WIRE0259 

HIRE0260 

HIRE0261 

HIRE0262 

HIRE0263 

HIRE0264 

HIRE0265 

HIRE0266 

VIRB0267 

HIRB0268 

BIRE0269 

HIRE0270 

HIRE0271 

HIRE0272 

HIRE0273 

WIRE0274 

HIRE0275 

WIRB0276 

WIRE0277 

HIRE0278 

WiaB0279 

WIRB0280 

WIRE0281 

WIRE02B2 

HIRE0283 

WIRE02a4 

HIRE0285 

HIRE0286 

HIRE0287 

WIBE0288 

HIRE0289 

HIRE0290 

HIRB0291 

WIRE0292 

HIRE0293 

HIRE0294 

WiaB0295 

HIRE0296 

WIREO‘297 

NIRE0298 

NI8F.0299 

HIRB0300 

WIRE0301 

HIfiE0302 

NIRB0303 

WIRE0304 

WIBE0305 
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I 


5 

1 

65  rOEfl»T(*  RAGIITOOE  OF  ELECTBIC  FIELD  = *, 

1PE10.3,' (FOLTS/HETER) •/  WIRE0306 

1- 

1*  THETBB  * '.IPEIO.!,' (OEGBBES)  V*  TMEMP 

= ^IPEIO.S,*  (DEGREES)  VHrRE0307 

2*  PHIP  * »PB10. 3, • <DK5FEFS) •////) 

IIIRE0308 

t: 

THB=>THE/BBOEG 

HrRE0309 

TMP>THP/BBOBG 

HIRE0310 

PHP»PRP/BBDEG 

HIRE0311 

CTE«DCOS (IHEJ 

HIRE0312 

f. 

CTP>OCOS  (THP) 

WIRE0313 

k 

CPP«DCOS (PBP) 

H1RE0314 

STE-DSII (THE) 

HIRE0315 

1 

STP^OSIR (THP) 

0IREO316 

S 

SPP>DSIB (PHP) 

W1RK0317 

r 

BIH«BH*CTE*STP 

HIBE0318 

P 

EZR>-Ea* (CTB*CTP*CPP-STB*SPP) 

WIRB0319 

B1H»-BH* (CTE*CTP*SPP*STE*CPP) 

MIRE0320 

BBl*BB«StP*SPP 

HIRB0321 

1 

BBT^BB^CTP 

VIRE0322 

p 

BB1*BB*STP*CPP 

IIIRE0323 

jS 

c 

HIRE0324 

f 

C»****FIEQOE«CI  OEFBIOEET  C4LC0I.AT1OES********* 

C 

WIRE0326 

1 

<«  COHTINOE 

0IREO327 

f 

RBBD(5.67,BBOa121)  F 

0IREO32B 

i 

67  FORHftT(Z10.3) 

HIRE0329 

'; 

BBTi*BB«F 

0IREO33O 

f 

BETAL-BBL«F 

HIRE0331 

s 

OSaDSIBCBBTftL) 

HTRE0332 

t 

V 

DC«OCOS(BBTBL) 

HIRB0333 

i. 

GO  TO  (68,74), FSO 

HIRE0334 

c 

0IREO335 

h 

c 

COHPOTB  THE  EQUIftLBNT  FORCING  FUNCTIONS 

FOR  UNIFORB  PLANE  MATE  0IREO336 

\ 

c 

BZCmTION 

0I8EO337 

\ 

c 

0IREO13B 

c 

COHPUTB  THE  1,1,  AND  2 CONPONENTS  OF  THE 

PHASE  CONSTANT  FOR  HIRB0339 

i 

c 

OHIFCRN  PLANE  HATE  EICITATION  AND  A FREQUENCY  OF  F HERTZ  VIRE03U0 

1 

c 

0IREO341 

r 

68  BI«BBX*F 

0IREO342 

s 

BT«BBT*F 

0IRRO343 

1 

BX«BBZ*F 

0rREO344 

1 

EBZL-CDBZP  (-U*BZ*L) 

0IREO345 

i 

BP«8ETAtBl 

0IREO346 

r 

BNsBITA-Bt 

0IRRO347 

1 

BPBL>COEIP(ZJ4BETAL) 

HIRB0348 

1 

BNBL-COBZP (-IJ4BETAL) 

W’IRR0349 

f 

BP>BPBL*I2 (ZERO, L, -BP) 

0IREO35O 

1 

ENsBRBL«E2 (ZBRO,L,BR) 

0IREO351 

1 

69  CO  TO  (70,72), TYPE 

0IREO352 

I' 

c 

0IRRO353 

1 

c 

CORPOTB  FORCING  FUNCTIONS  FOR  UNIFORH  PLANE  HAVE  EXCITATION  AND  HIRE0354 

1 

c 

TYPE  1 STROCTURES 

0IRBO38S 

i- 

c 

0IREO356 

1 

70  DO  71  I«1,N 

0IRBO857 

TI»?4(I) 

0IREO358 

r 

ZI*73(I) 

0IREO359 

1 

BrP8Z«BT4II»3E*ZI 

HIRE0360 

BBTPBZ<COBZP(-|J«BYPBZ) -ONE 

0IPEO361 

\ 

EJBZ«CDBXP (-XJ«BZ*ZI) 

HIRE03e2 

EJBIoCDEXP (-ZJ4BT4TI) 

HIRE0363 

? 

7tH>BZR*EBIPBZ/THO 

0IREO364 

t 

V2H»-XJ*71B 

0IBEO365 

1 

11  (I)«T1R*(RP»EN) 

HIRR0366 
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I 

t 

> 


i 

I 


i 

■f 


V 

1 

1 


|: 

'■i. 


I 

>■ 


d 


onoo  orjo  onnr. 


V2(I)=V2H*(EP-B»)  W1HE0367 

ETO(Z)  = (ETB*TI«'EZR*ZI)«B2<ZEBO,OH,-BrPBZ)  VIRB0368 

71  ETL  (I) =BT0 (I) *EBXt  WIBB0369 

GO  TO  96  VIBE0370 

BIBB0371 

COHPOTE  POBCIMG  PDNCTIOHS  FOB  OBIPOBtl  PLAHE  BATE  BICITATIOH  ABO  BIRE0372 
TYPE  2 STBOCTOBES  B1RE0373 

BIRE037« 

72  DO  73  I=1,I(  BIBE0375 

TI=y4(I)  BIRB0376 

Zr=V3(I)  BIRE0377 

SBY=DSIH  (BY*n)  BIHB0378 

BJBZ=CDEXP{-XJ*BZ*ZI)  BIRB0379 

V2«--'SXR*EJBZ*SBY  BIRE0380 

Vin=XJ*T2n  BIRE0381 

VI  (I)»T1fl*(EP»B»)  BIRB0382 

V2(I)»T2H*(EP-BH)  BIBE0383 

ETO  (1) *ETB*EJBZ*B2 (-tl, YI, BY)  BIRB038tt 

73  BTL(I)»BTO(I)*EBXI.  BIRE038S 

GO  TO  96  BIRB0386 

BXRB0387 

CORPOTB  THE  B00I7ALBRT  POBCIMG  PBRCTIOBS  POB  HONUNIPOBK  BXCITATIOMBIBE0388 

BIRB0389 

74  BPBL^CDBZP(IJ*BETAL)  BIBB0390 

BMBL=COBXP(-XJ«BETAL)  BXRB0391 

B8XTE(6,75)  BXRE0392 

75  PORBAT(///'  EXCXTATIOM  Sl/OBCE  XS  A MOMOMXPOBB  PIBLD'//>  BIRB0393 

GO  TO  (76,83,83),  TYPE  BXBB0394 


COnPOTB  THE  CONTRXBOTXOH  DOE  TO  THE  LORGXTOniMAL  ELBCTBXC  PIBLD 
POB  TBB  BEPEBEMCB  BXRE 

76  BEA0(5,77)  ML0,E0,T0 

77  POFKAT(X10,2(10X,E10.3>) 

MRITB(6,78) 

78  PORBAT(*  tOMGXTUDXHAL  ELECTRIC  PIBLD  OM  REPERBHCB  BIBB*/) 
BBITB(6,79) 

79  PORnAT(5X,*SPECXPXCATIOH  POXBT (METERS) •, 

15X,' ELECTRIC  PIELO  IMTEHSITY (VOLTS/MBTER) •, 5X, • PHASE (DEGREES) V/) 
IL^ZERO 
EL>E0 
TL=T0 

SOHC-ZBROC 

SOHS^ZEROC 

BRITE(6,80)  XL,E0,T0 

80  FORHAT(13X,1PE10.3,24X,1 PE  10.3,251,1  PE  10. 3) 

DO  82  I<1,IL0 

READ(5,81)  ZI,E1,TI 

81  PORNAT(3(E10.3,10X)) 

BRITB(6,80)  ZI,8I,TI 
ZP^XI 

EP  = EI 

TP=TI 

X0*XP-XL 

Hla(EP-EL)  /XD 

BI^(EL*XP-BP*XL)/XD 

MI= (TP-Tl) /(RADEG*XD) 

CI= (TL*XP-TP*XL) /(BADEG^XO) 

«M»NI-BETA 

MP«NI*BETA 

EJCI=CDBXP(XJ4CI) 


BIBB0390 

BXRB0391 

BIRE0392 

BIRB0393 

BIRB0394 

BIRE0395 

BIRB0396 

BIRB0397 

BIRB0398 

BIRB0399 

BIRB0400 

BISB0401 

BIRB0402 

HIRB0403 

BIBE0404 

BIRE0405 

BIRB0406 

BIRB0407 

BIRE0408 

BIBB0409 

BIBB0410 

HIRB04lt 

BTBE0412 

HIRE0413 

BIRB0414 

BIRB0415 

BIBB0416 

BIRE0417 

BIBE0(t18 

BIRE0419 

HIRC0420 

BIRE0421 

BIBB0422 

BIRE0423 

HIRB0424 

BIBE0425 

BIRB0426 

BIRR0427 
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S0HC=S0I!C*BJCI*P5*<BI*KPBL*E1  {XL,IP,IIH)  ♦HI*EIIBt*E1  (XL, IP, HP) 
UBI*BPBL*E2(XL,XP,Rn)»BI*ENBL*B2(XL,XP,HP)) 

S0aS=SaaS-XJ*P5*EJCI*(8I*EPBL*EI(Xt,XP,IIH)-HI*ENBL*E1(XL,XP,NP) 

1 «BI«EPBL*B2 (XL, XP, HH) -BI *ENBL«E2 (XL. XP.HP) ) 

XL=XP 

EL=EP 

82  TL=TP 
ELOC=SOHC 
ELOS^SaHS 
GO  TO  8« 

83  ELQC^ZBBOC 
ELOS=ZBBOC 

C 

C COBPDTE  THE  CONTRIBUTION  DOB  TO  THE  LONGITUDINAL  ELECTRIC  FIELD 
C FOB  THE  HIRES 
C 

8N  DO  95  1=1, N 

R£AD(5,8S)  HL0,EL0,TL0 

85  FOBnAT(I10, 2(101, BIO. 31) 

XL=ZERO 

HPITE(6,86)  I 

86  FORHAT(//‘  LONGITUDINAL  ELECTRIC  FIELD  ON  HIRE  •,3X,I2/) 
HBITE(6,79) 

HRITE(6,80)  XL,BL0,TL0 

BL=BL0 

TL»?L0 

SUBOZBROC 

SDBS=ZBROC 

DO  38  J=1,RL0 

READ(5,87)  XI.EI.TI 

87  FOBaAT(3(E10.3,10X)) 

HBITE(6,80)  XI.EI.TI 
XP=XI 

EP*EI 

TP»TI 

X0»XP-XL 

HI>(EP>EL) /XD 

BI*(EL*XP-EP*XL)/XD 

HI=  (TP-TL)/(RADEG*XD) 

CI= (TL*XP-TP*XL|  /(BADEG*XD) 

HH*NI-BETA 

HP=NI*BETA 

EJCI»CDEXP(XJ*CI) 

SUHC«SUHC*EJCI*P5* (HI*EPBL*E 1 (XL, X P, NB)  *HI*ENBL*E1 (XL, XP.HP) 
l4BX*BPBL*E2(ZL,XP,Hn) «Bt*ENBL*E2 (XL, XP,KP) ) 
SaHS»S0HS-XJ*P5*EJCI*(HI*EPBL*E1 (XL.XP.NH)  -HI* BNBL*E 1 (XL, XP, NP) 
1*BI*EPBL*E2(XL,XP,Nn)-BI*ENBL*E2 (XL.XP,HP)  ) 

XL=XP 

EL*EP 

88  TL=TP 

VI (I)«SDHC-EL0C 
V2(I)=S0HS-EL0S 
C 

C CORPUTB  THE  CONTRIBUTION  DUB  TO  THE  TRANSVERSE  ELECTRIC  FIELD 
C AT  X*0  FOB  EACH  HIRE 
C 

XL=ZBRO 
HRITB(6,89)  I 

89  FOBHf,T(//'  TRANSVERSE  ELECTRIC  FIELD  AT  X=0  FOR  HIRE  •,3X,I2/) 
BEAD(5,85)  NT0,EET0,TT0 

HRITB(6,79) 


HIRE0428 

HIRE0<J29 

HIRE0430 

HIRR0431 

HIRE0432 

HIRE0433 

HIRE0434 

HIRE0435 

HIRE0436 

WIRE0437 

HIRE04  38 

HTRE0433 

HIRE0440 

VIRE0441 

WIRE0442 

HIHE0443 

WIRB0444 

KIRE0445 

HIRE0446 

HIPEOUH? 

HIRB0448 

HIRE0449 

HIRE0450 

HIRB0451 

HIRB0452 

HIREOilbl 

VIRE0454 

HIRn0455 

HIRE0456 

HIRE0457 

HIRE0458 

HIRE045‘} 

WIRB0460 

VIRR0461 

HIRE0462 

HIRB0463 

HIRB0464 

HIRE0465 

HIRE0466 

HIRE0467 

HTRB0468 

HIRE0469 

HIRE0470 

HiaB0471 

HIBB0472 

HIRE047J 

HTRE0474 

HIRE0475 

HrRE0476 

HIRE0477 

VIRE0478 

H1RE0479 

HIRE0430 

HIRB048  il 

WIRE0482 

HIRE0483 

HIRE0484 

HIRE0485 

H1RE0486 

WIRE0487 

HIRE0438 
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n n n n n 


HRITE(6,80}  Xl,EETO,TTO 

EL=EBTO 

TL*TTO 

SDRO^ZESOC 

DO  91  J=1,HTP 

REkD(S,90)  XI,E1,TI 

90  FOBHXT(3(E10.3,10X)) 

HRITZ(6,80)  XI,rl,TI 
XP=XI 

BP=EI 

TP*TI 

XD=XP-XL 

M»  (BP-EL)  /ID 

BI»=(BL*XP-BP«L)/XD 

III>(TP-TL)/(>BDE6*XD) 

Cl* (TL*XP-tP*XL) /(RADEG*  XO) 

EJCI=CDBXP(XJ*CI) 

SURO*SUSO»BJCI« (HI *81 (XL«XP, HI) (ZL,XP,HI) ) 

XL*XP 
EL*  BP 

91  TL«TP 
ETO(i) >sono 

COHPDTB  THE  COHTBIBOTIOM  DOE  TO  THE  TRAMSTBRSB  ELECTRIC  FIELD 
AT  X-L  FOB  EACH  HIRE 

XL>ZERO 
HRITE(6,92)  I 

92  FOBnAT(//'  TRAHSTBRSE  ELECTRIC  FIELD  AT  X*L  FOB  BIBB  •,3X,I2/) 
RBA0(S,8S)  MTL,EETL«TTL 

HBIT2(6,79) 

HRITB(6,80)  XL,BBTL,TTL 

EL*BBTL 

TL*TTL 

SU  1L«ZBB0C 

DO  99  J«1,HTL 

BBAD(5«93)  XI,BI,TI 

93  FORHAT(3(E10.3«10X)) 

RRITE(6,80)  XI,BI,TI 

xp«x: 

EP*ET 

TP*TI 

XD*XP-XL 

BI*(EP-EL) /XD 

ni»(EL*XP-BP*XL) /XD 

HI=(TP-TL) /(RADBG*XD) 

CI= (TL*XP-TP* . ) /(RADEG* XD) 

BJCr.*COBXP(XJ*CI) 

SUHL*S08L* EJCI*(nr*El (XL,XP,NI)*BI*E2(XL,XP,HI)) 

XL*XP 

EL*BP 

94  TL*TP 
ETL(I)*SaNL 

95  COHTIBOE 

CORPOTB  THE  TERHIMAL  COBBEHTS 
FORH  THE  EQUATIONS 

96  IF  (LSO.vTQ.  12.0R.LS0.  EQ.22)  GO  TO  100 
DO  98  I*1,H 


H1RE0989 
HIBB0990 
HIRB0491 
HIBB0992 
HIRB0993 
HIRB0994 
HIBB0995 
HIIB0496 
fIRB0497 
HtBE0498 
HIRE0499 
HIBE0500 
HIRE0501 
HIRB0502 
HIBB0503 
HIRE0S04 
HIRB0505 
HIRE0506 
HIRB0S07 
HIBB0508 
HIRB0509 
BIRE0510 
HIBB0S11 
HIBE0512 
B1RE0S13 
HIPB0S14 
HI3B0S15 
HI 180516 
HItB0517 
HIBE0518 
HIRB0519 
HIREOS20 
HIBE0S21 
HIBE0522 
HIRE0523 
HIRB0524 
HIRB0525 
HIRB0526 
HIBB0527 
HIBE0528 
HIRB0529 
■IRE0530 
HIRB0S31 
HIRB0532 
HIRE0S33 
HIRB0534 
HIRE053S 
HIRB0536 
HIRB0537 
HIRB0338 
HIBE0539 
HIRB0540 
HIRE0541 
HIRE0542 
HIRB0S43 
HIRE0544 
HIRE054S 
HIRB0S46 
BIRB0547 
HISE0S18 

•'IRB0S49 


L6 


o o o n n n 


SDnO=ZB80C 
SDnL=ZBBOC 
DO  97 

A(r,J)=XJ*D5*H2{I,J) 

SOBO=SOflO+H3 (I,J) •V2 (J) 

97  SOBl,*S0nL+B3(I,J)^BT0(J) 

vi(r)=yi  (i)-Bn(r) 

HI  <1,I)*SOHO 
ETi  (I)*S0Hi 

98  A (1,1)  ^ A (1,1)  tDC*  (rO  (1,1)  *11  (I -I)) 

DO  99  1=1, N ' 

SOHO»OHEC 

SOHt'OBBC 

IF(tSO,EQ.21)  snH0=TL(I,I) 

IF (ISO. EQ.  11)  S0HL=TI.(I,I) 

’’  GO^TO^IO?*’’  ♦*J*SOHL*H1  (1 , 1)  ♦ OC*SOH0* ETO (I)  ♦XJ»DS»SOHL*BTI.(I) 

100  DO  102  1=1, K 
SOHO=ZEROC 
SOHI.=  ZEBOC 
DO  101  J=1,M 

A(I,J) =XJ*DS*H2 (I,J) iOC* (10 (I,J)  ♦Tt(I,J) ) 

SOBO=SOHO*H3(I,J)*?2(J)  ' 

101  S(JBL=SaBL«^B3  (I,J)  *BT0(J) 

T1  (I)*f1  (I)  *871(1) 

HI  (1,1) =SUH0 

102  BTS.(I)=SOHt 
DO  106  1=1, B 

IF  (ISO.  EQ.  22)  GO  TO  104 
SOHO=ZEBOC 
SOHl=2EHOC 
DO  103  J=1,M 

S0H0=SUH0*Tt(I,a)*Hi (J,J) 

103  S0HL=S0Ht*TL (I,J) •BTL(J) 

B(I)’>V1(I)  »XJ*SUH0*DC*BT0(I)  ♦XJ*DS*S0BL 
GO  TO  106 

104  SOHO=ZBBOC 
SOaL=ZEROC 
DO  105  J=1,N 
SOHO*SOHO*n(X,J)*»1  (J) 

105  SOHL=SDBL*TL(I,J)*ETO(J) 

*SUB0*IJ*h1  (I, I)  ♦DC*SUHL*XJ*OS*ETI,(r) 

106  CONTINUE 

SOLYE  THE  EQUATIONS 

107  CALL  LEQT1C(A,N,H,B, 1,N,0,8A,IB8) 

IBB=IBa-128 

B8ITE(6,108)  P 

108  FOBBAT(1H1,'  PBEQU ENr-  (h ertZ ) <-  'IPBIO.U  ///> 

IFdEF.SE.  1)  GO  TO  110  ■sio.v,/,/) 

HRITE(6, 109) 

109  F08HAT (///,•  ♦♦♦♦♦SOLUTION  BFROB^^^^+*  ///) 

GO  TO  121  ’ ' 

110  »RITE(6,111) 

111  FORHAT(16X,'8IBE',8I,'IOHlAHPS) •,4X,'IOA(OEGSEES)  • 8X 

1 ItHIAHPS) »,4X,*ILA(DEGREES) •///}  ' ' 


COHPOTB  AND  PRINT  THE  TERHINiL  CURRENTS 


DO  im  1 = 1, N 


HIRE0550 

8IBE0551 

8IBE0552 

WIBE0553 

8IBE0554 

WIRE0555 

8IRB0556 

8IRE0557 

8IBE0558 

8IBE0559 

BIBE0560 

H1BE0561 

BIBE0562 

8IRE0563 

8IRE0564 

8IRE0565 

BIRE0566 

HIRE0567 

HIRE0568 

HIRB0569 

HIRE0570 

8IRE0571 

8IfiE0572 

BIBE0573 

HIRE0S74 

BIRR0575 

HIRB0576 

HIRE0577 

BIRE0578 

8IRB0579 

HIRE0580 

8IBB0581 

81880582 

8IBB0583 

8IRE0584 

HIRB0585 

BIRE0586 

Hir!E0587 

8IRE0588 

8IRE0589 

8IRB0590 

HIRR0591 

8IBK0592 

HIRE0591 

HIRE0594 

HrRE0595 

8IRB0S96 

81880597 

8IRE0598 

WIRE0599 

HIRE0600 

HIRE0601 

8IRE0602 

HtRE0603 

HIRE0604 

8IRE0605 

HIRB0606 

WIRE0607 

HIRE0608 

WIRE0609 

8IRE0610 
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IP(LSO.B0. 11.OB.LSO.B0.21)  GOTO  113 

WIRE0611 

SUHO=ZBROC 

H1BE0612 

DO  112  J:>1,K 

VIBB0613 

112 

S0f10*SUB0*r0  (I,  J)  *B(J) 

ilIBB0614 

Bi(i)=suno 

VIBB061S 

GO  TO  114 

BIBE0616 

113 

Bi(I)=rO(I,I)*B(I) 

HIBE0617 

114 

COBTIHUE 

BIRB0618 

DO  120  1=1, K 

VIBE0619 

IP(LSO.E0.21.0B.LSO. EQ.22)  GOTO  116 

«IBE0620 

I0=B{I) 

i(IBK0621 

SUnO=ZBBOC 

VIRB0622 

DO  115  J = 1,II 

BIBE0623 

115 

suno=saHO»n3  (i, J) *va (J) 

VIRB062* 

lL=-XJ*(m  (I,I)»DS*ETL(I)) ♦DC*B(I)  ♦XJ*DS*SOnO 

HIBB0625 

GO  TO  118 

UIRB0626 

116 

tO=HA(I) 

B1RB0627 

SURO=ZBBOC 

BIRE0628 

DO  117  J=1,N 

HIRB0629 

117 

S0tl0  = S0B0*B3(I,  J)  *B(J) 

UZRB0630 

IL=-XJ*(B1 <I,I)»OS*ETL(I))  ♦DC*«A(I)»*J*DS*SnBO 

BIRE0631 

118 

I0B=CDABS(I0) 

IHBE0632 

ILB=COABS(IL) 

HIRB0633 

IOR=OREAL(IO) 

BIBE0634 

I0I=0IBAG(I0) 

HIRB0635 

ILR=DREAL(IL) 

iiZRE0636 

ILI=DIRAG(Il) 

WIBB0637 

IP(IOR.EO.ZBRO.AMO.IOI.EQ.ZERO)  tOR=OIIE 

liIRE0638 

IF(ILR.EQ.ZERO. ANO.ILI.EQ. ZERO)  ILR=ONE 

BIRE0639 

I0A=DATAN2 (I 01 , lOR) *RADEG 

HIRB0640 

ILA  = 0ATAN2 (I  LI , ILR) • BADEG 

HIRB0641 

HRITE(6,119)  I,I0B,I0A,ILR,ILA 

«IBE0642 

119 

rOBRAT(18X,l2,7x, 1PE10. 3,4X, 1PB10. 3, 9X, 1 PE  10. 3 ,4 1, 1 PB10. 3/) 

VIRB0643 

120 

CONTINUE 

HIRE0644 

GO  TO  66 

i)IBE0645 

121 

STOP 

RIBE0646 

END 

VIRE0647 
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appendix  B-1 


Delete  0062 
Card 


^version  of  WIBE  to  Sl„»l»  I-,.,.,.... 


Double  Precision 


0071-0072  Change  all 


REAL  *8 
COMPLEX  *16 


Single  Precision 


COMPLEX 


DCMPLX(1.D0,0.D0) 

DCMPLX(0.D0,0.D0) 

DCiPLXCO.DO.l.DO) 

DATAN 

DSQRT 


CMPLX(l.i:0,lj.L0) 

CMPLX(O.EO.O.EO) 

CMPLX(0.E0,1.E0) 


DSQRT 
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APPENDIX  B-1  (continued) 


Card 

Double  Precision 

Singl 

0314 

DCOS 

COS 

0315 

DSIN 

SIN 

0316 

DSIN 

SIN 

0317 

DSIN 

SIN 

0332 

DSIN 

SIN 

0333 

DCOS 

COS 

0345 

CDEXP 

CEXP 

0348 

CDEXP 

CEXP 

0349 

CDEXP 

CEXP 

0361 

CDEXP 

CEXP 

0362 

CDEXP 

CEXP 

0363 

CDEXP 

CEXP 

0378 

DSIN 

SIN 

0379 

CDEXP 

CEXP 

0390 

CDEXP 

CEXP 

0391 

CDEXP 

CEXP 

0427 

CEDXP 

CEXP 

0470 

CDEXP 

CEXP 

0505 

CDEXP 

CEXP 

0536 

CDEXP 

CEXP 

0632 

CDABS 

CABS 

0633 

CDAB3 

CABS 

0634 

DREAL 

REAL 

Precision 
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APPENDIX  B-1  (continued) 


Card 

0635 

0636 

0637 

0640 

0641 


Double  Precision 

Single 

DIMAG 

AIi-fAG 

DREAL 

REAL 

DIMAG 

AIMAG 

DATAN2 

ATAN2 

DATAN2 

ATAN2 

Precision 
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ZEROC 


= O.+jO. 
= O.+jl. 


V (velocity  of  light  in 
free  space) 


RADEu (conversion  of  radians 
to  degrees) 


2.997925  x 10 
m/ sec 


= 180. /tt 


Read  and  print; 

1 

Structure  type  (1,2,3) 

= TYPE 

Load  Structure  option 
(11,12,21,22) 

= LSO 

1 

Field  Specification 

option  (1,2) 

= FSO 

1 

Number  of  wires  (n) 

N 

1 

Relative  permittivity  of 

medium  (e  ) 
r 

= ER 

1 

i 

Relative  permeability  of 

medium  (u  ) 
r 

= MUR 

Line  length  (j^) 

= L 

j 

TYPE  = It, 

1 

t 

Radius  of  reference  wire, 

^wO’ 

= RWO 

i 

i 

I 

1 

TYPE  = 3: 

Interior  radius  of  cylindrical 

1 

shield 

= RS 

Read  and  print  entries  in  terminal 
impedance  (admittance)  matrices  at 
x=0,  Zq(Yq),  and  x=/,  Store 

Zp(Yp)  in  array  YO  and  Z^(Yj^)  in  array 
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Compute  the  contributions  to  and 
N due  to  the  longitudinal  electric 
field  along  the  reference  wire. 


(3 


Compute  the  contributions  to  M and 
N due  to  the  longitudinal  electric 
field  along  the  n wires.  Store  M and 
N io  arrays  VI  and  V2,  respectively. 


1 


(inc) 

Compute  the  entries  in  E^(0)  due  to 
the  transverse  electric  field  at  x=0. 
Store  in  array  ETO. 


(inc) 

Compute  the  entries  in  E^(:^)  due  to  the 
transverse  electric  field  at  x=^.  Store 
in  array  ETL. 


Fomi  the  equations  (3-2a)  or  (3-4a). 


Arrays  Vi,  V2,  ETO,  ETL  now  contain 
(inc)  (inc)  , (inc) 

M-Et(2),  N,  E^(0)  and  £^.(0), 

respectively.  Array  Ml  contains  on  its 

main  diagonal  the  entries  in  Z ^ N. 

•M  ^ 


Solve  the  equations  with  LEQTIC 


Compute  the  terminal  currents  at  x=0 
and  those  at  x=^ via  (3-2b)  or  (3-4b). 


vT 

% 

I 

I 

‘i 


Compute  and  print  the  magnitude  and 
phase  of  the  terminal  currents. 


APPENDIX  C 


Function  Subprograms 
El,  E2 

Program  Listings 


COMPLEX  PONCTIOM  B1*16(A.B,X) 

IMPLICIT  REAL*8(A-H,0-Z) 

COMPLEX* 16  XJ.OMEC 

DATA  THREE/3.DO/,PZ1/1.D-2/, TEN/10. DO/ 

COMMON  XJ,ZER0,T«0,0NE,0MEC 

B2=D*B 

A2=A*A 

XB=X*B 

XA=X*A 

BPA»B+A 

XBPA2=X*BPA/TB0 

BHA»B-A 

XBHA2=X*BHA/TKO 
IP(XBPA2.E0< ZERO)  GO  TO  1 
SBPA*DSIR(XBPA2) /XBPA2 
GO  TO  2 

1 SBPA*ONE 

2 IP(XBHA2.E<].ZERO)  GO  TO  3 
SBBA^DSIH (XBnA2) /XBHA2 

GO  TO  4 

3 SBMA^ONE 

4 IP  (XB. BQ.ZERO)  GO  TO  5 
SB30SIN(XB)/XB 

GO  TO  6 

5 SBi’ONB 

6 IP(XA.  EQ.ZBRO)  GO  TO  7 
SA«DSIN(XA)/XA 

GO  TO  8 

7 SA*OHE 

8 XR»-SBPA*SBHA*(B2-A2)/T«0*B2*S»-A2*SA 
IP  (X. BQ.ZERO)  GO  TO  13 
IP(0A8S(XA).LE.PZ1}  GO  TO  9 
XIA*A*(DCOS(XA}-SA)/X 

GO  TO  10 

9 XI  A«-X  A*  A2*  ( (ONE-X  A*  X A/T  EN)  /THREE) 

10  IP(OABS(XB).LE. PZ1)  GO  TO  11 
XIB»B* (DCOS(XB)-SB)/X 

GO  TO  12 

1 1 XIB*-XB*B2*( (ONE-XB*XB/TEN)  /THREE) 

12  XI^XIA-XIB 
GO  TO  14 

13  XI*ZERO 

14  E1-XH+XJ*XI 
RETURN 

END 


PNE10001 

PNB10002 

PHE10003 

PNB10004 

PHE10005 

PNE10006 

PMB10007 

PHE10008 

PMB10009 

PHEtOOlO 

PNE10011 

PNE10012 

PMS10013 

PNE10014 

PHE10015 

PRB10016 

PNE10017 

PHE10018 

PNE10019 

PIE10020 

PNE10021 

PNE10022 

PNE10023 

PNB10024 

PNB10025 

PRBI0026 

PNB 10027 

PNB10028 

PNE10029 

PNE10030 

PNE10031 

PNB10032 

PNB10033 

PNE10034 

FHE10035 

PRE10036 

PMB10037 

PNE1003B 

PNB10039 

PNE10040 

PNE10041 

PNE10042 

PNB10043 

PKB10.044 

FME1004S 
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APPENDIX  C-1 


Conversion  of  El  to  Single  Precision 


' Delete  Card  0002 


001 

El*16 

El 

0003 

C0MPLEX*16 

COMPLEX 

Of'04 

Change  all 

D's 

to 

E's 

0015 

DSIN 

SIN 

0019 

DSIN 

SIN 

0023 

DSIN 

SIN 

0027 

DSIN 

SIN 

0032 

DABS 

ABS 

0033 

DCOS 

COS 

0036 

DABS 

ABS 

0037 

DCOS 

COS 
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COMPLEX  rOMCTIOM  E2*16<A,B,X)  PHE20001 

IMPLICIT  RBIL*a  (i-H,0-Z)  FHE20002 

COHPLEX*16  XJ.OREC  PME20003 

COMHOH  XJ,ZBRO,TVO,ORE,ONBC  PBB20004 

DIP»B-R  FBE20005 

FA«X*DIF/T«0  F(IB20006 

FB«X*(B»A)/T«0  FMB20007 

IF(FA.BQ.ZBRO)  CO  TO  1 FRB2000S 

B2*DIF* (OSIR (FA) /FA) ♦CDEXP  (X J*FB)  FME20009 

GO  TO  2 FHB20010 

1 B.'«DIF*OMEC  FltE20011 

2 ClITIMOE  FRB20012 

RiTORM  FRB20013 

BRO  FKB20014 
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APPENDIX  C-2 


Conversion  of  E2  to  Single  Precision 


Delete  Card  0002 


0001 

E2*16 

E2 

0003 

COMPLEX* 16 

COMPLEX 

0009 

DSIN 

SIN 

0009 

CDEXP 

CEXP 

MISSION 

of 

Rotne  Air  Development  Center 
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RADC  plans  and  conducts  research,  exploratory  and  advanr'od 
developtnent  programs  in  command,  control,  and  comatmications 
(C^)  activities,  and  in  the  cJ  areas  of  information  sciences 
and  intelligence.  The  principal  technical  mission  areas 
are  conenunications , electroms jnetic  guidance  and  control, 
surveillance  of  ground  and  aerospace  objects,  intelligence 
data  collection  arid  handling,  information  system  technology, 
ionospneric  propagation,  solid  state  sciences,  microwave 
physics  and  electronic  reliability,  raaintainahility  and 
compatibility . 


